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Of the twenty-seven internal degrees of freedom of 
propane, all nondegenerate, twenty-two may appear as 
fundamental absorption bands. These bands fall into three 
symmetry classes, designated A,, B, and Be, and dis- 
tinguishable by their characteristic contours. Because of 
overlapping, however, it is impossible in many cases to 
determine their positions precisely. This is especially true 
in the regions of the C—H valence and deformation fre- 
quencies. Some ten or twelve fundamental bands may be 
identified with confidence as well as a number of com- 
binations. An A, band at 870 cm and a B, band at 748 


cm! have been partially resolved, the line spacing béing 
about 1.47 cm™ in agreement with predictions based upon 
electron diffraction measurements. The fine structure of 
the B, bands has not been observed (the predicted spacing 
is 0.5 cm~) but the interval between maxima of the P and 
R branches is approximately 26 cm as expected. With 
24 cm-atmospheres of gas no bands were observed between 
15u and 35y, although the symmetrical C—C deformation 
might be expected to produce a band of appreciable in- 
tensity within these limits. This frequency has apparently 
been observed in Raman spectra at 375 cm“. 





OME years ago the infra-red absorption 
spectrum of propane was investigated with a 
prism spectrometer by Bartholomé.'! He found 
eleven bands between 14y and 3.2yu, and eleven 
others, generally weaker, at wave-lengths below 
3.2u. They were not resolved sufficiently, how- 
ever, to determine the characters of the associ- 
ated motions. Several measurements of the 
Raman spectrum have also been reported,? and 
the distances between carbon atoms have been 
determined by means of electron diffraction.* 
This molecule belongs to the symmetry group 
C,, with the three carbon atoms forming an 
isosceles triangle whose vertex angle is obtuse. 
Twenty-seven internal degrees of freedom are 
indicated, without degeneracies. All the funda- 
* Now at Hwa Nan College, Nanping, Fuchien, China. 
(1933) Bartholomé, Zeits. f. physik. Chemie B23, 152 
2P. Daure, Ann. de physique 12, 375 (1929); S. 
Bhagavantam, Ind. J. Phys. 6, 595 (1932); K. W. F. 
a and F. Képpl, Zeits. f. physik. Chemie B26, 209 


*L. Pauling and L. O. Brockway, J. Am. Chem. Soc. 59, 
1223 (1937). 


mental frequencies are permitted in Raman 
scattering, and twenty-two of the twenty-seven 
should be infra-red active. 

Because of the relatively large moments of 
inertia, complete resolution of the rotational 
structure of the bands cannot be expected. 
Certain ones may be partially resolved, how¥CCr, 
and for the others band envelopes m Vyield 
valuable information concerning the iirections 
of oscillation of the electric moment, as has been 
shown by Dennison,‘ by Nielsen,’ and by Badger 
and Zumwalt.* The last named authors have 
described the three different shapes of band 
envelopes to be expected for oscillations in differ- 
ent directions within the molecule. For propane 

=—0.90 and p=22.7 in their notation.** When 
the direction of oscillation of the electric moment 


4D. M. Dennison, Rev. Mod. Phys. 3, 280 (1931). 

5H. H. Nielsen, Phys. Rev. 38, 1432 (1931). 

®R. M. Badger and L. R. Zumwalt, J. Chem. Phys. 6, 
711 (1938). 

** S=(2b—a—c)/(a—c) and p=(a—c)/b, where a, b, and 
c - reciprocals of the principal moments of inertia, with 
a>b>c. 
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Fic. 1. The pair of Bz bands 
near 13y, and the dcy bands 
near 7p. 
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«. Ci, 2. Contours of several fundamental bands including the partially resolved A, band at 870 cm™ and the 
Kox B, combination band at 864 cm™ (11.5). 
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ABSORPTION SPECTRUM OF PROPANE 


is along the twofold axis (the axis of the inter- 
mediate moment) the symmetrical motion falls 
in the class A;, which includes nine modes of 
vibration. These bands are without obvious Q 
branches. The seven different motions in sym- 
metry class B, are parallel to the plane of the 
carbon atoms but antisymmetric with respect to 
the twofold axis, hence the changes in electric 
moment are along the axis of minimum moment 
of inertia. The resulting bands have P, Q and R 
branches of comparable intensities, the Q branch 
being rather broad. The symmetry class Bz» 
includes six oscillations in which the changes of 
electric moment are along the axis of maximum 
moment of inertia, i.e., normal to the plane of 
the three carbon atoms. They give rise to bands 
with narrow and intense Q branches. A fourth 
symmetry class A» includes five symmetrical 
vibrations, none of which yield infra-red bands 
since they involve no change in electric moment. 

Another rough but convenient classification of 
bands places them in three groups depending 
upon which particular atoms are principally in- 
volved in the motion. Only limiting cases are 
considered, and it must be remembered that the 
actual oscillations are not so simple; in general 
they involve linear combinations of the limiting 
cases. The first group includes motions primarily 
associated with the CHg, the second includes the 
CH; oscillations, and the third those vibrations 
of the triangular structure formed by the three 
radicals each considered as rigid. This third 
group consists of nine torsional oscillations, 
three for each radical, and the three vibrations 
analogous to those of a symmetrical nonlinear 
triatomic molecule XY». The masses of the 
radicals being 15, 14 and 15, these last three 
frequencies should all lie below 1200 cm. 
Dadieu and Kohlrausch’ have predicted that the 
lowest one, corresponding to the symmetrical 
deformation, should occur at approximately 400 
cm~!. Subsequently Bhagavantam observed a 
Raman displacement at 377 cm=', and Kohl- 
rausch and Képpl reported 373 cm-!. There is 
in the Raman spectrum one sum band (1278 
cm~') which seems to involve this frequency, but 
the corresponding difference has not been ob- 
served. Bartholomé failed to find any infra-red 


7 A. Dadieu and K. W. F. Kohlrausch, Akad. Wiss. Wien. 
139, 165 (1930). 
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band near 375 cm, and our own results are 
equally negative. It is interesting to note that 
the dcc vibration in propane appears to be so 
nearly inactive, since the corresponding infra-red 
band is very strong in nonlinear triatomic 
molecules. 


EXPERIMENTAL 


We have examined the spectrum of propane 
between 1.354 and 35y with a prism grating 
spectrometer, using an absorption cell 25 cm in 
length. The windows were either of NaCl or of 
KBr. Since the instrument is not completely 
enclosed, atmospheric water and carbon dioxide 
interfere appreciably with the measurements in 
some regions. Three bands have been partially 
resolved, the contours for several others have 
been determined, and a number of new bands 
have been located. No absorption of measurable 
intensity was found between 134 and 35u. 
Figures 1 to 3 present the observed results. The 
effective slit widths and the amount of gas used 
in each case are indicated in Table I, while 
Table II lists the positions of the various bands, 
their apparent symmetry characters, and possible 
assignments. It must be emphasized that a 
completely unique scheme of assignments would 
require more information than is as yet avail- 
able. Moreover, in regions where several bands 
overlap it is not always possible to identify the 
components without ambiguity; hence the pre- 
cise positions must be uncertain in some cases. 
This is particularly true around 1400 cm~! and 
3000 cm-! where the C—H vibrations occur. 
Some assistance in the interpretation of various 
bands is provided by data already available upon 


TABLE I. 
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TABLE II. Classification of the propane bands. 








SYMMETRY CLASSES 
A2 Bi Bz 


ANALOGS 


ETHANE ETHYLENE CYCLOPROPANE 





2980 


1179 


3107, 3069 3090 
3020, 2988 
1343, 1444 1453, 1430 
948, 943 1020 
1015, 1035 
815 





2980, 2965 
2900, 2955 


1461, 1463 

1375,* 1380 

1170,* 831 
275* 





867 
870 
375 


1184, 865 








* Not observed directly. Italics indicate Raman measurements. Subscripts x and o indicate oscillations parallel to or perpendicular to the 


geometric axis of the group involved. 


the analogous motions in ethane, ethylene and 
cyclopropane. These are also included in Table IT. 


Most of the A, frequencies are determined by 
Raman observations only, and for the liquid 
state. Two of the B, class also appear in Raman 
scattering as rather broad lines. The A, infra- 
red bands, lacking Q branches, are difficult to 
locate precisely unless isolated, as at 870 cm. 
The methyl C—H valence vibrations apparently 
are complicated by resonance interaction, so 
that in class A, there are four Raman lines near 
3000 cm-', and in each of the classes B; and By 
at least two infra-red bands. These, of course, 
cannot be assigned individually to specific mo- 
tions. 

It is suggested that the Bz band at 748 cm~! 
arises from a rocking motion of the methyl 
radicals in a direction normal to the C; plane. 
For the corresponding motion parallel to the 
plane, the B, band at 922 cm is proposed. 
Such specific assignments are perhaps hardly 
justified since both motions involve oscillations 
of the methylene also: in fact the pair of B, 
bands at 1152 and 922 cm~, and the pair of A: 
bands at 1179 and 748 cm™, are each combina- 
tions of the two limiting cases to which they 
have been assigned. The two normal frequencies 


in either case may have been quite close together, 
with the observed separation due to interaction. 
The torsional motions 1, of the methyl radicals 
around their own symmetry axes (in and out of 
phase) are of particular interest, but have not 
been identified. Their frequencies are probably 
low, and in consequence the excited states should 
be well populated. The transitions corresponding 
to any of the 7, oscillations may have somewhat 
different frequencies when the initial state is one 
of the 7, excited levels. This would explain the 
broadening of the observed Q branches, par- 
ticularly that of the 922 cm band, which 
exhibits a second intense maximum at 925 cm". 
Another possible interpretation, which, however, 
seems less probable, is to assign the frequencies 
922 and 925 to the two infra-red 7, bands, and 
the frequency 748 to the 7, band of type Bs. 
The doublet type band at 870 cm~! may be 
assigned with confidence to the symmetrical 


TABLE III. Observed values of Avm for By bands. 
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27 
26 
26 
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valence vibration of the C; framework. Appar- 
ently a second band, either of class B; or Bs, is 
superposed upon it, the Q branch appearing at 
864 cm~!. In the neighborhood of 1936 cm! a 
band is observed which, though only partially 
resolved, appears to be of the doublet (A;) 
type, with a very irregular contour. This is 
probably a combination of the two B, bands at 
922 and 1053 cm~!. One other combination, 
observed as a Raman line, has already been 
mentioned; it is 1278 =922(B,)+375(A,). Two 
moderately intense bands remain to be con- 
sidered. They both have prominent Q branches, 
one at 1338 cm and the other at 720 cm”; 
and they appear to belong to type Be. No positive 
assertion regarding their origin can be made. 
Even the lower frequency seems much too high 
for the fundamental 7,. It also seems too intense 
to represent the difference transition 1468-748. 
A possible explanation is that these two bands 
are sum and difference combinations of a low 
frequency Az band (7,) with the B; band at 
1053. This would indicate a frequency of 333 
cm! for one of the 7, levels. 

A number of other infra-red absorption regions 
have been observed, mostly at higher frequencies, 
representing multiple and combination bands. 
These are displayed in Fig. 3. Although some 
resolution is indicated, it seems hardly feasible 
at present to assign band centers, or to make 
definite assertions regarding their interpretation. 

The Bz band at 748.0 cm, shown in Fig. 1, 
has been resolved into a fairly regular set of 
some twenty lines, almost equally spaced, as 
would be expected for a molecule having two 
nearly equal moments of inertia. By isolating 
the normal state a value of 1.475 cm has been 
obtained for h(a—d)/4x*, where a is the re- 
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ciprocal of the smallest moment of inertia, and d 
is the mean of the reciprocals of the other two. 
The A, band at 870 cm has also been partially 
resolved, with an average line interval in the R 
branch of 1.40 cm~!. These frequency differences 
compare favorably with the value 1.47 computed 
for A; and Bz bands using the C—C distances 
obtained by electron diffraction, and the C—H 
distance from methane. The corresponding line 
interval in the B, bands should be about 0.5 


TABLE IV. Positions of observed lines in partially 
resolved bands. 








BAND AT 
1936 cm™! 
P 


BAND AT 
870 cm"! 
P R 


BAND AT 
748 cm"! 
Pp R 





1914.9 
1916.8 
1918.7 
1922.3 
1924.5 
1926.0 
1927.7 
1930.6 


1937.6 
1940.5 
1943.5 
1945.0 
1949.0 
1951.9 
1954.0 


862.4 872.0 
864.0 874.2 
867.0 875.6 
877.1 
878.7 
880.2 
881.6 
883.0 


749.6 
751.0 
753.5 
754.7 
756.3 
757.7 
759.1 
760.3 
761.8 
763.4 


733.1 
734.7 
736.4 
737.8 
739.2 
740.9 
742.2 
743.5 
745.1 








cm~'!, which could be observed only under very 
advantageous circumstances. The interval Av,, 
between maxima of the P and R branches in B, 
bands may be computed approximately, follow- 
ing Gerhard and Dennison,® its value being 26 
cm~!. Table III lists the values of Av,, observed 
for various B, bands. In Table IV the positions 
are given for the observed fine structure lines in 
the bands at 748 cm and 870 cm~, and also 
for a number of the maxima appearing in the 
complex band near 1936 cm-. 


8S. L. Gerhard and D. M. Dennison, Phys. Rev. 43, 197 
(1933). 
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The Structure of Black Carbon 
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This paper is concerned with the structure of extremely small crystallites of carbon deposited 
upon silica by the pyrolysis of methane. Electron diffraction patterns have been obtained from 
more than 100 carbon films, from some by the reflection method, and from others by transmis- 
sion after removing the films from their supports with hydrofluoric acid. In some of the speci- 
mens the individual crystallites are strongly oriented and reflection patterns obtained from 
these are approximations to rotation patterns from single crystals. From these patterns we 
have concluded that the films are made up of pseudo-crystals in each of which the carbon atoms 
are hexagonally arranged in planes as they are in graphite, but successive parallel atomic layers 
are displaced so that no regularities exist other than the uniform separation of the rm and 
the regular arrangement of atoms in each of them. 





HE structure of what was formerly called 

amorphous carbon has been the subject of 
study by x-ray diffraction for many years, the 
first pattern from carbon having been obtained 
by Debye in 1916. Only powder patterns from 
randomly oriented small crystals have been 
studied.! From these it has been inferred that 
the structure of black carbon resembles that of 
graphite with discrepancies which increase with 
decreasing crystal size until, when the mean 
linear dimensions of the crystals are of the order 
of 20A, the pattern consists of only three dis- 
tinguishable rings and resembles that of an 
amorphous material.2 One can see that the 
number of measurable parameters is so limited 
that on the basis of such a pattern alone no de- 
ductions as to the arrangement of atoms in the 
smallest carbon crystallites can be confidently 
made. Nevertheless, structures have been as- 
signed to the crystallites of such carbons by a 
sort of extrapolation from the results obtained 
from other carbons whose structures degenerate 
progressively from that of graphite. 

We have obtained electron diffraction patterns 
from which the structures of such extremely 
small crystals can be reliably inferred without 
reference to any other carbon. In our experiments 
a marked degree of preferential orientation is 
observed, and the diffraction patterns are approxi- 


1X-ray diffraction patterns from oriented carbon 
crystallites were, however, obtained by M. Pirani and W. 
Fehse, Zeits. f. Elektrochemie 29, 168-174 (1923). 

2 The most thorough investigation is that of V. Hoffman 
and D. Wilm, Zeits. f. Elektrochemie 42, 504 (1936). 


mations to those which would be produced by a 
single crystal rotated about a principal axis. 
Such patterns yield evidence which is not avail- 
able in the previously reported patterns from 
unoriented crystallites. 

Representative electron diffraction patterns 
are reproduced in Figs. 1-3. The three films from 
which these patterns were obtained were pre- 
pared by deposition of carbon upon flat fused 
silica surfaces in a furnace containing methane 
and nitrogen* maintained at about 1000°C. The 
pattern of Fig. 1 was produced by electrons 
scattered by reflection from the surface of the 
deposited carbon. The thin films which gave the 
other two patterns were first removed by means 
of hydrofluoric acid from the silica plates on 
which they were formed, and each was then sup- 
ported across a narrow slit in such a way that 
an electron beam could pass directly through it. 
The pattern of Fig. 2 was produced with the 
film normal to the electron beam, and that of 
Fig. 3 with the film rotated from a position 
normal to the beam by an angle of 45° about a 
horizontal axis. 

On the original plate of the pattern of Fig. 1 
one can observe six uniformly spaced diffuse 
spots or arcs lying along the center line, The 
best estimate of the separation between these is 
s=6.4(0) mm, and from this the equivalent 
crystallographic spacing d can be calculated by 


3 The degree of orientation of the crystallites in the 
surface of the film which gave the pattern of Fig. 1 had been 
increased by mild oxidation in air, subsequent to the 
original deposition of the film. 
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Fic. 1. Electron diffraction pat- 
tern obtained by the _ reflection 
method from the surface of a film 
of carbon deposited upon fused 
silica by the pyrolysis of methane. 
(Reproduced 0.6 full size, LA = 2.34 
X 10-8 cm?.) 


sequently 
fluoric acid. 


means of the relation d=LX/s, where L is the 
distance from specimen to photographic plate 
and X is the electron wave-length. Substituting 
for LX the independently determined value 2.34 
X10-*§ cm’, one obtains d=3.6(6)A which is 
about 9 percent larger than ¢o/2=3.35A, the 
separation between adjacent hexagonal planes in 
graphite. From the result of this comparison it is 
natural to infer that the pattern of Fig. 1 was 
produced by graphite-like crystallites oriented 
with their hexagonal planes predominantly par- 
allel to the supporting quartz surface. If this 
inference is correct there must be features in 
Fig. 1 along a horizontal line through the primary 
beam position, insofar as the shadow of the 
specimen permits, which correspond to reflections 
having Miller indices of the form (hkO), and 
perhaps also features at other positions corre- 
sponding to general indices (hkl). Crystallo- 
graphic planes of the hexagonal system for which 
the third Miller index is zero make up a set 
having spacings inversely proportional to the 
square roots of the whole numbers 1, 3, 4, 7, 9, 
12,13--- ete., and the radii of electron diffraction 
rings corresponding to these separations must be 
directly proportional to these square roots. 
Actually one can observe in Fig. 1 five features 
lying approximately along the horizontal line 
through the primary beam position, and the 
radii of these, as closely as they can be measured, 
are proportional to the square roots of the 
integers 1, 3, 4, 7 and 9. The calculated edge of 
the hexagonal cell is a9 =bo9=2.4(6)A, which is 
just the corresponding value for graphite. There 
are, however, no traces whatsoever of features 


OF 


Fic. 2. Pattern by transmission 
at normal incidence from a carbon 
film deposited upon silica and sub- 
removed with hydro- 


BLACK CARBON 


Fic. 3. Pattern from a carbon film 
inclined 45° to the primary beam, 
by rotation about a horizontal axis. 


other than the (00/1) and (hkO) reflections. The 
diffraction pattern thus gives no evidence of 
regularities in the structure of the individual 
crystallites in this specimen of black carbon other 
than in the arrangement of carbon atoms in the 
hexagonal planes and in the separation between 
these planes. The relative intensities of the 
different (hkO) features correspond furthermore, 
as closely as we can estimate, to intensities to 
be expected from parallel but otherwise un- 
related planes of carbon atoms, and not to in- 
tensities from crystallites having the complete 
graphite structure; relative to a corresponding 
pattern from oriented graphite crystals the 
“arcs” corresponding to the integers 3 and 9 are 
much too weak, although, of course, much 
stronger than similar features which would be 
produced by parallel but otherwise unrelated 
planes of atoms having a close packed arrange- 
ment. The separations and relative intensities of 
the (hkO) features would have enabled us to 
deduce directly the arrangement and spacing of 
atoms in the separate planes, quite independently 
of knowledge of the graphite structure. 

These deductions in regard to the structure of 
the carbon film are borne out by a more detailed 
consideration of the diffraction pattern. One can 
deduce very simply the pattern to be expected 
from an array of parallel and uniformly spaced 
but otherwise completely independent planes of 
carbon atoms. The lattice reciprocal to such a 
set of cross gratings is made up of the (00/) 
spots of the graphite reciprocal lattice, and the 


surfaces of cylinders normal to the hexagonal 
planes and generated by (hk) lines which pass 
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(hk) 











/ Ls x 
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Fic. 4. Diagram illustrating the production of an (hk) 
diffraction feature of Fig. 1 by preferentially oriented 
cross gratings of carbon atoms. 


through the points (hk0) of the graphite re- 
ciprocal lattice. For the case of hexagonal planes 
parallel to the surface and to the primary beam, 
the sphere of reflection, which is approximately a 
plane, intersects these cylinders in straight lines 
normal to the shadow edge. If the various pseudo- 
crystallites are imperfectly aligned, with the 
hexagonal planes of many of them making small 
angles with the general surface, then the resulting 
diffraction lines will be distributed over an 
angular range about the normal to the shadow 
edge. The superposition of many such lines will 
be expected to give to each of the (hk) diffraction 
features the shape of a fountain, or narrow fan, 
appearing to diverge approximately from its 
intersection with the horizontal line through the 
primary beam position; one observes just this 
shape for those (k0) features of Fig. 1 whose 
form can be clearly made out. 

It is not difficult to calculate roughly the 
distribution of scattering intensity to be expected 
in one of the (hk) features of the pattern of 
Fig. 1. This distribution is quite approximately 
the same as it would be if the carbon planes were 
completely unrelated, and the calculation is 
carried out on that basis. Let x represent the 
angle between the general plane of the surface 
and the planes of a set of parallel cross gratings, 
and assume that the number of gratings having 
normals lying in the element of solid angle dw is 
proportional to f(x)dw. In Fig. 4 the plane of 
the paper is chosen to represent the sphere of 
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reflection and the horizontal line through the 
origin its intersection with the plane of the 
specimen surface. The lattice reciprocal to each 
cross grating is the set of (hk) lines described 
earlier, and the reciprocal lattices of many planes 
parallel to each other form coaxial (hk) cylinders. 
The reciprocal lattice cylinders of those sets of 
parallel planes, which are parallel to the primary 
beam, will intersect the plane of Fig. 4 in straight 
lines. The broken reciprocal lattice lines (hk) 
correspond to the two orientations x for which 
sets of cross gratings parallel to the primary 
beam contribute to the intensity at the point 
R, 6. The total intensity at R, 6 is made up of in- 
crements from all sets of cross gratings whose 
normals are parallel to the tangents drawn from 
R, 6 to the sphere of radius Ro = LX/d, where d is 
the spacing between the cross-grating lines which 
correspond to the feature under consideration. 
By assuming for the distribution function 
f(x) =exp (b cos? x) and integrating over all 
values of x corresponding to gratings whose 
normals are parallel to these tangents, one finds 
that the number of gratings contributing to the 
intensity at R, @ is proportional to 


+1 
(1—Ro?/R?) f exp [6B(x+A)? ]dx/(1—x?) ’ 
=" 1) 


where 
A =(R?—R,*)*(tan 0)/Ry and B= Ry?(cos? 6) /R°. 


One represents random orientation of cross 
gratings by setting b=0, and the number of 
gratings contributing to the intensity at R, 6 for 
this case is proportional to (1—R,?/R?)! as one 
sees by evaluating the integral. One can expect 
that the actual intensity at R, 6 for 640 will be 
given by the ratio of these two quantities multi- 
plied into the intensity at R, @ for random orien- 
tation. The latter as calculated by v. Laue? is 
E?/R(R?— R,*)*. Here E, the atomic form factor, 
can be replaced’ by (Z—f)/R? which, in the 
range with which we are concerned, is approxi- 
mately proportional to 1/R. In this way one 
finds that the intensity at R, @ will be approxi- 


4M. v. Laue, Zeits. f. Krist. 82, 133 (1932). This ex- 
pression is infinite at R= Ry, due to the assumed unlimited 
extent of each grating. We understand that Dr. B. E. 
Warren has carried through the v. Laue calculation for the 
case of randomly oriented gratings of limited extent. 

5 N. F. Mott, Proc. Roy. Soc. A127, 658 (1930). 
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mately, J= Expression (1)/R?(R?—R,?). In order 
to determine the value of the constant b we have 
estimated that, for the (10) feature, the intensity 
in Fig. 1 is about one-tenth as great at R= Ro, 
6=30°, as at R=Ro, 06=0°, and by numerical 
integration we have found that )=8 corresponds 
to a value of this intensity ratio of 9.4. Using 
this value (b=8) the intensity has been calcu- 
lated® at various points within the (10) feature, 
and the values obtained are plotted as a contour 
diagram in Fig. 5. (The numbers on the contour 
lines of Fig. 5 are J multiplied by 50.) Except 
for the region close to R= Ro, where the assump- 
- tion of unlimited extension of each grating intro- 
duces serious error, the calculated intensity 
distribution of Fig. 5 is similar to that which is 
apparent in Fig. 1. This seems to us to confirm 
the assumption that the various cross gratings 
of which each crystallite is composed are un- 
related except for their parallelism and uniform 
separation. (One cannot, however, determine 
from the experiments whether or not successive 
atomic layers are merely displaced laterally in 
random fashion from their positions in the 
graphite lattice, or are both displaced and 
rotated.) 

We believe that the considerable vertical 
breadths of the (00/) features in Fig. 1 arise 
from poor resolving power of the individual 
crystallites due to small dimensions parallel to 
their c axes. Some of the breadth of these features 
might, conceivably, be attributed to non-uniform 
separation of atomic planes, but this would 
result in progressive increase in breadth with 
increasing order of reflection. As closely as we 
can estimate the breadths at half-maximum of 
the various orders are about the same, and we 
therefore conclude that the planes in each crys- 
tallite are uniformly spaced and that the spacing 
is the same in different crystallites. From the 
vertical breadths we estimate that the mean 
dimension in the direction parallel to the ¢ axes 
is 8A, indicating that on the average the crys- 
tallites are made up of only three planes of carbon 
atoms. With this fact in mind the striking con- 
trast between the character of the (00/) and the 
(hk) reflections is especially significant, as it 


* We are indebted to Dr. L. A. MacColl for advice and 
assistance in the numerical integrations. 
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offers further confirming evidence that there is 
no regularity in the way the uniformly spaced 
planes of atoms lie upon each other in each 
crystallite, even as regards one plane and its 
next neighbor. 

This conclusion is supported by diffraction 
patterns obtained by the transmission method. 
Measurements of the features of the patterns of 
Figs. 2 and 3 indicate that all of them, as in 
Fig. 1, have indices of the form (#0) or (00/). 

The rings in the pattern obtained at normal 
incidence (Fig. 2) have radii proportional to the 
square roots of the whole numbers 1, 3, 4, 7, 9 
and 12 and thus correspond to cross-grating 
reflections of the form (hk). There are no traces 
of (00/) reflections, and we know therefore that 
the cross gratings in this film are strongly 
oriented, as is the case in the film which produced 
the pattern of Fig. 1. Microphotometer curves 
obtained from the (10), (11) and (20) rings of 
Fig. 2 are reproduced as Fig. 6, the ordinates 
being the differences between the logarithms of 
the light intensity J transmitted through the 
photographic plate and the logarithms of the 
intensity Jg corresponding to the estimated back- 
ground curve. These microphotometer curves are 
more satisfactory than similar curves previously 
obtained from x-ray diffraction patterns of black 
carbon, because strong (00/) diffraction rings are 
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Fic. 5. Calculated intensity distribution in a (10) diffraction 
feature of Fig. 1, plotted as a contour diagram. 
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Fic. 6. Microphotometer curves of the (10), (11) and (20) 
diffraction rings of the pattern of Fig. 2. 


present in the latter and overlap the (hk) re- 
flections. The maxima of Fig. 6 are asymmetrical 
in the manner to be expected from imperfectly 
oriented cross gratings,‘ and we cannot calculate 
the mean breadth of the cross gratings directly 
from the widths of the curves at half-maximum. 
A fair estimate can, however, be obtained by 
substituting, for the ring broadening in the 
Scherrer formula for crystal size, twice the 
difference between the radius of a particular 
ring at maximum intensity and that at half- 
maximum on the inner side. This will not yield 
exactly the mean linear dimension but the error 
will not be important. Estimated in this way, the 
widths at half-maximum, AR, of the (10) and (11) 
peaks of Fig. 6 are 2(0.54) mm and 2(0.65) mm. 
These correspond to average linear dimensions, 
e, of the cross gratings of 22A and 18A, respec- 
tively, as calculated from the Scherrer formula in 
the form e=Z\/AR. The reason for the dis- 
crepancy between these estimates is not clear; 
similar disagreement has been observed in x-ray 
determinations.? Microphotometer curves were 
taken upon another pattern through the (10), 
(11), (21) and (30) rings, and an attempt was 
made to account in some simple fashion for the 
variation in half-breadth which was observed. 
In this we have been unsuccessful. 

The diffraction features which appear on the 
pattern of Fig. 3 are (002), (004) and the (hk) 
reflections corresponding to the whole numbers 
1, 3, 4, 7, 9 and 12. The (00/) reflections are 
those which appear most strongly along the 
vertical center line, and the (hk) reflections along 
the horizontal. The fact that even along the 
horizontal line the (002) reflection appears 
weakly, whereas this reflection is entirely absent 
in the pattern of Fig. 2, is evidence that the 
crystallites in the film which produced the 


pattern of Fig. 3 were less sharply oriented than 
those in the film which gave the latter pattern. 
Except for this quantitative difference the orien- 
tation was, however, the same in the two films. 
The difference between the character of the 
(00/) and the (hk) reflections is brought out 
almost as clearly in Fig. 3 as in Fig. 1. One 
observes in Fig. 3 that the locus of (002) re- 
flections is a circle, whereas the (10) and the (11) 
features are not true circular arcs but flare out- 
ward above and below the center line. The 
contour diagram of Fig. 5 can be used to estimate 
very crudely the intensity distribution to be 
expected within these features. To accomplish 
this one imagines this diagram reflected through 
theline @=0° and both parts then rotated through 
360° about @=90° as an axis, resulting in an 
appropriate volume density distribution in re- 
ciprocal space. Intensities in the pattern of 
Fig. 3 are then represented by the intersection 
of this solid with a plane normal to the plane of 
the paper and inclined 45° to 6=0. The hori- 
zontal radial distribution through an (hk) feature 
of Fig. 3 will be represented by the contour lines 
at 6=0 in Fig. 5, and the vertical radial distribu- 
tion approximately by the contour lines at @=45° 
in the same diagram. One should realize, of 
course, that the diagram of Fig. 5 was calcu- 
lated by assuming cross gratings of infinite ex- 
tent, whereas the gratings are actually only about 
20A wide; furthermore, the diagram corresponds 
to more strongly oriented crystallites than was 
the case in the film which gave the pattern of 
Fig. 3. 

It has been pointed out that the degree of 
preferential orientation is different in different 
films (for example, those which produced the 
patterns of Figs. 2 and 3). More than a hundred 
films have been investigated and it has been 
discovered that the degree of orientation is 
considerable in some and slight in others. The 
patterns chosen for consideration here have been 
selected from those showing a high degree of 
orientation. 


SUMMARY 


We have shown that these films of pyrolyti- 
cally deposited carbon consist of pseudo-crystals 
each made up of 3 or 4 parallel and uniformly 
spaced but otherwise randomly disposed atomic 
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SPECTRA OF SIMPLE PORPHYRINS 


planes in each of which the carbon atoms are 
hexagonally arranged as in graphite. This struc- 
ture is just that inferred from x-ray studies’ of 


7For example, V. Hoffman and D. Wilm, reference 2; 
H. L. Riley, Sci. J. Roy. Coll. Sci. 10, 10-20 (1940). 
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similar carbons, but in our experiments evidence 
is obtained which is not available in the pre- 
viously reported patterns from unoriented crys- 
tallites and the deductions are made without 
reference to any other carbon. 
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Spectroscopic Studies of the Simpler Porphyrins 


III. The Absorption Spectra of ms-Tetraphenylporphine and a Series of Its 
Metal Complex Salts 


V. M. ALBERS AND H. V. KNorR 
The C. F. Kettering Foundation for the Study of Chlorophyll and Photosynthesis, Antioch College, 
Yellow Springs, Ohio 
(Received April 14, 1941) 


The molecular absorption coefficients of ms-tetraphenylporphine and the Cu, Ags, Mg, 
Zn, Cd, Hg, SnClz, Pb, MnCl, FeCl, Co and Ni complex salts have been measured over the 
visible region of the spectrum. The curves representing the absorption coefficients as a function 
of the wave-length are shown. The variations in the spectra of the metal salts of this one 
ms-tetrasubstituted porphine are much greater than those which have been found in the spectra 
of the different ms-tetrasubstituted porphines which so far have been studied. 


INTRODUCTION 


thon two previous publications in this series 


have been concerned with absorption 
spectra of synthetic porphines which were not 
metal complex salts. The first! was a study of 
the spectra of porphine and four ms-tetra- 
substituted porphines, of which ms-tetraphenyl- 
porphine was one. The second? was a study of 
the spectra of the two isomeric forms of those 
porphines for which the two isomers were 
available. At the time of the first investigation 
the ms-tetraphenylporphine was the most diffi- 
cult of the series to prepare and not enough of 
the substance was available to perform a 
separation of the isomers, consequently the 
material used in that investigation of the 
spectrum of ms-tetraphenylporphine contained a 
mixture of the two isomers. Since that time, a 
procedure for the synthesis of this porphine has 


1V. M. Albers and H. V. Knorr, J. Chem. Phys. 4, 422- 
425 (1936). 

*H. V. Knorr and V. M. Albers, J. Chem. Phys. 9, 197- 
203 (1941). 


been discovered* which yields quantitatively 
only the isomer of HCI number 13.5.* This is a 


Fic. 1. Structural formula of ms-tetraphenylporphine. 


3 Rothemund and Menotti, J. Am. Chem. Soc. 63, 267- 
270 (1941). 

* The HCl number of a porphyrin is defined as the 
concentration of HCI which will extract two-thirds of the 
porphyrin from an equal volume of ether solution on 
thorough shaking. 
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Fic. 2. Absorption spectrum of ms-tetraphenylporphine in 
benzene solution. 
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Fic. 3. Absorption spectrum of the ms-tetraphenylporphine 
Cu complex salt in benzene solution. 


very important discovery, since the separation 
of the isomers of these substances involves a 
tedious process of HCI fractionation from ether 
solutions. Accordingly, this porphine was chosen 
as one from which’a long series of metal complex 
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Fic. 4. Absorption spectrum of the ms-tetraphenylporphine 
Ag» complex salt in benzene solution. 


salts was prepared! in order to study the variations 
of the properties with the substituted metal. 

All of the preparations were purified by 
recrystallization and their composition verified 
by analysis. In the case of the magnesium 
complex salt there is still some question about 
the structure; the analyses were consistent with 
the structure in which the magnesium goes into 
the molecule as MgOH but this structure has 
not yet been definitely proved.‘ 

The formula for ms-tetraphenylporphine is 
shown in Fig. 1. The metal complex salts are 
prepared by substituting the metal for the 
hydrogens which are attached to the nitrogen 
atoms in the pyrrole rings. When the metal is 
monovalent, as is the case with the silver salt, 
two metal atoms replace the two hydrogen atoms. 
When the metal is divalent the two valences of 
the metal are linked to the two nitrogen atoms 
and, in the case of trivalent and tetravalent 
metals, the extra valences are usually satisfied 
with chlorine atoms. 


4 Rothemund and Menotti, J. Am. Chem. Soc. In proces: 
of{publication. 
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The spectra were measured with a recording 
spectroradiometer which has been previously 
described.’ The spectra were examined, in all 
cases, from 900 mu to 370 mu. Since none of the 
substances showed characteristic absorption 
above 700 my, only the region between 700 mu 
and 370 my was actually measured. 

The spectra were all measured for benzene 
solutions and the absorption coefficients were 
calculated from Beer’s law in the form of Eq. (1), 


I7/Ib=e-“?, (1) 


where I7/Io is the ratio of transmission of the 
absorption cell filled with solution to its trans- 
mission when filled with the pure solvent, c is 
the concentration of the solution in moles per 
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Fic. 5. Absorption spectrum of the ms-tetraphenylporphine 
Mg complex salt in benzene solution. 


°V. M. Albers and H. V. Knorr, J. Opt. Soc. Am. 28, 
121-123 (1938). 
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Fic. 6. Absorption spectrum of the ms-tetraphenylporphine 
Zn complex salt in benzene solution. 


10-5 


cients 


Molecular Absorption Coefficients x 1072 


Moiécuiar Absorption C 





Wavelength in Mu 


Fic. 7. Absorption spectrum of the ms-tetraphenylporphine 
Cd complex salt in benzene solution. 


liter, x is the thickness of the absorption cell in 
cm and a is the molecular absorption coefficient 
of the substance. 

In measuring these spectra with the spectro- 
radiometer, a 108-watt spiral filament lamp was 
used as a source to illuminate the entrance slit 
of the Hilger Quartz Monochromator. The 
absorption cell was placed between the exit slit 
of the monochromator and a photoelectric cell. 
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Fic. 8. Absorption spectrum of the ms-tetraphenylporphine 
Hg complex salt in benzene solution. 


The photoelectric cell current was amplified with 
a d.c. amplifier whose output was fed to a Leeds 
and Northrup Speedomax recorder. Records 
were made by driving the wave-length drum of 
the monochromator and the record paper in 
the recorder with synchronous motors operating 
from the same a.c. line. In order to obtain 
I7/Iy as a function of wave-length it was neces- 
sary to have a record made with the absorption 
cell filled with pure solvent and another record 
made with the cell filled with the solution. The 
deflections of the recorder pen were then meas- 
ured at 2.5-my intervals and the ratio dr/do 
=I7/Ip, where dr was the deflection when the 
cell was filled with solution and dy was the 
deflection when the cell was filled with solvent, 
since it has been shown that the deflection of 
the recorder pen at any given wave-length is 
proportional to the light intensity on the photo- 
electric cell at that wave-length.’ Since it is 
necessary that the temperature of the lamp 
filament and the gain of the amplifier remain 
constant throughout the making of the records 
for dr and do, the following procedure was used 
for each solution concentration. Four records 
were made, the first and fourth were made with 
the absorption cell filled with the solvent and 
the second and third were made with the cell 
filled with the solution. At the beginning and 
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Fic. 9. Absorption spectrum of the ms-tetraphenylporphine 
SnCl, complex salt in benzene solution. 


end of each record the zero position, i.e., the 
position of the pen when no light was allowed 
to pass through the system, was obtained. By 
this means any variation of the lamp or amplifier 
could be detected and, in addition, two inde- 
pendent determinations for dr and do were 
available at each wave-length, the mean values 
of which were used to compute the values 
of Tr/Io. 

The value of a was determined at 2.5-mu 
intervals between 700 mu and 370 mu. It was 
found that a concentration of from five to seven 
milligrams of material in ten cc of solvent was 
sufficient for the regions of lower absorption. 
A series of concentrations made by diluting a 
given volume of solution with an equal volume 
of solvent was made up and records for each 
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Fic. 10. Absorption spectrum of the ms-tetraphenylporphine 
Pb complex salt in benzene solution. 
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Fic. 11. Absorption spectrum of the ms-tetraphenylporphine 
MnCl complex salt in benzene solution. 
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Fic. 12. Absorption spectrum of the ms-tetraphenylporphine 
FeCl complex salt in benzene solution. 
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Fic. 13. Absorption spectrum of the ms-tetraphenylporphine 
Co complex salt in benzene solution. 
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Fic. 14. Absorption spectrum of the ms-tetraphenylporphine 
Ni complex salt in benzene solution. 


concentration were made until a solution of 
sufficient dilution was obtained to determine 
with maximum accuracy the value of @ in the 
region of greatest absorption. An analysis of the 
propagation of errors in calculating a from the 
values of I7/Iy indicates that only values of 
I7/Io between 0.15 and 0.50 should be used in 
order to prevent the errors involved in reading 
the deflections from producing excessive errors 
in a. Whenever the values of I7/Iy for two 
concentrations of solution fell within this range 
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at any one wave-length, both were used in 
determining the value of a and a mean of the 
two values of a was used. The values of a at the 
tops of the bands and in the regions between 
bands are accurate to about + one or two 
percent. The monochromator slits were operated 
at about 0.001-inch width. Under these condi- 
tions the two yellow lines of mercury are just 
resolved. 

Figs. 2-14 show the spectra of the substances 
studied. In all cases, the absorption in the violet 
end of the spectrum is much greater than that 
in the rest of the region so the curves are broken 
at the point where this heavy absorption sets in 
and the strong band is plotted on a scale 1/20 
as large. The ordinates for the part of the curves 
to the left of the break are indicated at the left 
edge of the graphs and those for the part of the 
curves to the right of the break are indicated 
at the right edge of the graphs. 


DISCUSSION 

In comparing these data with those previously 
obtained from ms-tetrasubstituted porphines, 
one is impressed by the fact that the variations 
in the spectra of the different metal salts of 
ms-tetraphenylporphine are much greater than 
those which have been found in the spectra of the 
different ms-tetrasubstituted porphines. Seven 
of the salts (Cu, Age, Zn, Hg, FeCl, Co and Ni) 
have spectra with two principal bands while 
five of them (MnCl, Pb, SnCle, Cd and Mg) 
have three principal bands. In addition to the 
principal bands all of these substances have one 
or more weak bands. In Table I the wave-lengths 
and absorption coefficients of all of the bands in 
the spectra of these substances between 700 mu 
and 370 mu are tabulated. 

Beer’s law apparently holds throughout the 
visible region for all of these substances, since ab- 
sorption coefficients for different solution concen- 
trations at the same wave-length agree very well. 
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Helium Repulsive Potential from Collision Cross Section Measurements* 


I. AMpuR AND H. PEARLMAN 
Research Laboratory of Physical Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 


(Received April 23, 1941) 


Total collision cross sections have been measured for fast helium atoms, having velocities 
corresponding to 300-1000 electron volts, scattered in room temperature helium gas. Using 
classical theory, attempts were made to find a repulsive intermolecular potential which could 
account for the experimental results. A hard sphere model or a Lennard-Jones type potential 
were both unsatisfactory, whereas exponential type potentials of the form V(r) =A exp (—ar") 
seemed promising. For the distances involved in the scattering experiments, 0.55 to 1.05A, the 
specific repulsive potential which predicted total collision cross sections within the experimental 


error was V(r) =11.3 exp (—4.63r!) X 10-" erg. 





HE elastic scattering of high velocity 

neutral particles by slow moving neutral 
particles affords a method for experimentally 
determining intermolecular repulsive forces. If 
the fast particles have energies of the order of 
100 electron volts or greater, scattering beyond 
about one degree results from interactions over 
very small distances where only repulsive forces 
are of significance. The calculations are further 
simplified by the fact that the scattering par- 
ticles, whose velocity seldom reaches one percent 
of that of the fast particles, may be regarded 
as initially at rest. 

In the present experiments, beams of helium 
atoms with energies from 300 to 1000 electron 
volts were scattered by helium gas at room 
temperature. The apparatus and experimental 
procedure have been previously described! in 
connection with the determination of total 
collision cross sections of high velocity hydrogen 
atoms scattered in room temperature molecular 
hydrogen. As in the earlier experiments, total 
collision cross sections for the helium atoms 
scattered in helium gas were computed from the 
ratio of axis beam intensities after passage 
through helium at two different pressures. The 
average angular aperture of the detecting system 
was, as before, 2.25 degrees. The accuracy of the 
cross sections, computed from 11 individual 
runs, is estimated at 15 percent. The values are 
listed in Table I as a function of beam particle 
voltage. 


* Contribution from the Research Laboratory of Physical 
Chemistry, Massachusetts Institute of Technology, No. 
478. Presented before the American Physical Society at the 
Cambridge Meeting, February 22, 1941. 

‘ Amdur and Pearlman, J. Chem. Phys. 8, 7 (1940). 


TABLE I. Total collision cross sections of helium atoms 
scattered in helium. 
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As pointed out by Massey and Mohr,’ it is 
permissible to use classical theory to describe 
collisions in which the beam particle is scattered 
through a relative angle greater than a critical 
value 90,. 0, is approximately equal to A/2ro 
where \ is the de Broglie wave-length of the 
colliding system and 7», the equivalent hard 
sphere radius, that is (cross section/z)'. The 
use of classical theory is definitely justified in 
the present calculations since the calculated 0, 
values are nearly constant at the various 
voltages, with a value close to 0.5°, whereas all 
measured scattering occurred at relative angles 
greater than 4.5° (2.25° in apparatus coordinates). 

Kennard? has derived the equations applicable 
to the classical scattering by a fixed scatterer of 
a beam of particles homogeneous in velocity, 
and has indicated the changes required when 
the scatterer, although initially at rest, is free to 
move after collision, as in the present experi- 
ments. In this case the equations remain formally 
unchanged if relative coordinates are used 
throughout. The mass of the beam particle is 
replaced by the reduced mass of the system; the 
angular deflection of the beam particle with 
respect to the beam axis, by the relative angle 
of scattering; and the initial beam velocity, by 


2 Massey and Mohr, Proc. Roy. Soc. Al41, 434 (1933). 
3 Kennard, Kinetic Theory of Gases (McGraw-Hill Book 
Company, Inc., 1938), Chapter III. 
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the relative initial velocity of the system. For 
the particular case under consideration, the 
relative initial velocity of the system is virtually 
the same as the initial velocity of the high 
energy beam particle, the reduced mass of the 
system is equal to half the mass of the beam or 
scattering particle and the relative angle of 
scattering is twice the angular deflection of the 
beam particle with respect to the beam axis. _ 
For deflections of 10 degrees or less, the 
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Fic. 1. Dependence of cross section on voltage as calcu- 
lated for various assumed potentials: A, 8.71 exp(— 4.637) ; 
B, 8.71 exp(—4.63r!); C, 11.3 exp(—4.63r3); D, 8.71 
Xexp(— 4.6373) (units, 107° erg). Circles represent ex- 
perimental points. 


relative angle of scattering is given by 


(1) 





*[V(ro) — V(r) ]rdr 
—=f (r?—r?)3 


where V(r) is the intermolecular potential 
expressed as a function of the distance r between 
centers of the beam particle and the scatterer; 
ro, the distance of closest approach of the center 
of the beam particle to the center of the scatterer ; 
u, the reduced mass; and 2, the relative initial 
velocity of the colliding system. The total 
collision cross section, for a given average 
angular aperture 9, is given by 


ar f bedbe 
r) 


where G(@) is the scattering function and do, 
the impact parameter defined by 


V(ro 
(:- r Nr = be’. (3) 


1 
2 Mv* 


= the’, (2) 











So= 2 f GO sin 6d0}| = 
) 





Relations (1), (2) and (3) are sufficient to 
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predict total collision cross sections, consistent 
with an assumed intermolecular potential, for 
comparison with experimental values. The actual 
method of calculation is as follows: 

(A) A potential function V(r) is assumed. 
This permits the calculation from Eq. (1) of 
points for a curve of 3yuv? versus ro for a constant 
value of 6 (in the present case equal to 0.0785 
radian, the relative average angular aperture 
of the detecting system). From this curve, it is 
possible to find the 7» values which correspond 
to the experimental atom gun voltages as 
expressed in terms of 3,2". 

(B) Using the assumed potential form, the 
value of the impact parameter bo, corresponding 
to a given combination of ro and $v", is computed 
from Eq. (3). 

(C) From this value of bg the total collision 
cross section, for a given relative average angular 
aperture and beam particle energy, is directly 
computed from Eq. (2). 

It is therefore possible, in principle, to predict 
a curve of So versus voltage for any assumed 
intermolecular potential. Two general forms of 
repulsive potential were assumed, for both of 
which the integral in Eq. (1) may be evaluated 
without great difficulty. The first of these was the 
Lennard-Jones type which may be represented 
by V(r)=C/r* where s is a positive integer. In 
this case, it is possible, as shown by Kennard,’ 
to integrate Eq. (1) and obtain an expression for 
6 in closed form. The second general form of 
potential used was V(r)=A exp (—ar"). In this 
case, it can be shown that if 2/n is a positive 
integer, an asymptotic expansion of Eq. (1) 
exists from which 0 may be calculated with high 
accuracy. The general method is as follows: 

If V(r) =A exp (—ar"), Eq. (1) may be written 


” A exp (—aro")rdr 





e=— 
pv? J ro (r?—19")3 





(= A exp re ar Enid 


= 792)! 


Integrating the first term directly, the second 
term by parts and simplifying, gives 








e= 


2roAan f exp (—ar")r"—dr 
r0 (r?—19°)} 
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If y is a new variable such that y=a(r"—79"), 
then 


exp (—y) dy 


2/n 7 
(eat) 1] 
ary” 


When 2/n is a positive integer [(y/aro") +1 
may be expanded into a polynomial, after which 
the denominator of the integrand may be 
expanded according to the binomial theorem. 
The final result is an infinite series whose terms 
are I functions multiplied by known numerical 
coefficients. The binomial expansion is valid only 
for (y/aro")<1, so that the series does not 
converge. However, by taking only the first 
four or five terms, 8 may be computed to one 
part in several thousand since the first part of 
the series is an asymptotic expansion for 0. 
This expansion is very satisfactory for values of 
ary" from 3 to 6, the range encountered in the 
present calculations. 
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Fic. 2. Comparison of repulsive potential deduced from 
present experiments with Buckingham modification of 
Slater-Kirkwood potential. 


The repulsive potential which was first tried 
was the Lennard-Jones type suggested by de 


Boer and Michels,’ V(r) =4.45/r2X10- erg 
(rin A units). The resulting theoretical curve of 
Se versus voltage for this assumed potential 
was unsatisfactory both with respect to average 
magnitude and slope. The predicted cross 
sections decreased slowly from 6.6A? at 300 volts 
to 5.4A? at 1000 volts. Since satisfactory agree- 


‘de Boer and Michels, Physica 6, 409 (1939). 
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ment could not be obtained with this potential 
form by any change of constants, it was decided 
to use “‘softer’’ potentials of the form V(r) 
=A exp (—ar"). The modified Slater-Kirkwood 
repulsive potential proposed by Buckingham,‘ 
V(r) =8.71 exp (—4.63r) X10" erg gave results 
which, although better than those obtained by 
using the potential of de Boer and Michels, 
were still unsatisfactory. Changing the value 
of in the generalized potential expression from 
1 to 3 and 3 while keeping A and a constant at 
8.71 and 4.63, respectively, did not materially 
improve the agreement, but did serve to point 
out that if the So curve resulting from the 
assumption, V(r) = 8.71 exp (— 4.637?) K10-" erg, 
were laterally displaced in the direction of 
increasing voltage, the new curve would fit the 
experimental points quite well. This lateral 
displacement merely involves changing the 
assumed potential by a multiplicative constant 
which eventually multiplies the predicted voltage 
consistent with a given Sp. Applying this 
procedure in the present case, it is found that 
the repulsive potential which predicts collision 
cross sections in agreement with experiment may 
be represented by V(r) =11.3 exp (—4.637') 
X10-'° erg. Figure 1 illustrates the agreement 
between the experimental points and the pre- 
dictions resulting from assumption of the various 
exponential type repulsive potentials. 

The present results indicate that, for separa- 
tion distances of 0.55 to 1.05A, the repulsive 
potential is considerably softer than the Slater- 
Kirkwood or slightly modified Slater-Kirkwood 
potentials which are usually used for larger 
distances. Figure 2 compares the potential 
deduced from the present experiments with 
Buckingham’s modification of the Slater-Kirk- 
wood potential, at the distances involved in the 
scattering of the 300- to 1000-volt helium atoms. 

Margenau® states that the original calculation 
of Slater and Kirkwood’ neglected second-order 
exchange forces which have the effect of softening 
the repulsion. Near the minimum (2.9A) the 
magnitude of the repulsive energy is decreased 
about 9 percent by taking into account these 


5 Buckingham, Proc. Roy. Soc. 168, 264 (1938). 
6 Margenau, Phys. Rev. 56, 1000 (1939). 
7 Slater and Kirkwood, Phys. Rev. 37, 682 (1931). 
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second-order effects. Although Margenau’s 
calculations do not apply to smaller distances, 
it appears that the correction would be very 
much larger at such distances. 

The results of this paper show that the 
measurement of total collision cross sections of 
high velocity particles can be used as a sensitive 
method for determining intermolecular repulsive 
potentials if the experimental results reveal a 
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large dependence of cross section on velocity. 
In such cases it is found that the greater the 
slope of the cross section versus velocity curve, 
the softer will be the potential required to 
reproduce the experimental points. An extreme 
illustration of this point is found in the collision 
cross section for rigid spheres, which has no 
velocity dependence and is the result of an 
infinitely steep repulsive potential. 
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A study has been made of the effect of an inert solvent 
on the efficiency of the primary process in the photolysis 
of oxalyl bromide and nitrosyl chloride. The oxalyl bro- 
mide was studied at the wave-lengths 2650, 3130, 3650, 
4047, and 4358A, the nitrosyl chloride at 3650, 4047, 4358, 
5461, and 5790A. In order to eliminate the back reaction 
between nitric oxide and chlorine, oxygen was introduced 
which converted the nitric oxide into nitrogen dioxide. 
The solvent used was carbon tetrachloride. Quantum 
vields for nitrosyl chloride showed only a slight trend 
ranging from 0.7 at 3650A to 0.5 at 5790A as compared to 
approximately 2 at all wave-lengths in the gas phase. 
Oxalyl bromide gave approximately the same quantum 


HE effect of a solvent on the efficiency of a 

primary photochemical process has been 
discussed by Franck and Rabinowitch.' Accord- 
ing to their view the quantum yield of the initial 
step should be lower in a solvent than in the gas 
phase because of (1) deactivation by solvent 
molecules, or (2) a higher efficiency of recom- 
bination due to the fragments of the original 
molecule being held together by a ‘“‘cage”’ of 
solvent molecules. The chief difficulty in applying 
these views to most systems lies in the fact that 
secondary processes are involved and it is im- 
possible to determine how much of an observed 
solvent effect is caused by a modification of the 


secondary reactions. In order to avoid such com- — 


plications we have studied two systems in which 
the fragments produced by photodissociation of 
the absorbing molecule do not react with the 


wa and Rabinowitch, Trans. Faraday Soc. 30, 120 
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yield in solution as in the gas phase at 2650A, i.e., ~ 0.9, 
and the values in solution decreased to ~0.3 at 4358A 
with most of the change occurring between 3130 and 3650A. 
The results are discussed in terms of the Franck-Rabino- 
witch hypothesis and it is pointed out that they are most 
readily accounted for by assuming a deactivating action 
of the solvent rather than a high efficiency of recom- 
bination. A new quartz capillary arc of high intensity is 
described. The equilibrium N2O,=2NOz has been studied 
in carbon tetrachloride solution and the equilibrium 
constant at 25°C found to be 6.4X10-5, concentrations 
expressed in mole fractions. 









solvent molecules and any secondary reactions 
are unaffected by the solvent or do not affect 
the yield of the primary process. The reactions 
are (1) the photolysis of oxalyl bromide and (2) 
the photo-oxidation of nitrosyl chloride. Both 
were studied in carbon tetrachloride solution and 
some experiments were performed with nitrosyl! 
chloride in hexane solution in which case action 
with the solvent had to be considered. 


EXPERIMENTAL METHODS AND APPARATUS 


The light source (Fig. 1) used in most of the 
experiments was a thick-walled capillary arc. 
Quartz capillary of about 1.5 mm wall thickness 
and internal diameter between 2.7 and 3.3 mm 
has been used on various arcs; a bore of around 
3 mm and a length of about 4.5 cm for the 
working part seemed to give the best results. 
The bottom tungsten lead is sealed into the large 
quartz tube with De Khotinsky cement, and the 
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arc filled with mercury to the upper bulb. The 
tube is pumped on to remove most of the gas 
trapped on the walls, though there is no need for 
great care in removing air, since gas bubbles do 
not tend to collect in the capillary. After the 
rest of the tube is filled with mercury and the 
screw clamp tightened on the rubber tubing, 
mercury may move slowly through it, but sudden 
oscillations are prevented. The upper part of the 
holder into which the arc is screwed may be 
filled with water, so that the rubber and upper 
part are kept cold. The arc was operated on 220 
d.c. (the top is connected to the negative side 
of the line) with a resistance in series. The re- 
sistance is adjusted to give a current of around 
10 amperes, and the arc struck by breaking the 
mercury column near the bottom bulb with a 
small flame. A brisk stream of cooling water 
should be turned on immediately and the jacket 
(a quartz tube for filter work—a metal tube 
with a quartz window for the monochromator) 
screwed up so that the arc is entirely surrounded 
by flowing water. During operation, the re- 
sistance in series is varied between 10 and 15 
ohms, depending on the arc’s characteristics and 
the intensity desired. Around a kilowatt was the 
average power used, though arcs have been 
operated on as much as 1500 watts for several 
hours. They will not operate well on 110 d.c. 
with little or no resistance in series. After forty 
or fifty hours of operation it is advisable to clean 
the capillary with hydrofluoric acid and repair 
any small cracks which may develop. None of 
these arcs have ever broken down violently 
except when the water supply failed. One was 
operated (with perhaps a dozen or more clean- 
ings) almost 1000 hours on about 900 watts. A 
tungsten lead which has been used for a long 
time in the bottom of an arc seems to help 
prevent blackening of the capillary. 

Zeiss monochromat filters A, B, and C were 
used for the 5790, 5461, and 4358A mercury 
lines, respectively, Corning 586 filter (8 mm) for 
3650A, and Corning 587 (3.9 mm) with Corning 
330 (noviol 0, 2 mm) for the 4047 A line. A large 
fused quartz prism monochromator was used to 
separate the 2650 and 3130A radiations. 

Nitrogen dioxide was prepared by heating C. P. 
lead nitrate. The gas was passed through phos- 
phorus pentoxide and freed from gaseous im- 
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purities by repeated freezings and pumpings. The 
solid nitrogen tetroxide was colorless. Apiezon 
grease “‘L’’ was not attacked noticeably by 
short exposures to nitrogen dioxide. A weighed 
ground joint with a small bulb attached was 
pumped out in the vacuum line. A known volume 
and pressure of nitrogen dioxide was condensed 
into the bulb which was then sealed off and 
weighed, along with the rest of the ground joint. 
Weights of nitrogen dioxide thus obtained 
checked to within less than 1 percent of the 
weights calculated from the volume-pressure 
measurements, using the equilibrium constants 
for nitrogen dioxide and nitrogen tetroxide of 
Verhoek and Daniels.* 

Chlorine was prepared from tank chlorine, 
which was bubbled into water at 0°C until a 
considerable amount of chlorine hydrate had 
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formed. On warming to 10°C the gas was evolved 
and passed through calcium chloride and phos- 
phorus pentoxide. The chlorine was distilled in 
a vacuum line, the first and last fractions being 
discarded. 

Nitric oxide was prepared by a modification 
of the method of Johnston and Giauque.* 


( — and Daniels, J. Am. Chem. Soc. 53, 1250 
1931). 

§ Johnston and Giauque, J. Am. Chem. Soc. 51, 3194 
(1929). 
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Nitrosyl chloride was prepared by a method 
similar to that of Taylor and Denslow.* 

Nitric oxide (in excess) and chlorine were 
allowed to react in a three-liter flask. When the 
chlorine was practically all used up the gases 
were repeatedly condensed in liquid air and 
evaporated until pressure measurements showed 
complete reaction. Then the gases were con- 
densed in liquid oxygen and non-condensible 
nitric oxide pumped off, the operation being 
repeated until no gas was uncondensed by liquid 
oxygen. Finally the nitrosyl chloride was dis- 
tilled, with the first and last fractions discarded. 
Some preparations were made in a system using 
greaseless glass valves; others were made in 
a system containing stopcocks greased with 
Apiezon grease ‘‘L.”” The grease, which was 
protected from light when in contact with 
nitrosyl chloride, was not affected visibly. No 
effect of stopcock grease on the gas photolysis 
was found. Weights of nitrosyl chloride were 
checked against pressure-volume measurements 
by the method used with nitrogen dioxide. These 
different values checked to within about } 
percent. 

Oxalyl bromide, prepared from oxalyl chloride 
and hydrogen bromide in a system using grease- 
less glass valves, was kindly supplied by Dr. 
Tuttle. Titration of the acid formed on hydrolysis 
of oxalyl bromide checked within about 3 percent 
of the pressure-volume measurements. Free 
halogen in the material amounted to less than 
0.1 percent. 

Sulfur-free carbon tetrachloride, saturated 
with chlorine, was exposed to sunlight for several 


4 Taylor and Denslow, J. Phys. Chem. 31, 374 (1927). 


ATWOOD AND G. K. 








ROLLEFSON 








days, but was found to produce acid on illumi- 
nation with chlorine even after such treatment. 
Finally chlorine dioxide was bubbled through and 
the solution illuminated with a 150-watt bulb 
for five days, chlorine dioxide being renewed as 
it was used up. Chlorine dioxide was prepared 
by the method of Bray.® Finally the carbon 
tetrachloride was washed with alkali and distilled 
from phosphorus pentoxide, the end fractions 
being discarded. 

Crude hexane was shaken with small amounts 
of concentrated sulfuric acid until the acid took 
up no color. Then the hexane was shaken several 
times with 0.1 N potassium permanganate in 10 
percent sulfuric acid, and a few times with 0.1 V 
potassium permanganate in 10 percent sodium 
hydroxide solution. The fraction distilling be- 
tween 65-70°C was used. 

Pressures of the corrosive gases were measured 
with a click gauge connected to a sulfuric acid 
manometer. Samples of the various gases were 
condensed and sealed off in small bulbs which 
could be broken by shaking in the reaction cells. 
Two cells were used—one cylindrical Pyrex cell 
40 mm in diameter and 30.4 mm long, the other 
a rectangular quartz cell 30 by 40 mm and 24.2 
mm long, used in the work below 3650A. 

The experimental set-up is shown in Fig. 2. 
Absorption of substances in the cell could be 
measured in two ways. Direct reading of currents 
from photo-cell P; with the cell in and out of the 
light path may be made, or both photo-cells may 
be used in a potentiometric arrangement. The 
beam reflected from the quartz plate to photo- 
cell P2 measures the incident light intensity ; the 
voltage of P2 is applied across the full resistance 
of the potentiometer (resistance boxes 2 and 3). 
The two methods gave identical results, but it 
was more convenient to use the first method 


TABLE I. Calibration of the sensitivity of the photronic cell. 














(MOLES/MIN. XMM) X106 
WaAVE-LENGTH, NITROSYL URANYL 
A CHLORIDE OXALATE THERMOPILE 
3650 0.57 0.49 0.59 
4358 1.22 1.44 1.08 
5461 5.0 — 6.8 
5790 7.8 — 7.3 





















5 Bray, Zeits. f. physik. Chemie 54, 463 (1906). 
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when the cell could be moved rapidly back and 
forth from the light path. When the cell could 
not be moved, it was necessary to use the latter 
method, since it corrects for fluctuations in the 
intensity of the arc. When absorption measure- 
ments were used for analytical purposes dense 
screens, cutting out about 97 percent of the light, 
were used, so that photochemical reaction during 
the measurement would be negligible. Analytical 
measurements were often made at one or two 
wave-lengths on a reaction initiated by light of 
another wave-length. 

The absolute number of quanta entering the 
reaction cell was determined from current outputs 
of photo-cell P;, which was calibrated against 
actinometers and a thermopile. Uranyl oxalate® 
was used as an actinometer between 2537 and 
4358A. This actinometer alone was used at 2650 
and 3131A. Nitrosyl chloride gas, assuming a 
quantum yield of 2.0, was used between 3650 
and 5790A, and those calibrations checked 
against a Moll large surface thermopile. Photo- 
cell sensitivities from these actinometers and the 
large surface thermopile calibrated against a 
standard lamp are compared in Table I. The 
nitrosyl chloride gas values were used in calcu- 
lating the nitrosyl chloride solution yields, while 


TaBLE II. Oxalyl bromide absorption coefficients in carbon 
tetrachloride solution. 





Wave-Lenctn, A 


2650 
3130 
3650 
4047 
4358 





~120 
30.9+1.7 
5.9+0.2 
1.78 (2.43-1.58) 








M =concentration of oxalyl bromide in moles/liter. 
/=length of cell. ° 
K =1/1M Xlogio Io/I. 








TaBLe III. Quantum yields for the photolysis of oxalyl 
bromide in carbon tetrachloride solution. 








AVERAGE 


0.89 
0.79 
0.37 
0.26 
0.29 


OBSERVED VALUES 


, 0.82, 0.91, 0.87 
, 0.81, 0.74 
, 0.34, 0.38, 0.38, 0.35 


Wave-LENGTH 


2650 
3130 
3650 
4047 
4358 











52, 3139 


a and Forbes, J. Am. Chem. Soc. 
930), 


PROCESS 


TABLE IV. 








X105 T =25°C 
Nnw, 


STOICHIOMETRIC MOLE FRACTION OF K= 
NOz X108 





5.65 
5.69 
5.55 
6.18 
6.15 
6:34 
6.42 
6.41 
6.40 
6.49 
5.68 
5.73 
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145.0 
173.9 





4250 1.4* 





- trations. 


NO. =74.6 
NO, = 2.5 


Henry’s law constants 
_ fugacity (atmos. ) T=25°C 


mole fraction 











* Cundall, J. Chem. Soc. 51, 1076 (1891); 67, 794 (1895). 


uranyl oxalate values were used for the oxalyl 
bromide quantum yields. 


PHOTOCHEMICAL DECOMPOSITION OF OXALYL 
BROMIDE IN CARBON TETRACHLORIDE 
SOLUTION 


The decomposition of oxalyl bromide was fol- 
lowed by bromine absorption measurements, and 
iodometric titration of the bromine liberated. 
Since oxalyl bromide has a much smaller absorp- 
tion coefficient than bromine at 4358A, light 
absorption at that wave-length was used for 
analysis in the work at all wave-lengths. Titration 
and absorption measurements gave checks within 
5 percent, except in runs at 2537 and 2650A. At 
these wave-lengths the carbon tetrachloride 
itself absorbs light and liberates oxidizing ma- 
terial. As the concentration of oxalyl bromide is 
increased, the discrepancy decreases, since the 
solvent then absorbs a smaller fraction of the 
light. If a correction is made for oxidizing ma- 
terial produced from the solvent the two 
methods of analysis are brought into good 
agreement, even at low oxalyl bromide concen- 
The photolysis of oxalyl bromide 
appears to be unaffected by the accompanying 
carbon tetrachloride decomposition, except inso- 
far as the solvent acts as an inner filter. The inner 
filter action of the bromine at the longer wave- 
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(Moles /Liter) x10? 
Fic. 3. 3650A. --—-—-— (A log Io/I); vs. stoichiometric 


NO2z conc.; ———— (A log Jo/Z)243 vs. NOCI conc.; 
------ A log Io/I vs. NOCI conc. decrease. 


lengths was eliminated by using initial rates for 
the calculation of quantum yields. 

Two runs were made at 2537A, using a low 
pressure, gas-filled mercury arc which furnishes 
over 90 percent of its radiation at 2537A. 
The set-up in these two experiments was different 
from that in the others. A lower limit for the 
quantum yield was about 0.6; the geometry of 
the system indicated an upper limit of 1.2. 

Errors in the 4358A measurements are larger 
than those at the other wave-lengths, since the 
absorptions are small and the quantum yield is 
low. There appears to be a slight change in 
absorption coefficient with concentration at 
4358A, though the errors in measurement at the 
low concentrations are large. Table II shows the 
absorption coefficients of oxalyl bromide in 
solution, Table III summarizes the quantum 
yield data. 

Illumination of a solution of nitrosyl chloride 
in carbon tetrachloride caused no measurable 
amount of decomposition. Experiments with 
mixtures of nitric oxide and chlorine in solution 
in the same solvent showed that the back reac- 
tion was so fast (at least 100 times as fast as in 
the gas phase) that the photostationary state 
which could be set up by our light source was 
not measurable. In order to eliminate this back 
reaction oxygen was introduced into the system. 
Control experiments showed that oxygen reacts 
very much faster than chlorine with nitric oxide 





G. &. 











ROLLEFSON 


in carbon tetrachloride. With oxygen present the 
net reaction on illumination is 2NOCI+0O, 
=2NO.+Cl.. This reaction can be followed 
readily by measurement of the increase in the 
absorption of light by the solution. The presence 
of any water in the carbon tetrachloride inter- 
fered with this measurement by reacting with 
the nitrogen dioxide. In order to eliminate this 
difficulty it was necessary to distill the carbon 
tetrachloride through fluffy phosphorus pen- 
toxide and condense it directly in the reaction 
vessel which was then sealed off. 

The nitrogen dioxide exists in the solutions in 
equilibrium with the tetroxide. From absorption 
measurements in very dilute solutions, in which 
the dissociation of the tetroxide is essentially 
complete, it was possible to determine the absorp- 
tion coefficients for the dioxide at the wave- 
lengths 4358 and 3650A. Since only the dioxide 
absorbs at 4358A, it is possible to obtain the 
data necessary for the calculation of equilibrium 
constants by measuring the absorption and the 
total amount of dissolved oxides at higher con- 
centrations. The constants obtained in this way 
are tabulated in Table IV. 

In order to follow the course of the oxidation 
of nitrosyl chloride by light absorption it is 
necessary to set up a relationship between the 
change in the light absorption and the amount 
of the nitrosyl chloride which has reacted. The 
change in log J)/J may be written 


A log Io/I=ki ACyo, +ReACci, +RsACyoci- 


From the stoichiometry of the net reaction 
ACci,= —F4Cyor:. The correlation between 
ACyo, and ACyoc is not so simple because of the 
equilibrium between the dioxide and tetroxide. 
However, if we speak of the total tetravalent 
nitrogen as a stoichiometric Cyo, then the ACyo, 
corresponding to this is numerically equal to 
ACxoc). Thus the above equation can be reduced 
to a function of ACyoc). The simplest way to use 
this expression for analytical purposes is to 
prepare a graph of A log Io/I against ACyou. 
Each term of the equation is plotted; the las. 


‘ two give straight lines but the first gives a curve 


because of the equilibrium between the dioxide 
and tetroxide. The contribution of each term for 
a given ACyoc; is read from the graph and the 
three are combined to give a point for the 
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desired curve. The curve obtained with \4358A 
could be used only for small values of ACyoc: 
because A log Io/J changed very little at the 
higher values. The curve obtained for \3650A 
was quite satisfactory and was used in most of 
the analyses with checks against the one for 
\4358A. The agreement between the two was 
found to be satisfactory. Figure 3 shows the 
curves for \3650A. In this figure the contribu- 
tions of the second and third terms have been 
combined into the single straight line. 

The quantum yields for the oxidation reaction 
are listed in Table V. It was found that unless 
the oxygen pressure over the solution was con- 
trolled considerable variation was noted in the 
results. The effect of varying the oxygen pressure 
is shown in the values for 3650A. High oxygen 
pressures apparently have some inhibiting action 
although oxygen must be present for the reaction 
to go. At the other wave-lengths the oxygen 
pressure was fixed at one atmosphere. 

The recombination of the nitric oxide and 
chlorine atoms formed in the photolysis of 
nitrosyl chloride may be prevented by carrying 
out the reaction in a solvent which reacts with 
chlorine atoms. Hence some experiments were 
performed with hexane as a solvent. No nitrogen 
dioxide was formed in this system under any 
conditions and the absorption decreased as the 
reaction proceeded. The reaction could not be 
run to completion since the solution became 
cloudy due to the separation of some hexane- 
insoluble product. Some experiments were tried 
with a solution of hexane in carbon tetrachloride 
and in these the reaction was carried to com- 
pletion; the absorption dropped to zero. The 
quantum yields for these experiments are con- 


tained in Table VI. 


TABLE V. Quantum yields for the oxidation of nitrosyl 
chloride in carbon tetrachloride solution. 





WaAvE- 
LENGTH 





OXYGEN 


PRESSURE OBSERVED VALUES AVERAGE 


3650 





0.70 
0.66 
1.5 ‘ 0.56 
2 F r 0.47 
1 y 76 0.72 
1 
1 


0.2 atmos. 


0.47 
0.50 
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TABLE VI. Quantum yields for the photolysis of nitrosyl 
chloride in hexane solution. 





WAVE-LENGTH OBSERVED VALUES 


3650 0.91, 0.87, 0.96, 1.19, 1.10 1.0 
4358 0.69 0.69 
0.2 
0.5 


AVERAGE 





5790 7, 0.30, 0.33, 0.32, 0.37 0.32 
3650 0, 0.60 0.55 
Solvent 0.05 mole 
fraction of hexane 
in CCl, 





DIsCUSSION 

According to the hypothesis of Franck and 
Rabinowitch! the quantum yield of simple reac- 
tions such as are considered here should be less in 
liquids than in the gaseous state. Our results 
show that this is true at all wave-lengths in the 
case of nitrosyl chloride and at long wave- 
lengths with oxalyl bromide but not at short 
wave-lengths where our results are essentially 
the same as found by Tuttle in the gas phase.’ 
If the decrease in yield of the primary process is 
caused by a high efficiency of recombination of 
the molecular fragments brought about by a 
“cage” action of the solvent molecules it might 
be expected that the yield would vary with 
wave-length since at high frequencies the large 
kinetic energies possessed by the separating 
fragments would favor escape from the cage. The 
excess energy in the case of nitrosyl chloride 
amounts to 41 kcal. at 3650A and 11 at 5790A.° 
The products are the same at all wave-lengths 
except possibly that the chlorine atom may or 
may not be raised to the P, state, a difference of 
2.5 kcal. The estimate of the excess energy for 
the oxalyl bromide depends upon what values 
are chosen for the carbon-carbon and carbon- 
bromine bonds. If we assume that the primary 
process is the same between 3650A and 4358A, 
the range in which an effect was noted, the vari- 
ation in energy amounts to 13 kcal. Inspection 
of our results shows that in carbon tetrachloride 
solution there is only a slight trend in the yield 
with nitrosyl chloride and with oxalyl bromide 
the only change of importance occurs in a rather 
narrow wave-length range. Nitrosyl chloride in 
the presence of hexane shows a much greater 
variation. In this case a reaction with the solvent 


7J. E. Tuttle, Ph.D. Thesis, University of California, 
40. 
8 Dixon, Zeits. f. physik. Chemie B24, 328 (1934). 
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occurs and it is possible that the excess energy 
supplied at the shorter wave-length is required 
to supply an activation energy for such a reac- 
tion. A point in favor of this view is supplied by 
the experiments in the dilute solution of hexane 
in carbon tetrachloride. Under these conditions 
some of the excess energy could be taken up by 
the carbon tetrachloride before a reaction with 
hexane could occur and actually we find a 
quantum yield only half as great as in pure 
hexane at the same wave-length. 

The results which have been presented are 
more readily accounted for by assuming a deac- 
tivating action of the solvent rather than a high 
efficiency of recombination. Such a deactivation 
need not involve a change of electronic state in 
the photoactivated molecule but only the 
removal of vibrational energy. Usually the ab- 
sorption of light produces a strongly vibrating 
molecule; if vibrational energy is removed from 
such a molecule until the amount left is less than 
that required for the dissociation of that par- 
ticular electronic state the molecule will remain 
intact. If the potential energy curve which could 
be drawn to represent the excited state has no 
minimum or only a very slight one it may not be 
possible to remove enough energy to prevent 
dissociation. Apparently such a situation is 
found with hydrogen iodide in hexane solution.® 
If we consider oxalyl bromide from this stand- 
point the difference in behavior at long and at 
short wave-lengths corresponds to excitation to 
two different electronic states. The one produced 
at short wave-lengths has no minimum in its 
potential energy function, hence the solvent 
cannot prevent the dissociation and the quantum 


® Warburg and Rump, Zeits. f. Physik 47, 305 (1928). 
They found the quantum yield the same in the gas phase 
and in solution at 2820A. At shorter wave-lengths their 
= were complicated by absorption of the solvent they 
used. 
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yield is the same in solution and in the gas phase. 
At long wave-lengths two-thirds of the molecules 
in solution lose so much vibrational energy that 
they cannot dissociate and the quantum yield 
drops accordingly. In the case of nitrosyl chloride 
between two-thirds and three-fourths of the 
molecules are deactivated. The slight trend in 
the results is such that it could be said that with 
the larger energies there is a little less chance of 
the energy being lowered below the dissociation 
limit. According to Goodeve and Katz!® the 
potential energy curves for the upper states of 
the nitrosyl chloride involved in the absorption 
process at all the wave-lengths we have used 
except 5790A are repulsion curves. However, 
these curves intersect one which does have a 
minimum and transfers of the type which are 
involved in predissociation occur between the 
states. If this set of energy levels is assumed to 
be correct then the deactivating action of the 
solvent would involve removal of vibrational 
energy plus this transfer to an excited state with 
a minimum in its potential energy function. 
Ultimately the molecule is deactivated further to 
the lowest state but the dissociation is prevented 
in the first action of the solvent. 

The arguments which have been presented do 
not exclude the possibility that a great many 
cases may exist in which a low quantum yield in 
solution is due to recombination. That effect has 
probably been overemphasized in the literature. 
A deactivation by the removal of vibrational 
energy is better suited to account for some of the 
specific effects observed and the small or absent 
trends with wave-length. It also offers a pos- 
sibility of predicting the effect of a solvent if any 
information is available concerning the potential 
energy function of the excited state. 


1 Goodeve and Katz, Proc. Roy. Soc. London A172, 432 
(1939). 
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The equations for sound waves passing through liquid mixtures are given. The amount of 
unmixing and the increase in sound absorption are calculated and found to be too small for 


observation. 





W. WILLARD! has found that the ab- 
e sorption of ultrasonic waves in liquid 
mixtures may have a maximum for a certain 
concentration. In a discussion of Willard’s paper 
at the May, 1940, meeting of the American 
Physical Society, the writer proposed an explana- 
tion, which Willard has kindly mentioned in his 
paper, namely that sound waves produce a 
partial unmixing.? 

Unfortunately, the detailed calculation shows 
that both the unmixing and the increase in 
absorption are much too small to be of interest, 
and therefore only a short report is given here. 

The liquid mixture is assumed to consist of 
two substances of molecular weight M, and M2, 
and concentrations (mole per cc) ¢; and Co. 
The particle velocity is called “4; and up». 

Let a plane sound wave, of frequency v and 
pressure amplitude P travel through the medium 
in the z direction. Let yu; and pe be the chemical 
potentials of the two components in the mixture. 
The equation of motion of the component one 
is then for one mole. 


~) 
0z T 


kr 


Ou, 


Ci +co.0T T; 





(1) 


Ci Oz (uy; —us) 


Here kr is the coefficient of the thermal 
diffusion, as defined by Chapman.’ The last term 
is the force of friction resisting the relative 
motion of the two components. 

Introducing the mole fraction x 


x=€;/(C1 +62), 
this can be written, together with the equation 


'G. W. Willard, Phys. Rev. 57, 1057 (1940); J. Acous. 
Soc. Am. 12, 438 (1941). 

? Y. Rocard, J. de phys. et rad. 1, 426 (1930). 

*S. Chapman, Proc. Roy. Soc. 177, 38 (1940). 


for the second component 
Ou, RT kr\ Op Om Ox 
M,—=- ( Vi-—-# ‘) —_—— — 
ot Cp x 702 Ox 02 


2RT 
+—(1—x)(wi—wu2), (2) 
D 


OU» RT kr \op 
M,—= - (145 -B )~+ 
ot Cc,» 1—x/Jdz 1—x Ox dz 


x Op, Ox 


2RT 
———(u,— Ue) x. 
D 


(2') 


Here, V;, V2 are the partial mole volumes in the 
mixture, 8 the expansion coefficient and c, the 
specific heat at constant pressure of 1 cc of the 
mixture, so that 


dT /dz=(TB/c,)dp/dz. (3) 


On the right side of (2) the first term is the 
driving force of the pressure gradient, the second 
that of the thermal diffusion, the third that of 
the concentration gradient and the last that of 
friction. The normal effect of viscosity has been 
omitted, as at first approximation the effects are 
additive. D is defined so that the flow due to a 
concentration gradient is 


D Oui 
RT(i—x) dz 
Furthermore, the energy equation is written 


oT dp aT dHax 
¢y-—=TB—+K——— —. (4) 
at dt ds ax at 


The first term on the right is the useful work 
done, the second is the heat evolved in connection 


with thermal diffusion, the last one is due to a 
possible deviation from the law of ideal solutions, 
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representing a (negative) heat of dilution. Heat 
conductivity and slow energy exchange terms 
have been omitted. 

However, in the integration we neglect, in the 
energy equation, the last term, and in the equa- 
tion of motion the effect of the concentration 
gradient. 

We call the sound velocity W, the value 
without unmixing Wo, 2a’ the increase of the 
absorption coefficient 2a, defined so that the 
intensity diminishes proportional to exp (—2az), 
Ax the change of x at its maximum. 

We also introduce the notation 


pa=M, Vi, pe= M;2/ V2. 
Then 


9 


K 1 | a 
-— pot (14 ) (vx=0), (5) 
Wo c Woe p 


2 
Pp 
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2a’ 2r? M?M27c1c2 D p2— pi 





2 WaeMit+A—x)M2 RT pips 





po—pi RTB_ kr 1—-x <x 
Pipe cp x(li-x)\ M, Mz, 


P wDysi-x x\! 
Ax =x(1-—x)— a ms +) 
Ww? RT M, Ms 


po—pi RTB kr li—zx «£ J 
x|=—" -—— ran -(—+— . () 
Pip2 Cp x(1 —x) M, M, 
Equations (6) and (7) are given in first approxi- 
mation. The phase of Ax is 90° ahead of the 
pressure. 
Numerically, 2a/v* is, for a reasonable case, 


about 10-'8, Ax/x, for P one atmosphere pres- 
sure, 1-107!8p. 
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A statistical mechanical theory of fusion based upon the use of local free energies is pre- 
sented. An integral equation is formulated for the distribution function of average density in a 
region occupied by a system of molecules. Periodic solutions characteristic of a crystalline phase 
are found for certain ranges of values of a set of parameters depending upon temperature and 
volume. When the parameters decrease below certain critical values, all terms of the Fourier 
series representing the distribution function vanish with the exception of the constant term. 
A uniform density distribution characteristic of a fluid phase is then obtained. The melting 
parameters of argon at several pressures are calculated with the aid of the theory and compared 


with experiment. 


INTRODUCTION 


N recent investigatious of the fusion process 
by the methods of statistical mechanics, two 
different points of view have been advanced. 
Lennard-Jones and Devonshire have proposed a 
very interesting theory based upon the idea of 
an order-disorder transformation in the distribu- 
tion of molecules between normal and interstitial 
sites of a basic lattice.“* Born* has proposed a 
1 J. E. Lennard-Jones and A. F. Devonshire, Proc. Roy. 
Soc. A170, 464 (1939). 


2G. H. Wannier, J. Chem. Phys. 7, 810 (1939). 
3M. Born, J. Chem. Phys. 7, 591 (1939). 


theory of a different type in which the liquid- 
solid transition is attributed to the vanishing of 
the rigidity modulus of the solid and a resulting 
instability of the lattice under shear. Although 
the two points of view may not be mutually 
exclusive, that of Lennard-Jones and Devonshire 
is intuitively the more attractive. From the 
diffraction of x-rays, it is known that a crystalline 
solid is characterized by a periodic density 
distribution extending over domains of macro- 
scopic dimensions. In the scattering of x-rays by 
liquids, there is no evidence of long range 
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crystalline order, but only of a rapidly damped 
local order in the density distribution around 
each molecule. However, the Lennard-Jones- 
Devonshire theory leaves something to be 
desired in the details of its formulation, since it 
allows the disordered phase a residual long range 
order rather than requiring a complete destruc- 
tion of the lattice in the fusion process. This is 
due to the limitations of the formalism of the 
order-disorder theory, which is designed for the 
solution of problems involving a finite set of 
configurations such as superlattice formation 
in alloys. 

We shall present a theory of fusion similar in 
its physical aspects to that of Lennard-Jones 
and Devonshire. It is based upon the use of 
local free energies in a classical canonical 
ensemble. With aid of a device based upon the 
introduction of coupling parameters into the 
potential of intermolecular force, an integral 
equation is formulated for a distribution function 
oi(R) specifying the average density at a point 
R in a region occupied by a system of molecules. 
The kernel of the integral equation contains, in 
addition to the potential energy of a molecular 
pair, a correlation function g(R,R’) which 
determines the local order in the density at a 
point R’ established by a molecule situated at 
the point R. The correlation function is approxi- 
mated by the radial distribution function, to 
which it reduces in a liquid phase. For density 


distributions with the period of a specified space’ 


lattice, the integral equation is shown to be 
equivalent to a set of transcendental equations 
for the Fourier coefficients of the distribution 
function p,(R). As parameters in the set of 
transcendental equations there enter the Fourier 
transforms of the potential energy of a molecular 
pair and the correlation function. For certain 
ranges of values of the parameters, it is found 
that all terms of the Fourier series representing 
the distribution function p;(R) vanish except 
the constant term. When p;(R) reduces to a 
constant, the average density becomes uniform 
with no residue of crystalline long range order. 
A disordered density distribution of this type is 
characteristic of a liquid phase. For values of 
the parameters exceeding certain critical values, 
the transcendental equations determining the 
Fourier coefficints of p;(R) have non-vanishing 
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solutions, and a non-trivial periodic density 
distribution characteristic of a crystalline phase 
exists. The stable periodic structure is selected 
as the one of minimum work content. Since the 
parameters determining the distribution function 
pi(R) depend upon temperature and volume, 
areas in the temperature-volume plane corre- 
sponding to the liquid and solid phases are 
determined. Since it is convenient to calculate 
the distribution functions under the restraint of 
macroscopic homogeneity, areas corresponding 
to coexistence of the liquid and solid phases are 
determined with the aid of the thermodynamic 
criteria of heterogeneous equilibrium. 

The melting parameters of argon are calcu- 
lated with the aid of the theory and compared 
with experiment. An interesting inference from 
the theory is the absence of a critical point for 
fusion for certain crystal structures, among them 
the face-centered cubic structure characteristic 
of argon. 


THE DISTRIBUTION FUNCTIONS 


In a system of N molecules confined within a 
volume v, the average molecular density at a 
point R in this region is equal to (N/v)p,(R), 
where pi(R;)dv/v is the probability that any 
specified molecule 7 occupies dv at the point R,. 
According to classical statistical mechanics, the 
distribution function is determined by the local 
free energy W,(R) of the system when a molecule 
is fixed at the point R. 


pi(R) =exp | —(1/kT)[Wi(R) — Wo J}, 
exp [— W,(Rz: . -R))/kT ] 


=o f o- f exp (- Vy RT )dv,: 2 -dvn—n; (1) 


where kT is the product of Boltzmann's constant 
and the temperature, Vy is the potential of 
intermolecular forces, and W,(R,, ---R,) is the 
local free energy of the system when a set of n 
specified molecules k---/ are fixed at the points 
R,, -:-Ri. 

A crystalline phase is characterized by a 
periodic distribution function p,(R) with the 
period of a certain space lattice possessing 
maxima at specified points within the elementary 
cell A of the lattice. A liquid phase, on the 
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other hand, is characterized by a _ constant 
distribution function equal to unity everywhere 
in the region v. Coexistent phases would corre- 
spond to distribution functions p:(R) with 
physical discontinuities at the phase boundary. 
A theoretical investigation of the properties 
of the distribution function p:(R) might be based 
upon a direct evaluation of the integrals, Eq. (1) 
defining the local free energies Wi(R) and W,. 
Since this is an almost prohibitively difficult 
procedure in condensed phases, we shall employ 
a device, previously proposed by one of us,‘ 
which involves the introduction of coupling 
parameters into the potential of intermolecular 
force. We assume that Vy can be expressed as 
a sum of the potential energies V(|R,—Rz|) of 
molecular pairs, each term of which depends 
only on the distance |R,—R,| between the 
members of the pair. 
Vy= V(|Ri—R.)). 


k<l 


(2) 
We may define a potential Vy(é), 

Vy (&) = Vn rtéiVi, 
DY V(|R.—R.|), 


k<l$i 


Vi=X V(|Re—Ri|), 
k 


(3) 


Vy-i1= 


corresponding to a hypothetical system in which 
all terms in the potential energy associated with 
a specified molecule 7 are multiplied by a 
parameter £;. By allowing &; to assume values 
between zero and unity, the coupling between 
molecule 7 and the other molecules may be 
varied from zero to its value in the actual 
system. Corresponding to the potential Vy(é), 
we define the local free energies 


’ —W, (Rx, -+-Ri, &)/kT 
exp [ ( ty &:)/kT ] (4) 


=f or f exp [- Vy(&:)/RT Jdv1- ° *dvy—n, 


where the integral extends over the coordinates 
of all molecules except those of the set of n 
fixed molecules k---/. According to Eq. (1), 
the distribution function p;(R, £) in the partially 
coupled system becomes 


pi(R, £) =exp { — (1/kT) [WA(R, &) — Wo(é) J}, (5) 
‘J. G. Kirkwood, J. Chem. Phys. 3, 300 (1935). 
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where the parameter & is understood to refer to 
a specified molecule fixed at the point R. For 
the determination of the distribution function 
pi(R, £), and for the calculation of the thermo- 
dynamic functions of the system, we shall need 
to consider the following distribution functions 
for molecular pairs,® 


p2(R, R’, £) 

=exp { —(1/kT)[W2(R, R’, £) — Wo(é) J}, 

 p2(R, R’, £) 

=exp {| —(1/kT)LW2(R, R’, t) — Wi(R, £) J}, 

g(R, R’, )=exp {—(1/kT)LW2(R, R’, é) 
+Wo(é)— Wi(R, §)— Wi (R’, £) J}, 


where é is understood to refer to a molecule at 
the unprimed position R. The probability that 
the partially coupled molecule is situated in dv 
at R and another specified molecule in dv’ at R’ 
is (1/v*)po(R, R’, )dvdv’. The probability that a 
specified molecule is situated in dv’ at R’ when 
the partially coupled molecule is fixed at R is 
(1/v)pe(R, R’, —)dv’. From Eqs. (5) and (6), 
the following relations between the distribution 
functions are evident, 


p(R, R’, £) = pi(R, £)pi(R’)g(R, R’, E ’ 
 po(R, R’, f)= pi(R’)g(R, R’, ), 


(6) 


.where it is understood that & has the value 


unity when not indicated explicitly. Properly, 
pi(R’, £) should appear in Eqs. (7), but it may 
be shown that p;(R’, &) differs from p;(R’, 1) by 
terms of negligible order when £ does not refer 
to the molecule fixed at R’. We shall refer to the 
distribution function g(R, R’) as the correlation 
function, since it determines the correlation, 
(fi(R) fo(R’) w—(fi(R))a(fo(R’))w, of any two 
functions f,(R) and f2(R’) of the positions of 
two molecules. The correlation, 


1 v v . 
1J 2/Av — (J 1) av. J 2) Av = i(R 2(R’ 
(Fufdn—Adakfern=— ff ful) fa0R) 


XLg(R, R’) —1]ei(R)pi(R’)dvdv’, 


(8) 


would vanish if g(R, R’) were identically unity 


5 The distribution functions are similar to those defined 
by i Mayer and E. Montroll, J. Chem. Phys. 9, 2 
1941). 





STATISTICAL MECHANICS OF FUSION 


everywhere. In a liquid phase for which p;(R) 
is unity everywhere, the correlation function 
g(R, R’) reduces to the radial distribution func- 
tion go(|R’—R|) of the liquid, which has been 


extensively investigated from the experimental 


standpoint by x-ray scattering. By virtue of 
their definitions, the several distribution func- 
tions satisfy the following conditions, 


1 v 
J o(R, £)dv=1, 


v 


— 
f e(R, R’, {)p(R’do’=1, (9) 


v 
pi(R,0)=1; g(R, R’, 0)=1; 
()) po(R, R’, 0) = i(R’). 
We shall now undertake the formulation of 
an integral equation for p;(R). Differentiation 
of both sides of Eq. (5) with respect to ~ and 


reference to Eqs. (3), (4) and (6) lead to the 
result, 


d log pi(R, £) 


N 
= — MV R, Av V Av_ty 
pon (R, £))w—(V(E))w] 


1 v 
(OV, D)w=- f V(|R’—R]) 
Uv 
X “ po(R, R’, £)dv’, 


(VOm=— f fviR-Rp 


X p2(R, R’, &)dudv’. 


(10) 


Integration of Eq. (10) and substitution from 
Eq. (7) yields 


log \p,(R) = f K(R, R’)p:(R’)dv’, 


K(R, R’) = —(N/vkT) V(|R’—R|)L(R, R’), 
(11) 


L(R, R’) = f o(R, R’, dé, 


log = —(N/vkT) f (V (8) ad é. 


The constant \ is fixed by the normalization 
condition, Eq. (9). If the correlation function 
g(R, R’, £) is known for all values of & between 
zero and unity, L(R, R’) can be calculated and 
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the kernel K(R,R’) of the integral equation 
constructed. A formal investigation of the 
properties of the correlation function leads to an 
integral equation for this function, similar in 
form to Eq. (11), with a kernel involving p;(R) 
and in addition the correlation function for sets 
of three molecules. Since the simultaneous 
integral equations obtained by a continuation 
of this process are extremely complicated and 
their solution prohibitively difficult, we shall 
base the further development of the present 
theory on approximations to g(R,R’, £) sug- 
gested by certain reasonable physical hypotheses. 

The long range order characteristic of crystal- 
line structure is determined by the distribution 
function p:(R). The local order imposed by a 
molecule fixed at a point R on its environment 
is determined by the correlation function 
g(R, R’), which may conveniently be regarded 
as a function of R and R’—R. The average 
density at a point R’—R from a molecule fixed 
at R is proportional to g(R, R’—R)p»;(R’). Thus 
the dislocation of the lattice by a molecule fixed 
at an arbitrary point R is determined by the 
correlation function. On physical grounds, one 
is led to attribute the following properties to the 
correlation function 


g(R, R’—R)=1, 


lim 
R’—R!-— © 


(12) 
lim 


plim se R’—R, §)=0; &>0. 


The first of these relations implies that the 
correlation in position between two molecules 
vanishes at large intermolecular distances. The 
second implies the existence of large repulsive 
forces at small distances, excluding molecular 
interpenetration. The first of the relations (12) 
requires careful examination. In a three-dimen- 
sional system, it is true in the liquid phase, as 
we know from the properties of the radial 
distribution function to which it reduces. It can 
also be shown to be true in the well-ordered 
crystal with harmonic lattice vibrations. How- 
ever, it is not valid for the linear chain or for 
the surface grating, in which the correlation in 
position between two molecules remains finite at 
infinite separation.* It is interesting to remark 


* We are indebted to Professor P. Debye and Professor 
J. E. Mayer for illuminating discussions on this point. 
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that the approximation of g(R, R’—R) by unity 
everywhere is equivalent to assuming an Ein- 
stein crystal. 

From another point of view, the correlation 
function may be regarded as a_ functional 
g(R’, R,[]) of the distribution function p; in 
the sense of Volterra.* Though we shall not make 
extensive use of this property, we remark again 
that for the disordered phase, p;(R) =1, 


g(R’, R, [1 ]) =go(|R’—R]), (13) 


where go(|R’—R}) is the liquid radial distribu- 
tion function, which depends only on the 
distance |R’—R|. We shall base the further 
development of our theory on the approximation 
of g(R’, R, [p]) by g(R, R’, [1]), that is by the 
radial distribution function go(|R’—R]|) of the 
disordered liquid phase, perhaps metastable. 
The existence of undercooled liquid phases makes 
this approximation physically meaningful. Its 
adequacy from a quantitative standpoint is 
difficult to estimate by theoretical criteria, but 
can be judged a posteriori in the applications of 
the theory which are to follow. With this 
approximation, the integral equation (11) be- 
comes 


log Xo.(R) = f K.(|R’—R|)pi(R’)dv’, 


Ko(|R’—R|) = —(N/vkT) 
x V(|R’—R|)Lo(|R’—R|). (14) 


La(|R’—R|)= f go(|R'—RI, Bde 


where due to the short range character of 
V(|R—R’|) the integral may be extended to all 
of space, except in the immediate vicinity of 
the boundaries of the region v. 

We shall be interested in those solutions of 
Eq. (14) with the period of a specified space 
lattice. Among the possible lattices and distribu- 
tion functions, the one corresponding to a mini- 
mum value of the work content at constant 
temperature and volume is to be selected as 
appropriate to thermodynamic equilibrium. The 
periodic distribution function and its logarithm 

6 V. Volterra, Theory of Functionals and of Integral and 


Integro-Differential Equations (Blackie and Son, Ltd., 
1930). 
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may be expanded in the Fourier series, 


+o 
pi(R)= > s(h) exp (2rth-R), 
hi, he, hg 


+0 
log \pi(R)= > gq(h) exp (2mth-R), (15) 
hi, h2,h3 


h=/,bi+h2b2+Asbs, 


where fy, he, h3 are integers and b,, be, bs are 
the basic vectors of the reciprocal lattice of the 
assumed structure. Substitution of the Fourier 
series (15) into the integral equation (14) yields 
the relations 


q(h) =a(h)-s(h), 
ah) = (2/h) f rKy(r) sin (2rhr)dr, (16) 
h=|h|. 


For periodic solutions, the normalization condi- 
tion, Eq. (9) becomes. 


1 A 
-f p(R)dv=1, (17) 
A 


where the integration extends over the ele- 
mentary cell A of the lattice. Substitution of the 
second of the series (15) into Eq. (17) yieldst 


A 
A= (1/a) f exp [> gq, exp (2rth-R) ]dv, (18) 
h 


which determines A as a function of the variables 
gn. Calculation of the Fourier coefficient s, from 
the second of the series (15) gives 


1 4 
n=— f exp (—2zih-R) 
AA 
Xexp [> gn exp (2rth’-R) dv 
h’ 
(19) 
s,=0 log /dqn*, 
where q,* is the conjugate of g,, equal to ¢-x, 


since both p;(R) and its logarithm must be real. 
Combination of Eqs. (16) and (19) gives 


qr = anGn(Qo, rey ), 


Gi(qo ++) =9 log X/0qx*, 


+ Note—For simplicity in notation a(h), g(h), and s(h) 
will occasionally be written as an, gn and Sn. 


(20) 
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an infinite set of transcendental equations for 
the Fourier coefficients g, of the logarithm of 
the distribution function. Thus the determina- 
tion of real periodic solutions of the integral 
equation (14) is equivalent to determining the 
solutions of the set of transcendental equations 
(20) as functions of the transforms a(h) of the 
kernal Ko(r). From Eqs. (20) and the normal- 
ization condition (17) we remark that 


qo= Qo, So=1,. (21) 


where we write go and so for gooo and Sooo. By 
inspection, one finds that Eqs. (20) always have 
the solutions g,=0 for h>0, according to which 
pi(R) has the constant value unity everywhere, 
corresponding to a disordered phase. We shall 
not undertake at the present time the difficult 
mathematical investigation of the infinite set of 
equations (20). However, when the transforms 
a, vanish for h greater than some finite ho, the 
equations reduce to a finite set of transcendental 
equations. It may then be shown that for speci- 
fied lattices, non-trivial periodic solutions p,(R), 
with non-vanishing Fourier coefficients s, for 
h>0O, exist for sufficiently large values of the 
transforms a, and correspond to lower work 
content than the solution p;(R)=1. Since the 
transforms a, are functions of temperature and 
volume, the Fourier coefficients s, of p,(R) 
become discontinuous functions of volume, all 
except Sooo Vanishing when the volume exceeds a 
critical value v* at a specified temperature. We 
shall postpone the detailed discussion of the 
approximate solutions of Eqs. (20) to a later 
section of the article. 


THE THERMODYNAMIC “FUNCTIONS AND 
HETEROGENEOUS EQUILIBRIUM 


We shall undertake the formulation of ex- 
pressions for the thermodynamic functions of 
the system in terms of the Fourier coefficients 
s, of the distribution function p,(R), under the 
restraint of homogeneity according to which the 
periodic solution of the structure of minimum 
work content at the specified average density 
and temperature is assumed valid throughout 
the region v occupied by the system. The lowest 
work content may actually be realized with 
different periodic solutions in different domains 
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of the region v, corresponding to the coexistence 
of two or at most three distinct phases. Although 
this situation can be treated explicitly by 
minimizing the work content with respect to 
the masses of the several homogeneous parts, 
it is more convenient to employ the thermo- 
dynamic criteria of heterogeneous equilibrium 
in conjunction with the theoretical expressions 
for the thermodynamic functions calculated 
with the restraint of homogeneity. 

We refer to Eq. (18) of the article ‘Statistical 
Mechanics of Fluid Mixtures’ by one of us.* 
From this equation and the last of Eqs. (11) 
of the present article, we find for the molar 
Gibbs free energy F, identical with the chemical 
potential in a one-component system, 


F/RT= —log w+ F(T) RT, 
(2amkT)} 


F°(T)/RT= —log —- 
Nh* 


where v is the molar volume, N is Avogadro's 
number, m the molecular mass and / Planck’s 
constant. The parameter \ is specified as a 
function of the Fourier coefficients g, by Eq. (18) 
of the present article. The mean energy E of 
the system is 


E/RT=3+(N/2kT)V, 


V=(1 of fir R’—R!)g(R’, R) 


X pi(R)pi(R’)dedv’, (23) 


where V is the average potential energy of a 
molecular pair. Introducing the Fourier series, 
Eq. (15), into the integral of Eq. (23) and 
approximating the correlation function by the 
radial distribution function go(|R’—R!}), we 
obtain 
E/RT= —(3) } Br! Sh 24%, 
h 

(24) 
B.= —(2N/ohkT) f rV(r)go(r) sin (2rhr)dr. 

0 


If we define 6,(£) as 


Bu(e)=—(2N/ohkT) fr V(rga(, ) sin (rind. 
a 20 
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J. 


we remark from Eqs. (14) and (16) that the 
transforms a, and §,(£) are related in the 
following manner 


a= f B,(é)dé. (26) 


In order to obtain the equation of state, we 
differentiate Eq. (22) with respect to volume at 
constant temperature and refer to Eqs. (16) 
and (18) 


v /op\ 1 1 1 FAQr2 

-—(— --+-5-(—) (27) 
RT\007 7 v 2r arr dvJ or 

We shall assume that the stable solutions of 

Eqs. (20) are continuous functions of the volume 

except at a point v* and that g,=0 for v>v* 

and h>0O. Integration of Eq. (27) under these 


conditions with a minor transformation based 
upon partial integration yields 


pv/RT= 1 —AQ*+3 Eu | sa]? 
h 


J O(ayv’) ” ao 
+e f (1/v’)? : | Sa | ‘w| -of —d'" (28) 
* y’? 


, 
v v rd 





Q=log  —4 Yo’an|sn|/2; A’ =A exp (—ao), 
h 
AQ* = 2(v* —0) — 2(v* +0), 


where the prime denotes the omission of the 
term for 4=0 in the sums and the discontinuity 
AQ in pu/RT is determined by the requirement 
that the work content F—pv be continuous at 
v*. For volumes not differing greatly from 2%, 
the integrals 


v ra) / 
of (4/of |g, |2ae’ 





Qh 
dv’ 
can be neglected without great error since they 
can be shown to be bounded in absolute value 
by [an—(v*/v)a,*}|\s,|2 if av is monotonic 
function of the volume. 

If we denote by F’, A’ and E’ the excesses 
(F—Fy), (A—Ao), and (E—£,) of the free 
energy, work content, and energy over their 
disordered values Fy, Ao and Ep at the specified 
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temperature and volume, we may write from 
Eqs. (22), (24) and (28) 


F’/RT=—log \’, 
A’'/RT=—log VN +4 Dd’ an|sn|?+A2*, (29) 
h 


E’/RT=—}3 L Br | S| *, 


Sam AQ*, 





p'v/RT= —} + ty Qh Sh 
h 


where the last two formulas involve the neglect 
of the volume integrals between v* and v in 
Eq. (28). 

The conditions for heterogeneous equilibrium 
between the liquid phase (disordered) and the 
crystalline phase (ordered), equality of the 
chemical potentials and pressures in the two 
phases, combined with Eqs. (29) lead to the 
relations 


log = —AF,/RT, pb’ =Apo, 
AF,=F,(T, v,) —F,(T, Vs); (30) 
Apo=poT, v1) — pol T, 2s), 


where Fy and pp are the free energy and pressure 
of the undercooled disordered phase. Since }’ 
and p’ are functions of the g,, which depend 
upon temperature and volume through the 
transforms a,, Eqs. (30) provide two relations 
between v; and v, at a specified temperature, 
from which, in principle, the liquid and solid 
volumes may be determined, and from them the 
other fusion parameters calculated. The incre- 
ments in the thermodynamic functions of the 
disordered phase may be approximated in the 
following manner, 


AFy= — Av/ko, 
Apo= —(1/ko) log v/v, 
Tyo Vs Vi Av 
ABytpav=—"+""[ og ~——|, (31) 
Ko Ko Vs Vs 
Av=v1—2,, 


if the compressibility xo and the thermal expan- 
sion coefficient yo are treated as constant in 
the volume interval Av. With the aid of Eqs. 
(29) and (31), Eqs. (30) may be put into the 
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somewhat more convenient form, 
log \’ = Av/xoRT, 


Vs jae V1 
Q=A0N*+— ——log (— 
KoRT Vs 


9) 


Us 


Q=log -3 hy an | Sn |? 
h , 


The entropy of fusion, AS, equal to AH/T, is 
calculated from Eqs. (29) and (31) to be 


AS/R=3 dD’ Ba! sn}? +ASo/R, 
h 


Av Vs Vv) Av 
— log ( — J-—— : 
0 Kok T Vs Us 


K 


(33) 


Yo 
AS, /R= 


Eqs. (32) and (33) will be employed to determine 
the fusion parameters Av and AS in a later 
section of the article. 


THE TRANSFORMS a, AND 6, 


The transform §; defined in Eq. (24) may be 
computed by numerical integration from the 
potential energy V(r) of a molecular pair and 
the experimental radial distribution function 
go(r), extrapolated if necessary to the under- 
cooled liquid region. The calculation of 8,(&) 
and a, is considerably more difficult, since it 
demands a knowledge of the radial distribution 
function go(r, €) for intermediate degrees of 
intermolecular coupling not realized in the actual 
system. However, certain approximate results 
can be stated without entering into a detailed 
theory of the radial distribution function. 

It is convenient to employ the impenetrable 
sphere approximation to intermolecular repul- 
sion, at small distances, according to which the 
molecules are assigned a diameter ado. This may 
be accomplished by introducing a potential of 
intermolecular force, 

Vir) = Vo; 


0<r<ao, (34) 


where Vo is a positive constant which is allowed 
to approach infinity in the partition function 
and the distribution functions derived from it. 
If we divide the interval of integration in Eq. 


25) defining 8,(£), we may write 


Br (E) = roBn (E) +B (8), 


—2N 7” 
By (£) — f rgo(r, &) sin (2xhr)dr, (35) 
Un 0 
—2N fr” 
6, () =—— rV(r)go(r, &) sin (2xhr)dr, 
vhkT 


2 


where 79 is equal to Vo/kT and V(r) is the 
potential of intermolecular attraction in the 
interval ag? r<«. Examination of the phase 
integrals defining the correlation function, here 
approximated by the radial distribution function, 
shows that it depends upon the coupling pa- 
rameter £ in the form go(r, 70, £/7). We shall 
need to consider the following limits 


lim golr, Ero, E/T) =o(r, £/T); §E>0 


lim 


TID 


(36) 
go(r, t, t/T0) =go0(?, T). 


The function go(r, £/T) corresponds to full re- 
pulsive coupling and partial attractive coupling 
of one molecule of the pair with the other 
molecules of the system, whereas go(7, 7) corre- 
sponds to partial repulsive coupling and vanish- 
ing attractive coupling. Substitution of the 
6.(&) of Eq. (35) into Eq. (26), change of variable 
in the first integral and passage to the limit 
Vo yields 


_— ) 
an =a +an,™, 


—2N 
a,O= 


f rL©(r) sin (2xhr)dr, 
vh 0 


e 4) 


_rV(r)L™(r) sin (2xhr)dr, 
vhkT Jay 
(37) 


ay) =_—_— 


Lo(n)= f go(r, 7)dr, 
0 


1 
Lin(r)= flr, &/ Tak. 
If go(r,&/T) for r>aq is assumed to be an 
analytic function of &/7T, go(r,&/T) differs 
from the radial distribution function go(r) for 
complete intermolecular coupling, £=1 by terms 
of order 1/7 and we may write, remembering 
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Eq. (26), 


a," = B,{1+0(1/T)}, (38) 


where the neglected terms of order 1/7 represent 
the effect of intermolecular attraction on the 
radial distribution function. Since 6, is of the 
order 1/T, we may write 


an, =a, +8,+0(1/T*), 


where a, depends upon both volume and 
temperature, but reduces to a function of volume 
only, for a system of impenetrable spheres with 
vanishing intermolecular attraction. 

In order to get a clearer idea of the properties 
of the function L°(r) from which the transform 
a, is calculated, we shall examine the thermo- 
dynamic functions of the disordered phase alone. 
From Eqs. (22) and (28) we may write for the 
work content 


A,/RT = —log (e)+0 f (ao/s*)do+ FY/RT. (40) 


(39) 


If on the other hand the work content is ex- 
panded in the moments of the potential of 
intermolecular attraction, we obtain 


Ao/RT = —log (evs) — Bo/2+O0(1 eas 
vp=(Zy)", 

where Zy° is the phase integral appropriate to 

impenetrable spheres with no attractive forces, 

and vy is a free volume of the type introduced 

by Eyring’ or by Lennard-Jones and Devonshire® 


in their theories of the liquid state. Comparison 
of Eqs. (40) and (41) yields 


d log (v/v) 
d log v 


1 O(Bov) . 1 
-1(28°) s5.40(4) 
2 ov T T? 


From Eqs. (14), (15) and (39), we obtain 
4nN % 


v 0 


ay= —log (v/vs) — 


v 
r?L (r)dr=log (-) 


vs 
d log (vs/v) 1 (Bor) 
+- 


0d logv 2 ov 


7H. Eyring and J. O. Hirschfelder, J. Phys. Chem. 41, 249 
(1937). 

8 J. E. Lennard-Jones and A. F. Devonshire, Proc. Roy. 
Soc. A163, 53 (1937). 


(43) 
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If L°(r) is approximated by a constant in the 
sphere do, the transform a, becomes 


a(v, T) 
a, = 


(2rhao) 2 


o(v, T)= 3| log (-) + 
Us 


As we shall later see, the term a, in the 
transform a, leads to the prediction that a 
system of hard spheres without attraction must 
crystallize at sufficiently small volumes. 


sin 2rhao 
cos 2rhay—- 


2rhay 
(44) 
0 log (v;/v) 1 9(Bov) | 
2 ov J” 





0 log v 


APPROXIMATIONS FOR THE FACE-CENTERED 
Cusic LATTICE 


We shall investigate the solutions of Eq. (14) 
for a cubic lattice with a unit cell of volume 
A=4v/N, having an average content of four 
atoms. The lattice vectors ai, a2, a3 are orthog- 
onal and of the same length a equal to A}. The 
reciprocal vectors b;, be, bs are also orthogonal 
and of equal magnitude 1/a. We denote by x, y, 
and z the components R-b;, R-be, and R-b;. 
If we approximate the Fourier transform a, by 
zero for values of h greater than or equal to 2/a, 
we conclude from Eqs. (20) that all g, vanish 
for h>v3/a, and the infinite set of transcendental 
equations reduces to a finite set determining the 
twenty-six variables g, in the distribution 
function* 

—>+1 


log N’pi(x, y,2)= Dd’ g(h, ho, hs) 
hi, he, hg 


=-1 


Xexp [2ri(hixtheyt+hsz) ], (45) 


the constant term in the Fourier series being ab- 
sorbed in the parameter )’ defined as \ exp (— 4). 
Since the transforms a, depend only on the 
absolute value of h, it is evident that Eqs. (20) 
have solutions depending only on the absolute 
value of h, in which case the twenty-six variables 
q(h) reduce to three q’’, g’ and q corresponding 
to values of h of 1/a, V2/a and v3/a, respectively. 
Since the origin to which p;(x, y, 2) is referred is 
arbitrary, from one set of solutions qg, other scts 
g, exp 27th: Ro), with Ro an arbitrary constant 


* The vanishing of the coefficient g, for h > 2/a is duc to 
the neglect of the a, and does not imply the vanishing of the 
corresponding Sp. 
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TABLE I. The transcendental equation. 





| log »’ 

| 0.506 =log ’(q*) 
0.5738 
0.7819 
1.0203 
1.287 
1.578 
1.892 
2.223 
2.572 


0.9730 =a* 0.362 =s* 
0.9734 
0.9783 
0.9897 
1.0082 
1.0325 
1.0633 
1.0991 
1.1400 
1.1848 .7385 2.936 
1.232 .7609 


3.312 
| .7794 3.696 


84q 
2.82 =8q* 


3.00 
3.50 
4.00 
4.50 
5.00 
5.50 
6.00 
6.50 





0.3853 
4472 
.5052 
.5579 
.6053 
.6466 
.6824 
7127 


7.00 
7.50 
8.00 


vector, may be constructed. If we place the 
origin at the point (0,0,0) in any cell and 
employ the solutions q’’, g’ and q, the distribution 
function, Eq. (45), becomes 


1 
pi(x, y; Z) =i ¥, 2) f'(x, y; 2) f(x, y; z), 


f(x, y, 2) =exp [6q’"(cos 2rx 
+cos 2ry+cos 272) ], 
(46) 

f'(x, y, 2) =exp [12q'(cos 27x cos 2ry 

+cos 2ry cos 272 
+cos 272 cos 27x) ], 

f(x, y, 2) =exp [8q(cos 27x cos 27y cos 272) J. 
The factor f(x, y,z) has the symmetry of a 
simple cubic structure, f’(x, y, z) that of a body- 
centered structure, and f(x, y, z) that of a face- 


centered structure. The transcendental equations 
(20) reduce to 


a’ A log X’ 


'—.. - 
12 oq 


’ 


ad log 


ee 
8 dq 


= J | J | J ‘ite paften9) 


X f(x, y, z)dxdydz, 


where a’, a’ and a are the values of the Fourier 
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transform a, for 1/a, v2/a, and v3/a. Since it 
may be verified that 0 log \’/dq’’ and @ log \’/dq’ 
vanish for all values of g when g’ =q’’ =0, special 
solutions of Eqs. (48) are obtained with q’=q"’ 
=(0. Corresponding to these solutions, is the 
distribution function 


pi(x, y, 2) =(1/d’) exp [8q(cos 2rx 


Xcos 2ry cos 272) |, 


1 1 1 
= f f i) exp [8g cos 27x cos 27y 
0 0 0 


cos 2rz |dxdydz, 


(48) 


which, for positive g, has congruent maxima at 
the points (0,0,0), (3,4,0), (4,0,3), and 
(0, 3, 3), that is at the lattice points of a face- 
centered structure of cell A. Properly, it should 
be demonstrated that the face-centered distribu- 
tion, Eq. (48), corresponds to the lowest work 
content consistent with the specified temperature 
and volume, before accepting it as the equi- 
librium distribution. For the present, we shall 
assume this to be true without proof. A proof 
would have to be based upon the specific form 
of the potential of intermolecular forces from 
which the transforms a, are calculated. For the 
distribution function, Eq. (48), the transcen- 
dental equations (47) reduce to a single equation 


q=aG(q), 
G(q) = 5 (0 log d’ 0g) q”'=1'=0- 


(49) 


Development of the exponential in the expression 
for \’, Eq. (48) leads to the following expressions 


for \’ and G(q) 


; « [(2n)!}? : 

\’(q) = a Tah q°", 

q 2 (2n+2)[(2n+2)! ]*g* , 

Gq) =— > ———._ (50) 
(q) > (nb)! ( 


8r’ n=0 
from which \’ and G may be computed as 
functions of g. For a less than 0.973,’ it is found 
that the only real solution of Eq. (49) is g=s=0. 
For @ greater than 0.973, there are two non- 
vanishing positive solutions, the larger of which 


n 


® Note-—Due to certain modifications in the theory, this 
critical value of @ differs slightly from the one stated in our 
preliminary report, J. G. Kirkwood and E. Monroe, J. 
Chem. Phys. 8, 845 (1940). 
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TABLE II. Reduced radial distribution function for .liquid 
argon from data of Eisenstein and Gingrich. 








4xx?g0(x) 4ax2go(x) 


1.96 17.88 
16.64 25.73 
21.86 32.30 
14.30 38.10 
16.21 44.27 
18.14 














represents the equilibrium distribution, since it 
may be shown to correspond to the least work 
content. Corresponding to each positive solution 
q of Eq. (49), there is also a negative solution 
—q. Since the negative solutions merely imply a 
change in origin from the point (0, 0, 0) to the 
point (4, 3, 3), we need only consider the positive 
solutions. In Table I, the stable solution gq of 
Eq. (49) and certain important functions of q 
are presented as functions of the parameter a. 
Subject to the approximations which have 
been described, Eqs. (32), expressing the condi- 
tions for equilibrium between the liquid phase 
and a face-centered crystalline phase, become 


log \’(q) =Av/KoRT, 
og \’(q) /Ko (51) 


Av Vv) 
Fae 
Us Vs 


2 9) 


q 0s 
log \’(q) -—— = 0.0043 + 


a Ko 
and the entropy of fusion, Eq. (33), becomes 
AS/R=4Bs?+ASo/R, (52) 


where g, 8 and s denote the values of gn, 8, and sp 
for h=v3/a at v., the 6, for h22/a being 
approximated by zero as in the case of the ay. 
From Eq. (24), we may write 


4 «o 
a eae f sVivaledra) in (i6an)dx, (53) 
° V3RT Jo 
x=r/ro; ro=a/v2, 


where 7p is the distance between the lattice 
points of the face-centered structure of cell A 
equal to 4v/N. An alternative expression for 


AS/R is 
AS/R=48Bs?—ABo/2+pdr, 
ABo = Bo(v1) — Bolas), 


where Bo is the value of the transform #) for 


(54) 
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h=0. With the aid of Eq. (24), it may be written 
in the form, 


—4rv2 fr” ee (55) 
Bo=———  ]_ x°V(rox)go(rox)dx. (55) 
0 eT q O\F 0 


Were a known as a function of volume and 
temperature, Eqs. (49) and (51) would provide 
three relations from which to determine gq, v, 
and v, at a specified temperature. From the 
values g and Av so determined and values of { 
and 6, AS could then be calculated with the aid 
of Eq. (53), or Eq. (54). Although it is possible 
to calculate Bo and 8, deficiencies in the theory 
of the liquid state allow only rough estimates 
of the parameter a. Using a free volume of the 
Eyring type it is found with the aid of Eqs. (39) 
and (44) that a@ is of the order of magnitude 
unity at volumes intermediate between v, and 2, 
in the case of argon, and thus approximates the 
critical value 0.973 in the fusion range. In the 
absence of adequate theoretical estimates of a, 
we can, however, calculate the fusion parameters 
Av and AS if the volume of the liquid v; is taken 
from experiment, by solving Eqs. (49) and (51) 
simultaneously for g, a and v,. An estimate of 
the entropy of fusion alone may be obtained if 
both v; and Av are taken from experiment and 
Eqs. (49) and the first of Eqs. (51) solved 
simultaneously for g and a. Both methods will 
be employed in the next section for the calcula- 
tion of the fusion parameters of argon. 


THE FusION PARAMETERS OF ARGON 


The liquid-solid transition in argon provides 
a favorable case for testing thedries of fusion. 
Not only has the transition been extensively 
studied from the experimental side, but also the 
intermolecular force between argon atoms is of 
a simple type and has been investigated in some 


TABLE III. The integrals Bo and 8. 








v BoT 
(CC/MOLE) 


21.3 
22.6 
23.9 
25.2 
26.6 
28.1 
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detail. In order to calculate the fusion parameters 
of argon on the basis of the present theory, 
information concerning the potential energy 
V(r) and also of the radial distribution function 
go(r) is necessary. For V(r), the potential energy 
of a pair of argon atoms, we use the familiar two 
power expression due to Lennard-Jones, as 
modified by Rushbrooke"” in accordance with 
the calculations of Corner." In our notation, 
this potential is 


Vi 3.27310" 6.487 10° 


k (ax) "1-4 ig (ax)® 





(56) 


x=r/to; ro=a/v2 


with the parameter a expressed in Angstrom 
units. Eisenstein and Gingrich” have determined 
the radial distribution function of liquid argon 
at 90°K under a pressure of about 5 atmos. from 
measurement of the angular distribution of 
intensity in x-ray scattering. In view of the 
small compressibility of liquid argon, the density 
of the liquid under the conditions of their 
measurements may be taken as 1.374 g/cc or 
0.02085 atoms/A* as measured by Mathias, 
Onnes and Crommelin® for liquid argon under 
its Own vapor pressure at 90°K. In order to 
extrapolate the Eisenstein-Gingrich function to 
other densities, we assume that the radial 
distribution function is a function only of the 
reduced radial distance r/ro without explicit 
dependence on density. Although this assump- 
tion lacks theoretical justification, we may 


TABLE IV. Fusion parameters for argon. 








AS/R AS/R AS/R 
CALC. CALC. CALC. 
(Eg. (52) (Ea. (52) (Ea. (54) 
WITH Av WITH Av WITH Av 

CC/MOLE CC/MOLE (CALC.) OBS.) OBS.) 


83.9 1.68 3.53 3.25 1.74 1.78 2.02. 
119.7 110 188 62 .70 1.50 1.63 
183.2 ‘71. «92 33—Sts«CB’ ‘97 «1.06 


Av Av 
AS/R OBS. CALC. 


KELVIN OBS. 











i G. S. Rushbrooke, Proc. Roy. Soc. Edinburgh 60, 182 
(1940). 

1]. Corner, Trans. Faraday Soc. 35, 711 (1939). 

2 A. Eisenstein and N. S. Gingrich, Phys. Rev. 58, 307 
(1940). We are indebted to Dr. Gingrich for a private 
communication placing at our disposal certain data, some 
of which are used in Table II. The radial distribution 
function employed here does not represent the optimum 
function obtainable from their data (N. S. Gingrich, Phys. 
Rev. 59, 290 (1941)), but small variations do not greatly 
affect the transforms Bo and 8. 

‘’ Mathias, Onnes and Crommelin, Proc. Acad. Sci. 
Amst. 15, 667 (1913). 
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TABLE V. Observed properties of liquid argon at the 
melting point. 


Tm 
DEREES K 


83.9 1 
119.7 1650 
183.2 5350 


xo 
(KG/CM?) ~! 


2.34 10-4 
6.1 X10~5 
2.1 K10-5 


exjeness (pacanss) “1 
0.709 
.661 
.631 


Pm 
KG/CM? 








Sources: At high pressures: v1 from Bridgman's" ‘‘best values for the 
volume. .of the amorphous phase. .on the melting curve’’; yo by extra- 
polation from his Table IV; xo from yo, his vi and dfm/d7Tm, using the 
relation, xo =yo(d7m/dpm) —(1/vi) (dvi/d pm). At one atmosphere pressure: 
vt, yo by extrapolation of liquid density data; «xo from T, vi, yo and 
(Cp—Cv) for the liquid (Cp —C» =4.8 cal./mol/degree—Eucken and 
Hauck, Zeits. f. physik Chemie 134, 161 (1928)). 


expect it to be reasonably adequate for extrapo- 
lation over small density intervals. The radial 
distribution function go(x) calculated from the 
data of Eisenstein and Gingrich is presented as 
a function of the reduced interatomic distance 
x in Table II. The integrals’ 8 and Bo of Eqs. 
(53) and (55) calculated by numerical integration 
with the interatomic potential (56) and the 
data of Table II are presented as functions of 
the molar volume v and the lattice parameter a 
in Table III. Since the energy of vaporization 
AE, of the disordered phase is equal to ByRT/2, 
an experimental check on the calculated value 
of Bo is possible. For liquid argon at 90°, AE, /RT 
has the observed value 7.14, with which the 
calculated value of 89/2 equal to 7.17 compares 
favorably. 

The theoretical fusion parameters of argon, 
AS and Av, are presented in Table IV, calculated 
with the aid of the best available data (Table V) 
on the density, compressibility and thermal 
expansion coefficient of the liquid and compared 
with the experimental values at several pressures. 
The fusion parameters were calculated in three 
ways. Both AS and Av were calculated with the 
aid of Eqs. (49), (51) and (52). In addition AS 
alone was calculated with use of the observed Av, 
Eq. (49), the first of Eqs. (51) and with each of 
Eqs. (52) and (54). At one atmosphere, com- 
parison is made with the data of Clusius and 
Weigand" and at the higher pressures with the 
data of Bridgman.’* At one atmosphere the 

44 Note-——We may anticipate that the radial distribution 
function of K. Lark-Horovitz and E, P. Miller, Nature 146, 
459 (1940), or the interatomic potential of O. K. Rice, J. 
Am. Chem. Soc. 63, 3 (1941) would lead to somewhat higher 
values of the integrals 6 and Bo. 

15 K, Clusius, Zeits. f. physik. Chemie B31, 457 (1936); 


K. Clusius and K. Weigand, ibid. B46, 1 (1940). 
146 P, W. Bridgman, Proc. Am. Acad. Sci. 70, 1 (1935). 











TABLE VI. Reduced melting parameters. 








Sus- Te Ve 





STANCE DEGREES CC/MOLE T»m/T- v1/%e v/Ve AS/R 
Ne 44.4 41.7 0.553 0.388 0.0513 1.64 
A 150.6 75.2 .556 .376 .0470 1.68 
Kr 209.4 92.1 554 371 .0487 = 1.69 
Xe 289.7 113.7 557 375 .0493 1.71 
Ne 126.0 90.1 .502 358 .0286 1.24 
CH, 190.6 99 A475 .340 0272 = =1.37 
CCl, 556.3 275.9 A451 333 .0144 1.16 
A (CALC.) .0432 1.78 








References: critical data:—except Kr—lInt. Crit. Tab. (Pickering), 
Kr (reference a); latent heats and 7m: Ne, A, Kr—(reference 15; 
Xe (reference b) ; Ne (reference c); CCl4 (reference d) ; CH4(reference e) ; 
vt: A (reference 13); Ne (reference f); Kr (reference a); Ne Int. Crit. 
Tab.; Xe CHa (see reference 15); CCla, Int. Crit. Tab. (extrapolated) ; 
Av: A, Kr, Xe, CHa (reference 15), Ne (reference 15), Ne (reference g), 
CCky, Int. Crit. Tab. : 

® Mathias, Crommelin and Meihuizen, Comptes rendus 204, 630 
(1937), or Proc. Acad. Sci. Amst. 40, 259 (1937). 

> K. Clusius and L. Riccoboni, Zeits. f. physik Chemie B38, 81 (1937). 
(933) F. Giauque and J. O. Clayton, J. Am. Chem. Soc. 55, 4875 

). 

4D. R. Stull, J. Am. Chem. Soc. 59, 2726 (1937). 

e K. Clusius, Zeits. f. physik. Chemie B3, 41 (1929). 

f Mathias, Crommelin and Onnes, Ann. d. Physik 19, 231 (1923). 

« W. H. Keesom and J. H. C. Lisman, Physica 1, 737 (1934). 


calculated fusion parameters agree remarkably 
well with the observed values. At the higher 
pressures the entropy increment agrees moder- 
ately well, the volume increment rather poorly, 
with the observed values. In view of the ap- 
proximations which it is necessary to introduce 
into the theory in order to arrive at explicit 
numerical results, no more than semi-quantita- 
tive agreement with.experiment is to be expected. 
We may therefore conclude that the experimental 
confirmation of our theory is adequate. 

It has long been a matter of speculation 
whether a critical point for fusion analogous to 
that for vaporization exists. Our theory predicts 
no critical point for the face-centered lattice if 
the parameter a is a monotone increasing 
function of density, unbounded at infinite 
density. This property is consistent with the 
approximate Eq. (44), since with increasing 
density the free volume tends to zero and a 
becomes infinite as close packing of rigid spherical 
molecules is approached. The basis for this 
conclusion is to be found in the fact that the 
liquid-solid transition is of the first order at 
constant volume. If at an arbitrary temperature 
a volume exists for which the parameter a 
attains its critical value, 0.973, the pressure 
jumps discontinuously to a lower value, and by 
the criterion of equal areas, the pressure-volume 


J. G. KIRKWOOD 





AND E. MONROE 





isotherm must possess a flat. It is to be remarked 
that no experimental evidence for the existence 
of a critical point for fusion exists, even in 
Bridgman’s extensive investigations of the 
liquid-solid transition at very high pressures. 

Some general theoretical conclusions con- 
cerning the fusion parameters may be based 
upon a law of corresponding states. As has been 
pointed out by Pitzer,!’ the partition function 
and all quantities derivable from it can be 
expressed as functions of the reduced tempera- 
ture 7/T., and reduced volume v/v, alone, if 
the potential of intermolecular force is a two 
parameter function of the distance, 7, and »v, 
being the critical temperature and critical 
volume. The two-power Lennard-Jones potential 
with fixed exponents is of this type. If such a 
potential is assumed, we would expect the 
transforms a, and £, of the present theory to 
depend only on the reduced temperature and 
reduced volume. The equations of the theory 
would thus involve only the reduced variables. 
As a consequence, the entropy increment AS 
and the reduced volume increment Av/v, of 
fusion at the triple point would have to be 
identical with those of argon at its triple point 
for all substances which have potentials of 
intermolecular force differing from that of argon 
only in the magnitudes of the coefficients of the 
two terms in the Lennard-Jones formula. If the 
triple-point pressures are of the order of magni- 
tude of one atmosphere, the comparison may 
equally well be made at one atmosphere, since 
the fusion parameters are not sensitive to 
changes in pressure of this order. Under these 
assumptions our calculations for argon should 
be applicable to the other rare gases. A summary 
of the reduced fusion data for these and for 
some other substances composed of approxi- 
mately spherical molecules is given in Table VI. 
Since all of the substances included in the table 
except nitrogen crystallize in face-centered cubic 
lattices, and all enjoy uninhibited molecular 
rotation at the melting point, the observed 
deviations from the law of corresponding states 
are doubtless attributable to departures from 
the assumed form of the potential of inter- 
molecular force. 

17K. S. Pitzer, J. Chem. Phys. 7, 583 (1939). 
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The Heat Capacity of Certain Halomethanes 
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A semi-empirical method of calculating the heat capacities of the halomethanes is presented. 
It is based on additive contributions of the bonds and angles making up the molecule. Values 
obtained for eighteen molecules at three temperatures by this method are compared with those 
calculated by statistical means. The average deviation is 9.9 percent. 





N recent years the halomethanes have become 

increasingly important in industry; yet there 
exist reliable thermodynamical data for only a 
few of them.' On the other hand, these com- 
pounds are being extensively studied by Raman 
and infra-red spectroscopy. In some cases where 
thermal data are available the two have been 
combined to yield complete thermodynamic 
functions in the temperature range in which the 
rigid rotator-simple harmonic oscillator approxi- 
mation is valid. Heat capacity and entropy 
have been calculated for some molecules for 
which no thermal data are available.2 However, 
there exist a considerable number of halo- 
methanes for which the fundamental vibrations 
are uncertain or unknown, and hence accurate 
calculations are not yet possible. Thus it seems 
worth while to investigate to what extent the 
existing data can be made to yield results of 
engineering value. The present investigation is 
concerned with a semi-empirical calculation of 
the heat capacity of certain halomethanes. 


THEORETICAL CONSIDERATIONS 


In all molecules of this type rotation has been 
excited to its classical value before room temper- 
ature has been reached; hence one is concerned 
only with the vibrational contribution. The 
vibration of a polyatomic molecule is a compli- 
cated process in which in each normal mode the 
motion of every atom is generally involved. 
Nevertheless, there exist certain modes in which 
the frequency is largely determined by the 


'R. C. Lord, Jr. and E. R. Blanchard, J. Chem. Phys. 4, 
707 (1936); Egan and Kemp, J. Am. Chem. Soc. 60, 2097 
ae Messerly and Aston, J. Am. Chem. Soc. 62, 886 

940). 

*See references in the review article by E. Bright 
Wilson, Jr., Chem. Rev. 27, 17 (1940). 


motion of only a few atoms. Thus in the molecule 
CXYZT there is a vibration, associated with 
the stretching of each valence bond, which 
remains nearly independent of the molecule in 
which the bond happens to be.* The remaining 
five vibrations result from the deformation of 
the six angles and are more or less characteristic 
of them. Thus one may write 


4 5 
OC =4R+D Cri +L C%), (1) 
1 1 


where the first term is the contribution of 
translation and rotation (plus R for C,—C,), 
and C%; and C%; are Einstein functions for the 
stretching and deformation vibrations, respec- 
tively. By choosing values of v; which are most 
nearly consistent with a number of molecules it 
is possible to approximate the second term with 
considerable accuracy. Table I gives a summary 
of the characteristic C—X vibrations in some 
halomethanes and the values chosen for each 
bond.* 4 

The deformation vibrations are, in general, 
more nearly functions of all the angles rather 
than any one, and hence the same kind of 
treatment cannot be followed. However, the 
third term may be satisfactorily approximated 
by expanding it into a summation over the six 
angles, weighted according to the total number 
of angles. Equation 1 then becomes 


4 5 
C*,=4R+>. ae p n;C%s;, (2) 
1 i 


3G. R. Leader, ‘‘Studies in molecular structure,’’ Ph.D. 
thesis, University of Minnesota (1940) and others. 

*Ta-You Wu, Vibrational Spectra and Structure of 
Polyatomic Molecules (China Science Corporation, Shanghai, 
1939). . 


527 





G. GLOCKLER AND W. 


F. EDGELL 


TABLE I. Characteristic C—X frequencies. 








C-H C-F 


C-Cl 





CHCl; 
CHCI1.F 
CHBr.F 3017 

CHBrCl, 3020 
v= 3000 


3020 cm™ 
3019 


CH3F 
CH.CIF 
CHC1,F 

CC1;F 

v= 1050 


CH;Cl 
CH:2BrCl 
CHBr:Cl 

CBr;Cl 
v=730 








where n; is the number of j angles and the 
summation is over all different angles. C%; is 
the Einstein function for the hypothetical char- 
acteristic deformation frequency 6;, associated 
with each angle, which will be evaluated later. 
Thus, by Eq. (2) the heat capacity of any 
halomethane is solely determined by the bonds 
and angles which make up the molecule. The 
success of this treatment is indicated by the 
values in Table III. 

This discussion does not apply to molecules 
of the type CX, since they have two stretching 
and two deformation vibrations separated widely 
enough to invalidate the above statements. This 
puts little restriction upon the general usefulness 
of Eq. (2). 

Bennewitz and Rossner® have successfully 
applied a similar treatment to molecules con- 
taining C, H and O in which a deformation 
frequency was assigned to each bond rather than 
angle. This type of approximation is valid 
whenever the forces resisting deformation arise 
primarily from the directed nature of the valence 
bond or whenever the bond appears only in 
angles of the same kind as in those molecules 
treated by Bennewitz and Rossner: The forces 
between the atoms not joined by valence bonds 
are of considerable importance in halomethanes 
with several halogen atoms. This is clearly 
shown by the fact that the normal coordinate 
treatment of halomethanes with three or more 
halogen atoms has been successful only when 
based on a central force field, i.e., assumption 
of forces only along the lines joining the atoms 
in the molecule. Thus it is not surprising to find 
that the modification presented above is more 
successful in treating molecules containing 
several halogen atoms. 


5 Bennewitz and Rossner, Zeits. f. physik. Chemie B39, 
126 (1938). . 


APPLICATION TO THE HALOMETHANES 


In order to evaluate the third term in Eq. (2) 
it is necessary to calculate the hypothetical 
deformation frequencies to be associated with 
each angle. In view of the meager experimental 
data available these constants have been ob- 
tained from heat capacities evaluated statisti- 
cally.*7 They are given in Table II, together 
with their contribution to heat capacities at 
several temperatures. For temperatures not 
included in Table II the contributions of the 
various frequencies may be found in tabulated 
values of the Einstein function.2 The standard 
state is the usual one of unit fugacity and zero 
pressure. Corrections to finite pressures have 
been adequately discussed.” 8 

The most useful application of Eq. (2) is to 
molecules for which no data at all are available. 
Thus one may calculate the heat capacity of 


TABLE II. Contribution of angles and bonds to heat capacity. 








CY or C'% 
373.1 


0.0738 
0.738 
1.738 
0.438 
1.503 
1.746 
0.825 
1.849 
1.622 
1.789 
0.0026 
0.595 
1.073 
1.284 


ANGLE OR 
BonD 


HCH 
HCCl 
CICCl 
HCF 
FCF 
FCCI 
HCBr 
BrCBr 
FCBr 
CICBr 
C-—H 
C-—F 
C-Cl 
C-—Br 





298.1 


0.0186 
0.451 
1.614 
0.214 
1.292 
1.625 
0.530 
1.777 
1.452 
1.688 
0.0002 
0.330 
0.778 
1.017 


5 OR v 


1893 cm™ 
942 
332 

1191 
483 
326 
883 
242 
410 
293 

3000 

1050 
730 
610 


473.1°K 





0.213 
1.052 
1.826 
0.740 
1.668 
1.832 
L132 
1.901 
1.751 
1.860 
0.0188 
0.910 
1.344 
1.507 








6 (a) D. P. Stevenson and J. Y. Beach, J. Chem. Phys. 6, 
25 (1938); (b) Geo. Glockler and Walter F. Edgell, J. 
Chem. Phys. 9, 224 (1941); (c) Walter F. Edgell and Geo. 
Glockler, J. Chem. Phys. 9, 484 (1941); (d) C®, for CH2Br2 
was taken from reference 6(c) rather than 6(a) for reasons 
cited in 6(c). 

7R. D. Vold, J. Am. Chem. Soc. 57, 1192 (1935). 

8 Fugassi and Rudy, Ind. Eng. Chem. 30, 1029 (1938). 





HEAT CAPACITY OF CERTAIN HALOMETHANES 


TABLE III. Comparison between semi-empirical and statistical heat capacities.* 














373.1°K(CAL./°K/MOLE) 
C®, CALc. C% Stat. % DEv. 
11.1 ei.3 
14.1 14.1 
17.4 
10.1 
11.9 
14.6 
13.1 
15.6 
16.3 
20.4 
11.5 
14.7 
18.5 
17.1 
17.7 
16.7 
18.1 
20.0 


473.1°K(cAL./°K/MOLE) 
C%® CaLtc. C%» Stat. % DEV, 


12.8 —2.3 
15.9 
19.0 
11.8 
13.8 
16.7 
14.9 
17.5 
18.1 
22.0 
13.2 
16.3 
19.9 
18.7 
19.3 
18.4 
19.5 
21.6 


298.1°K(cAL./°K/MOLE) 


SUBSTANCES C% Catc. C®% Stat. % DEv. 
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* The statistically calculated heat capacities are from the following sources: references 3, 6, 8. 


CHBrF; at 373.1°K to be 


C°, =7.948+ (0.0026+ 1.284+2 K0.595) 
+5/6(0.825+1.503+2X0.438+2 1.622), 
C°,=15.8 cal./°K/mole. 


In Table III a comparison is made between 
values calculated by Eq. (2) in which constant 
stretching and empirical deformation frequencies 
have been used and those calculated by the more 
rigorous statistical method involving experi- 
mentally determined frequencies. A mark (—) 
on the percent deviation column indicates that 
this value was used in evaluating the deformation 


frequencies in Table II. The largest deviation 
(excepting CH;F) is 2.7 percent, while the 
average deviation is +0.9 percent. It has been 
shown’ that heat capacities calculated statisti- 
cally and corrected to finite pressures by the 
usual methods are accurate to +5 percent. 
Table III indicates that the use of Eq. (2) and 
Table II will yield results accurate to 8 percent 
in the temperature range T=250 to 550°K. 
This is more reliable than most of the meager 
data available at present. 

It should be emphasized that Eq. (2) does not 
apply to molecules of the type CX4. 
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The assumption of an immobile hydrous multilayer at a wall of fixed potential was used in 


deriving equations relating zeta potentials, concentrations of univalent salt solution, wall 
potential, and thickness of the immobile layer. Calculations were made using data taken from 
published papers of various experimenters. Values of 8A tc 63A were obtained for the thickness 
of the immobile layer. The values of wall potential and film thickness were considered suf- 


ficiently consistent for each set of data to justify the original assumption. 





HE theory of Freundlich! on the variation of 

electrokinetic (¢) potentials with concen- 
tration of electrolyte involved the hypothesis of 
an immobile layer of electrolytic solution at the 
solid-liquid interface. More recently, Miiller? 
has considered the Freundlich rigidity layer as 
‘probably monomolecular” and treated the inter- 
face involved in ¢ potential phenomena as only 
an adsorbed monolayer containing ions.’ It is the 
purpose of this paper to show evidence in favor 
of a rigid multilayer. 

The ¢ potential may be defined as the electrical 
potential at an interface which is responsible for 
the experimental phenomena of streaming poten- 
tial, cataphoresis, and electro-osmosis. Zeta po- 
tentials are obtained from the experimental data 
by means of a formula by Helmholtz and von 
Smoluchowski.! The formula for streaming po- 
tential (E) is E=¢DPR/4rn where D is the di- 
electric constant, P the pressure forcing liquid 
through a capillary, R the specific resistance of 
the liquid, and 7 the viscosity. 

Abramson and Miiller* derived: an equation 
o=2a(c)'sinh ¢/8 for the charge, o, at the 
interface required to produce the ¢ potential in a 
concentration, c. (a and 6 are constants.) A simi- 
lar derivation utilizing different boundary con- 
ditions results in equations for charge density 
and interface potentials based upon the postulate 
of a rigid multilayer. 

Let ¥ be the potential at a distance X from a 


1H. Freundlich, Kapillarchemie (Akademische Verlags- 
gesellshaft m.b.H., Leipzig, 1922), p. 342. 

2H. Miiller, Cold Spring Harbor Symposium Quant. Biol. 
= Vol. 1, p. 7. 

A. Abramson and H. Miller, Cold Spring Harbor 

Pie Quant. Biol. (1933), Vol. 1, p. 29. 

4 Graetz, Handbuch der Elektricitat und des Magnetismus 
(Barth, Leipzig, 1921), Vol. II, p. 366. 


530 





flat surface in contact with a dilute solution of 
simple salt (see Fig. 1). Yo is the potential of the 
wall with its adsorbed ions. wy is the Nernst 
potential of the wall. Yo will be influenced by 
changes in yy and by adsorption of non-potential 
determining ions having high specific adsorption 
potentials. It will be assumed that monatomic 
ions of low polarizability have a_ negligible 
specific adsorption potential and therefore yy 
and yo are not changed by the addition of small 
amounts of alkali halides. Outside the yo plane 
only Coulomb forces will be considered. 

The potential, y, will decrease with distance, 
X, from the wall as a result of a diffuse double 
layer or Gouy “ion atmosphere”’ in a manner as 
shown in Fig. 1. Let ¢ be the potential at a 
distance, ¢, from the wall, ¢ being the thickness 
of the solution which is assumed not to take part 
in viscous flow during streaming potential deter- 
minations. Then, according to the Boltzmann 
distribution equation, the excess charge density 
II of ions in solution near the charged wall will be: 


Il = —nzeLexp (sey /KT)—exp (—ze¥/KT)], (1) 


where »=number of ions per unit volume of 
solution 
z=valence of the ions (only symmetrical 
valence considered) 
€=unit charge 
K=Boltzmann’s constant 
= absolute temperature. 


More compactly, 
Il = —2nze sinh exp /KT. (2) 


II can be eliminated from the equation by the 
use of Poisson’s differential equation (for a flat 
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surface), 
O*y dr 
—= —— J]. (3) 
0x? D 
Then 
Op S8anze | 
—— sin 


ax? D 


_< (4) 
KT 


ezy Oy 


a ——., (5) 


Ory OY = 16rnze | 
KT 0X 


— sin 


ax2aX  D 
Integrating, 
OW \? 162nKT esy 
(—) = — cosh ——+/, 
ax D KT 
where J is the constant of integration, or 


oy ( 1l6qrnKT exy 


‘OX 


b 
—_— cosh <+1) ; (7) 
D KT 


If dy /AX is negative, approaching zero at infinite 
distance where y is zero, 


oy 16anKT ezy ; 
—=-|— (cosh ~-1)| . (8) 
ox D KT 
Letting (8rnz*e?/DKT)!=Debye’s « and 
a=e2/KT 
aoy 


— = —v2«dX. (9) 
(cosh ay—1)? 





This equation can be rearranged into a standard 
form for integrating. 


a sinh aydy 
(cosh ay — 1)!(cosh? ay—1)! 


= v2«dX. 





(10) 





a sinh ayoy a 
i == fv2ax, (11) 
(cosh ay —1)(cosh ay+1)! 


(cosh ay+1)!—v2 1 
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When X =0, Y=yo and 
(cosh avo+1)!—v2 
Pe (cosh ae-+1)!+V2 
cosh apo /2—1 cosh ay/2—1 
—In " 
cosh a~o/2+1 cosh ay /2+1 
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or 


(13) 
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When X =1t, Y=¢ and 


cosh apo/2—1 
2«t=In — — 
cosh apo 2 of 1 


— (14) 
cosh af/2+1 


o, the charge density on the wall is the negative 
of the integral of the excess charge density from 
the wall to infinity. Or, from Gauss’ law, 
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and Eq. (8) (15) 


m-1)] (16) 
KT 
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o =—— sinh —, 
K 2 


(17) 


o = 35,300c} sinh 19.5yo (18) 


for water solutions if 7=298.2 and D=78.8. 
Formula (18) is strictly true only if the charge 
giving rise to the potential Wo lies in the yp» plane. 
Now if Wo is not a function of concentration, it 
may be eliminated from Eq. (14) in terms of ¢ 
for two concentrations. 


2t(k2—k) 
(cosh a¢,/2—1)(cosh af2/2+1) 


n— -———— - 
(cosh af;/2+1)(cosh afo/2 
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TABLE I. Calculations from published zeta potential data 








ELEc- 


EXPERIMENTER TROLYTE METHOD SOLID 


MOLARITY VoLts VoLts tIN 
c X10? t yo cm X106 





Schénfeldt! KCl 


Electro-osmosis 
oe 


Velisek and 
Vasicek? 


Baborovsky and Streaming 
Burgl? , Potential 


Monagham and 
White! 


Pyrex 

(capil- 

laries) 

DuBois and Glass 
Roberts® shi a 


Streaming 
Potential 
“ec 


Ceramics 
oe 


Electro-osmosis Pyrex 
i. (spheres) 
“cc 


0.03 0.0327 0.0391 
05 0309 .0390 0.29 
07 .0298 .0390 

10 015 053 

1 035 054 115 
316 042 054 

31.6 010 042 
10 .024 056 .075 
3.16 036 .060 
-040 053 
10 .006 052 

1 .029 061 
oh 041 051 

1 0057 028 
75 .0066 .026 
o .0104 032 
29 0113 025 
a 0156 026 
05 .0189 027 
O1 .0228 027 
005 .0262 029 
1 040 135 
a 066 097 
01 095 111 

1 053 107 
Al 083 108 
01 .098 107 
mn .049 .130 
05 052 .104 
01 .079 113 
| 040 048 
OS 044 050 
01 ‘ .048 








1N. Schénfeldt, Zeits. f. Elektrochemie 37, 734 (1931). 


2 J. Velisek and A. Vasicek, Collection of Czechoslovak Chem. Commun. 4, 428 (1932). 
3 J. Baborovsky and B. Burgl, Collection of Czechoslovak Chem. Commun. 3, 563 (1931). 


4 Monagham and White, J. Phys. Chem. 39, 935 (1935). 
5 R. DuBois and A. H. Roberts, J. Phys. Chem. 40, 543 (1936). 


Using Eq. (19) and published data of various 
experimenters, the value of ¢ was calculated. The 
average value of ¢ was substituted back in Eq. 
(14) to determine Yo. The average value of Wo 
was then used in Eq. (18) to calculate o. The 
results are shown in Table I. 

As is evident from the table, it is possible by 
proper choice of a value for ¢ to fit these equations 
to experimental zeta-concentration data. The 
failure of experimental data to fit exactly the 
derived equation may be due in part to errors in 
the data itself and partly to the assumption of a 
constant Wo in deriving the equations used. 

It is of interest to note that in data from 
Monagham and White’s work, yo is the same for 
both methods of measurement while both ¢ and ¢ 
are different. It may be that the flow of liquid 
had a more destructive effect upon the water 
layer in the capillary than in the diaphragm of 


Pyrex spheres. Also, the data from Velisek and 
Vasicek for a ceramics diaphragm show nearly 
the same value (0.05 volt) for sodium chloride, 
potassium chloride and potassium iodide, while 
again ¢ and ¢ differ for the three salts. 

At least to some degree, the assumption of a 
rigid multilayer which incloses part of the 
counter ions of the ion atmosphere at an interface 
makes it possible to account for observed varia- 
tions of ¢ with concentration of simple ions. 
Evidence which would seemingly contradict such 
an assumption may be found in a paper by 
Abramson and Miiller.* They applied the formula 
(20), o=2a(c)! sinh ¢/8, which corresponds to 
Eq. (17), to determine the charge at the interface 
of solid and solution necessary to produce the ¢ 
potential. For a large variety of electrolytes, 
their curves relating surface charge (c) to con- 
centration (c) for concentrations up to 0.01 
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molar closely resemble a Langmuir adsorption 
isotherm for adsorption of gas molecules on an 
inert surface. The adsorbed charge on the nega- 
tive glass surface proves to be negative, whereas 
counter ions in a layer near a negative surface 
would have a predominantly positive charge. 

However, several observations throw some 
doubt on the validity of this evidence. As 
Abramson and Miiller pointed out, at concentra- 
tions above 0.01 molar ‘‘¢ might conceivably 
decrease.”” According to the formula for o, the 
charge would be zero for zero ¢ at any finite 
concentration. Such a decrease in o at higher 
concentrations is foreign to a Langmuir iso- 
therm. Moreover, in the specific case of potassium 
chloride adsorption on glass, Abramson and 
Miiller’s o—c curve shows an increasing negative 
charge on the negative glass surface whereas a 
large difference in specific adsorption potentials 
of two ions so nearly alike in ionic weight and 
volume is not to be expected. Finally, a careful 
analysis of the equation for o will reveal that o 
may consist very largely of charges other than 
ions of the added electrolyte. 

The equation ¢=2a(c)! sinh ¢/8 involves the 
assumption that the entire charge on the wall 
is in the single ¢ potential plane. Verwey® con- 
siders the charges on the wall to be scattered in a 
finite layer, consisting of potential-determining 
ions in the wall, adsorbed ions on the wall and, 
perhaps, counter ions inclosed in an immobile 
water layer along the wall. If the only force 
holding ions in or near a glass surface were elec- 
trostatic (obviously not strictly true), if the 
dielectric constant of the solution were un- 
changed near the surface, and if all the potential 
determining ions of a glass electrode were in one 
plane, an equation similar to Eq. (17) could be 
derived to give ov =4nze/x sinh ayy/2. Using 
Abramson and Miiller’s constants, this would 
become oy = 2ac} sinh py /B. In this case yy would 
be the Nernst potential of the glass as an elec- 
trode, oy would be the charge in the potential 
plane, and c would be as before the concentration 
of ions of all kinds in the bulk solution. The 
Nernst potential of glass in pure water would be 
expected to be the same sign and perhaps larger 
in magnitude than the ¢ potential of glass in 


5E. J. W. Verwey, Chem. Rev. 16, 363-411 (1935). 
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pure water. It follows from the above equation 
that an increase in concentration of the solution 
must result in an increase in oy if py is to remain 
constant. Unless the electrolyte added in ¢ 
potential measurements contains potential-de- 
termining ions, Yy should remain nearly constant. 
According to the equation, the o—c curve for a 
constant yy would be a parabola of increasing 
negative charge resulting from a transfer of 
potential-determining ions. The reason for an 
increase in oy in this case is clear if the electrode 
and its double layer are considered to be two 
plates of a charged condenser. If the plates are 
moved closer together, a greater charge on the 
plates is necessary to maintain the original 
potential. An increase in the concentration of a 
solution results in a greater x, the reciprocal of 
which is Debye’s measure of the effective thick- 
ness of a double layer. That is, an increase in ¢ 
results in a decrease in the effective thickness of 
the double layer. 

Even though the assumptions necessary to 
derive an equation for cy are not strictly valid, 
it is evident that the o in the ¢ potential plane 
which could represent all the charges in the 
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Fic. 2. Variation with concentration of effective 
charge density. 


finite wall surface is in a large measure repre- 
sentative of the adsorbed potential-determining 
ions rather than ions of the added electrolyte. 
The first approximation of the charge of ad- 
sorbed potassium and chloride ions might be 
obtained by subtracting Abramson and Miiller’s 
curve (Eq. (20)) from a true parabola (Eq. (17)), 
(see Fig. 2). The resulting o for adsorbed counter 
ions would be positive instead of negative and 
would comprise the part of the positive ion 
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atmosphere in the immobilized layer between the 
yo» and ¢ planes. This ¢ might not necessarily 
follow a Langmuir adsorption isotherm. 


SUMMARY 


(1) Equations have been derived for the rela- 
tions between wall potential, zeta potential, and 
concentration of univalent salt solutions, assum- 
ing an immobile hydrous layer of finite thickness 
was attached to a wall of fixed potential. 

(2) When published ¢ potential-concentration 
data for simple salt solutions are substituted in 
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an equation relating film thickness to ¢ potential 
for a constant wall potential, values of 8A to 63A 
are obtained for the rigid layer thickness. 

(3) The results show that ¢ potential-concen- 
tration curves for simple univalent ions can be 
explained without attributing such large specific 
adsorption potentials to the ions as the o—c 
curves of Abramson and Miiller might indicate. 

(4) The hypothesis of a rigid or viscous water 
multilayer at a solid-water interface, though not 
proved, is at least tenable in the light of its 
ability to fit certain ¢-potential data. 
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The spreading phenomena on water of drops of mineral 
oils containing a wide variety of organic acids and amines 
having but one polar group have been studied under con- 
ditions where dissociation at the oil-water interface could 
most readily occur. Particular emphasis is placed on the 
homologous series as an aid in separating the variables 
involved. Various phenomena of edge diffusion are 
described and a qualitative theory is presented. Edge 
diffusion is shown not to be a serious source of error 
in measuring the cross-sectional areas of the adsorbed 
molecules if there are over 13 carbon atoms per straight 
chain. The different types of film spreading and break-up 
effects are described and the various physical and chemical 
factors involved are discussed. The changes in the spread- 
ing effects due to varying the pH and to the presence in 
the water of Ca++, Cut+, Pb**t, Fet**, Al*+**, Latt+, and 


1. INTRODUCTION 


T has long been known that highly purified 
mineral oils such as medicinal petrolatum do 

not spread on the clean surface of water but form 
lenses having well defined contact angles with the 
water. Such lenses were studied by Hardy,! 
Coghill and Anderson,” Lyons,’ and more recently 


1W. B. Hardy, Proc. Roy. Soc. A86, 610 (1912). 

2 W. H. Coghill and C. O. Anderson, Tech. Paper 262, 
U. S. Bureau of Mines, 1924. 

8C. G. Lyons, Trans. J. Chem. Soc. (London) 623 
(1930). 


Th**** on the interfacial film are described. It is shown 
that only the long chain monobasic saturated acids form 
rigid films when affected by metallic ions, while other 
acids are either condensed somewhat or are unaffected. 
Short chain metallic salts of these fatty acids are often 
oil soluble and have brief lifetimes of adsorption. An 
experimental proof is given that dissociated long chain 
acids and amines have infinite lifetimes of adsorption. 
Evidence is given for concluding that the lifetime is short 
for molecules having less than 14 carbon atoms. The effect 
of Cat* on the spreading of oleic-stearic acid mixtures is 
given in some detail and also a theory of the results. 
Reasons are given for concluding that the adsorbed acid 
molecules form two-dimensional film complexes with 
hydrocarbon molecules derived from the oil. 


by Langmuir.* Harkins and Feldman® have 
shown that a drop of a pure liquid would not 
spread on water if its spreading coefficient Fs 
were negative; and in agreement with this, 
values of Fs for mineral oils, calculated from 
surface tension data, range usually between —11 
and —14 dynes per centimeter. 

When such oils are either heated or radiated 
with ultraviolet light in the presence of air they 

4]. Langmuir, J. Chem. Phys. 1, 756 (1933). 


5 W. D. Harkins and A. Feldman, J. Am. Chem. Soc. 44, 
2665 (1922). 
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develop the ability to spread on water and the 
colored films due to the spreading of oxidized 
engine oils are a familiar sight on wet city streets 
and around waterfronts. Langmuir® pointed out 
years ago that such spreading was probably due 
to the presence of polar molecules developed in 
the oil as oxidation products, and that acid 
molecules must be particularly effective. Later 
optical properties of thermally oxidized lubri- 
cating oils were described by Blodgett’ and their 
value for use as “‘indicator”’ or ‘‘piston’’ oils was 
subsequently shown by Langmuir and Schaefer.* 

An extensive literature exists’ covering various 
physical and chemical properties of monolayers 
adsorbed at the air-water interface, but compara- 
tively little is known about the monolayers ad- 
sorbed at the oil-water interface. Yet adsorption 
phenomena at the oil-water interface are techno- 
logically and biologically more common. 

In recent years there has been an increasing 
amount of research on the interfacial adsorption 
of polar molecules initially in solution in one of 
the two liquid phases, employing interfacial 
tension techniques, but due in part to the experi- 
mental difficulties encountered in separating the 
numerous variables involved, much remains to 
be learned before there will be a satisfactory 
understanding of the fundamental phenomena 
involved. Particularly pertinent to the subject 
matter presented here is the well-known study of 
the adsorption of short chain acids by Harkins 
and his pupils,'®' the work on the effect of pH 
and buffers by Peters! and by Danielli," and of 
Alexander and Teorell' as well as of Alexander 
and Schulman!® on the equation of state of 
interfacial films. 

In the past decade in three stimulating papers, 


*T. Langmuir, Trans. Faraday Soc. 15, 62 (1920). 

7K. B. Blodgett, J. Opt. Soc. Am. 24, 313 (1934). 

‘I. Langmuir and V. Schaefer, J. Am. hon Soc. 59, 
2400 (1937). 

*N.K. Adam, Physics and Chemistry of Surfaces (Oxford 
University Press, 1938), second edition. 
- W. D. Harkins and H. King, J. Am. Chem. Soc. 41, 970 
(1919). 

“ W. D. Harkins and H. McLaughlin, J. Am. Chem. Soc. 
47, 1610 (1925). 

” R. A. Peters, Proc. Roy. Soc. A133, 140 (1931). 

3 J. F. Danielli, Proc. Roy. Soc. B122, 155 (1937). 

A. Alexander and T. Teorell, Trans. Faraday Soc. 35, 
727 (1939). 

46 4. Alexander and Schulman, Trans. Faraday Soc. 36, 
960 (1940). 
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Langmuir and Blodgett*!®'’ have presented their 
results on the properties of interfacial films of 
stearic acid obtained by quantitative studies of 
the spreading on water of drops of a solution of 
that acid in mineral oil. They also described the 
remarkable sensitivity of the oil drop spreading 
to the presence of polyvalent metallic ions as well 
as to the effect of the pH of the aqueous sub- 
strate. They concluded that every stearic acid 
molecule was eventually adsorbed at the oil- 
water interface at alkaline pH, while at acid pH 
the solubility in oil of the stearic acid prevented 
spreading until a critical concentration was 
exceeded. 

Ford, McBain and Wilson'® described the 
spreading of oleic acid in mineral oil and reported 
that the edge of the disk was preceded by a 
monolayer of the acid. This was ingeniously done 
by the use of rapid motion pictures of the 
spreading drop of oil. 

Stenstrom and Vigness'® described experiments 
on the spreading of oil drops containing stearic, 
oleic and linoleic acids on alkaline water and 
found that oil drops containing either oleic or 
linoleic acids spread on alkaline water, but unlike 
stearic acid would also spread on alkaline water 
containing calcium ions. 

Allen, Grant and Burdon”® have investigated, 
using a motion picture camera, the rate of 
spreading on water of oil drops which had pre- 
viously been irradiated with ultraviolet light. 
They found that in the latter stages of spreading 
the diameter of such an oil disk was accurately 
proportional to the square of the time. 

Of particular interest here are Mercer's”! ex- 
periments on the spreading on alkaline water of 
oil drops containing stearic acid. It was demon- 
strated that during the spreading of an oil drop 
there was shed or given off from the boundary 
of the disk a monolayer of stearic acid. Raising 
the temperature was found to cause this effect 
to become more prominent. 


16 J, Langmuir, J. Frank. Inst. 218, 143 (1934). 

177, Langmuir, Science 84, 382 (1936). 

18 Ford, McBain and Wilson, Kolloid Zeits. 78, 1 (1937). 

19 W. Stenstrom and I. Vigness, J. Phys. Chem. 43, 531 
(1939). 

20 Allen, Grant and Burdon, Trans. Faraday Soc. 33, 
1531 (1937). . 

#1 E. H. Mercer, Proc. Phys. Soc. London 51, 561 (1939). 
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Recently Askew and Danielli®? have studied 
the aqueous state of the adsorbed interfacial 
monolayers using a barrier arranged at the 
liquid-liquid interface, but due to experimental 
difficulties this more fundamental approach has 
been limited to a few materials. However, they 
came to one conclusion of interest here, that is, 
trimethyl cellulose and alpha-amino palmitic 
acid had adsorption areas per molecule consider- 
ably higher at the brombenzene-water interface 
than at the air-water interface. Murray® has 
investigated the effect of varying the polar 
content and nature of the oil used in the spread- 
ing of oils with emphasis on their use for mosquito 
control. Measurements were made of their 
spreading pressure and durability after spreading. 


2. OUTLINE OF WORK TO BE PRESENTED 


The plan of attack here adopted for the study 
of the spreading of oil drops on water was first 
to work entirely with one pure mineral oil while 
varying the nature of the polar additive com- 
pound dissolved in it (with emphasis on the 
employment of the homologous organic series as 
a tool to aid in unraveling the various chemical 
and physical factors involved in the mechanism 
of spreading), and second, to study each of the 
various factors so found as they are affected by 
altering the nature of either the mineral oil or 
the pure hydrocarbon compound employed. 

Accordingly, unless otherwise stated, the min- 
eral oil used will be understood to be a USPX1 
light petrolatum manufactured by L. Sonneborn 
Sons Company, Inc., of New York City, and 
having a density varying from 0.839 to 0.841 and 
an absolute viscosity of 19.72 centipoises at 25°C. 
This oil does not spread on water. Employing 
both a Cenco duNouy interfacial tensiometer and 
a drop weight apparatus and applying the 
Harkins and Jordan corrections,” surface tension 
measurements on this oil gave for the interfacial 
tension oil-water at 25° the value of 53.5 dynes/ 
cm, and for the surface tension oil-air 29.8 
dynes/cm. Using 71.97 dynes/cm for the surface 


22 F, Askew and J. Danielli, Trans. Faraday Soc. 36, 785 
(1940). 

231. Murray, Bull. Entomological Research 30, 173 
(1939). 

24 W. Harkins and H. Jordan, J. Am. Chem. Soc. 52, 1751 
(1930). 
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tension air-water, one obtains for the Harkins 
spreading coefficient, F,, the value of —11.3 
dynes/cm. 

If one regards the spreading of oil drops as a 
means of studying interfacial films, it is evident 
that it is restricted to studying those polar com- 
pounds capable of being dissolved in the mineral 
oil used. As mineral oils oxidize rapidly at 
temperatures around 100°C, and as it is so often 
true that a polar compound dissolves little in 
oil unless heated close to its melting point, it 
has not been found possible here to study addi- 
tive compounds whose melting points exceed 90° 
to 100°C. It is also a restricted technique in the 
sense that with it the films can only make their 
presence noted when the surface film pressure is 
equal to or greater than the compression due to 
surface tension forces inherent in the existence 
of a floating drop, i.e., for the above-mentioned 
mineral oil it is 11.3 dynes/cm. 


3. TECHNIQUE 


The oil drops were spread on water held in 
Pyrex trays while Pyrex barriers were used for 
scraping clean the surface of the water. All 
aqueous solutions were handled entirely in Pyrex 
containers which had previously been well cleaned 
with chromic acid and aged in hot distilled water. 
The distilled water was collected from a metal 
still which was designed to give unusual protec- 
tion from spattering and entrainment, and had a 
solid tin condenser thus largely eliminating con- 
tamination from traces of copper so serious in 
film research. For some experiments an all- 
quartz still was used to redistill the water. 

The amount of oil in each disk was measured 
by weighing it in Pyrex droppers using a chaino- 
matic balance good to 0.1 milligram; hence the 
error of measuring the quantity of oil used in 
forming a spreading film was never greater than 
0.5 percent for the smallest disk studied and pro- 
portionately less for the larger disks. 

All measurements, unless otherwise noted, 
were made with the water at 25°C+0.2° while 
the room air temperature was usually maintained 
within 0.5° of that. 

The hydrophil tray was arranged inside of a 
metal box having a glass top in order to prevent 
air currents from distorting the circular spreading 
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oil disks. Two thin metric rulers arranged at 
right angles close to but above the plane of the 
water in the tray permitted measurement of two 
mutually perpendicular diameters of the oil disk 
whenever air or water currents caused a slight 
distortion from a circular shape. 

The obtaining of the wide variety of polar 
compounds necessary for these experiments was 
a most serious problem, particularly since in the 
course of the early experiments it was found that 
certain of the phenomena studied were sensitive 
to the presence of organic polar impurities. 

The pentadecanoic, nonadecanoic and eico- 
sanoic, also a remarkably complete series of 
branch chain acetic acid derivatives having from 
14 to 20 carbon atoms per molecule, and the 
chaulmoogric acid and the 13-hydroxy stearic 
acid, many of them being of unusual purity, 
were furnished this laboratory by Dr. Roger 
Adams of the Department of Chemistry of the 
University of Illinois. The oleic and elaidic acids 
were of exceptional purity (for example, the 
oleic acid was colorless and odorless and had an 
iodine number of 88.3) were newly prepared by 
Dr. J. P. Kass of the Department of Botany of 
the University of Minnesota. The ricinelaidic 
acid was a purified specimen prepared by Pro- 
fessor St. Elmo Brady of the Department of 
Chemistry of Fiske University and has been 
described elsewhere.*®> The phenylstearic, ethyl- 
phenylstearic, xenylstearic, and tetrahydronaph- 
thylstearic acids were obtained from Messrs. 
Stirton, Peterson and Groggins of the Bureau of 
Agricultural Chemistry and Engineering, De- 
partment of Agriculture, and have been described 
elsewhere.2* The normal octadecyl, hexadecyl, 
tetradecyl, dodecyl, and the mixed octadeceny] 
and octadecendienyl amines were obtained from 
the Chemistry Division of Armour & Company, 
while the heptadecyl amine was prepared in this 
laboratory using pure heptadecyl amine hydro- 
chloride made with great care by the Varlacoid 
Chemical Company of New York City. The 
cooperation and generosity of these individuals 
and organizations are gratefully acknowledged. 

The ricinoleic acid was obtained freshly pre- 
pared from the City Chemical Company of New 


*° St. Elmo Brady, J. Am. Chem. Soc. 61, 3464 (1939). 
6 Stirton, Peterson and Groggins, Ind. Eng. Chem. 32, 
1136 (1940). 
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York City, the trioctyl amine was obtained from 
the Sharples Solvents Corporation, and the 
dioctyl amine was obtained from the Carbide and 
Carbon Chemicals Corporation. The remaining 
chemicals used for the work presented here were 
newly obtained Eastman Chemical Company 
preparations. 


4. EpGeE DIFFUSION 


The escape of stearic acid molecules from a 
spreading oil drop previously described'*”! was 
found to be a common phenomenon. Mercer’s 
technique of using aluminum powder, though 
suitable for photographic purposes, was discarded 
in order to avoid the consequent metallic con- 
tamination of the water and use was made of 
either grease-free talc or clean, precipitated and 
washed sulfur powder which was spread on the 
water from a previously flamed platinum mesh 
bag. The escaping molecules push the talc par- 
ticles away from the disk forming thereby a 
circular ring which encloses the escaped mono- 
layer. The greater the rate of escape of the polar 
molecules, the more rapid the rate of increase of 
diameter of the ring of talc. 

If the escaped monolayer is insoluble in water, 
the talc ring does not contract when the oil disk 
ceases emitting polar molecules. If there is no 
reaction between the adsorbed molecules and 
polyvalent metallic impurities, the edge loss 
continues as long as there are adsorbed molecules 
present. If the escaped monolayer is soluble in 
water, the ring rarely spreads more than a few 
centimeters from the disk and remains so. There 
is then an equilibrium between the rate of escape 
of the molecules from the oil disk and their rate 
of subsequent solution in the water. This can 
readily be shown by dropping a few talc particles 
in the space between the ring and the oil disk 
and observing the still rapid motion of the talc 
away from the edge of the disk even though the 
circular talc boundary still maintains itself at a 
fixed distance from the disk. 

The existence of rapid edge diffusion can often 
be detected by observing the difficulty of making 
two oil disks join each other. If there is edge 
diffusion they will appear to repel one another, 
sometimes with considerable speed. 

Experiments with the normal saturated acids 
from hexanoic acid (Cs) to eicosanoic acid (C29) 
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using equimolal concentrations in the petrolatum 
demonstrated that each showed edge diffusion, 
but whereas in the case of the Coo acid the effect 
was barely detectable, it was so pronounced for 
the Cs, C; and Cs acids as to be the most evident 
spreading effect caused by their presence in oil. 
In fact, if a weight concentration of from 10~‘ to 
10-* was used, the oil drop did not spread on the 


water but skated about while the nearby talc ° 


particles were propelled violently away from the 
disk. Now, when a drop of pure hexanoic or 
octanoic acid was held in the air just above the 
surface of the talc-covered water, the acid vapor 
from the drop also caused spreading of the talc. 
It therefore might be concluded as in the classic 
study of the propulsion of camphor particles on 
water” that the edge diffusion effect of such acids 
dissolved in oil drops is entirely due to the vapor 
spreading from the oil disk through the air 
followed by condensation on the water. However, 
upon employing decanoic acid (Cio) or higher 
homologs, no vapor spreading could be ob- 
served by the simple tests described above, and, 
consequently, it was concluded that this mecha- 
nism, though contributory in the case of the 
short chain length acids, was of minor conse- 
quence with the longer chain length molecules. 
It was found that equimolal solutions in oil 
of alcohols, esters, ketones, amines, phenols, phos- 
phates, and ethers all showed edge diffusion 
effects which decreased in intensity as the chain 
length of the hydrocarbon portion of the mole- 
cule increased. In fact, the effect for all the 
straight chain molecules containing the same 
number of carbon atoms per molecule seems 
roughly to be of the same intensity whatever the 
polar group. In a polar molecule containing 
several branch chains the length of the longest 
branch and not the total number of carbon atoms 
in the molecule seemed to determine the intensity 
of edge diffusion. Thus, either dioctyl or trioctyl 
amine caused a much more intense edge loss than 
did hexadecyl amine and either dibenzylacetic 
or di-n-heptylacetic acid much more than 
hexadecanoic acid. Again using a series of 
branched chain acetic acid derivatives, each 
having 17 carbon atoms per molecule ranging 
from methyltetradecylacetic acid to heptyloctyl- 


3 Marcelin, Solutions Superficielles, (1931), Chap- 
ter A. 
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acetic acid, it was found that the rate of edge 
diffusion increased rapidly as the length of the 
largest branched chain decreased. Many more 
such illustrations could be given. 

Although the greatest intensity of edge loss is 
shown by short chain molecules, and although 
it is true that the shorter the chain length the 
lower the melting point and the higher the vapor 
pressure of a molecule, it is not true that a 
simple melting point rule can be used for pre- 
dicting relative rates of edge loss, for the mo- 
lecular configuration is too important. Thus oleic 
acid with an 18 carbon chain has a melting 
point below 20°C while dibenzylacetic acid with 
16 carbon atoms per molecule has a melting 
point around 81°C. Yet the former exhibited a 
barely noticeable edge loss while the latter 
showed a most pronounced effect. Here the low 
melting point of oleic acid is due to its cis- 
configuration, which, however, does not cause 
increased edge diffusion, while the intense edge 
loss of dibenzylacetic acid is due to its being a 
branch chain aromatic molecule having only 7 
carbon atoms in its longest chain. However, it is 
true that the rate of edge loss increases as the 
melting point decreases for any homologous 
series of molecules containing only one straight 
chain per molecule. 

It is concluded that edge diffusion is essentially 
a kinetic phenomenon, a 2-dimensional evapora- 
tion occurring across a lineal boundary and 
requiring the presence of two liquids, one more 
polar than the other. Its mechanism in the case 
of long chain polar molecules presumably in- 
volves the adsorption of the polar molecule at 
the oil-water interface, followed by thermal 
agitation while adsorbed there and a certain 
amount of associated diffusion in the interface, 
and finally in the case of the adsorbed molecules 
nearest the rim of the oil disk, a more frequent 
breaking away of the hydrocarbon ends of the 
molecules from the oil during which time the 
polar end of the molecule continues to adhere to 
the water phase. The longer the chain length the 
slower the rate of two-dimensional diffusion in 
the interface and the less the frequency of 
breaking away of the hydrocarbon chains of the 
polar molecules nearest the rim or edge of the 
disk, since the rupture of a number of hydrogen 
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bonds is necessary for such an event to occur. 
It is possible that this effect varies with the 
nature of the hydrocarbon in the oil used. 


5. Some Factors CONTROLLING THE 
STABILITY OF OrL FILMS 


It was found from a comparative study of the 
breaking up of oil films that if the amount of the 
polar molecules dissolved in the oil was sufficient 
just to cause spreading as a thin disk, i.e., a 
weight concentration of 10-5 in case of the long 
chain molecules and 10-* in case of the short 
chain ones, the mode of breaking up of the disk 
was qualitatively indicative of the greatest chain 
length of the polar molecule. Compounds with 
but 6 to 9 carbon atoms per molecule spread in 
a second and at once burst into hundreds of tiny 
drops, each of which exhibited a considerable 
contact angle with the water while skating about 
for a period 20 seconds or less, depending on 
chain length. But if there were approximately 
from 10 to 13 carbon atoms per molecule, the 
disk spread in a few seconds and then broke up 
by a streaming of small flat disks from the 
periphery of the central disk leaving a constantly 
contracting central disk. These disks initially 
spread to a larger diameter, broke off more 
slowly, and after breaking off, contracted more 
slowly, the greater the chain length of the 
molecule. Finally, if there were 14 or more 
carbon atoms per molecule, the disk spread more 
slowly and to a greater diameter (exhibited 
surface oil flow effects described below) and 
broke up more slowly. These break-up effects can 
be of remarkable regularity and are fascinating 
to observe. - 

The mechanism operating to cause such a 
variety of break-up effects, of course, cannot be 
simple. However, the use of a wide range of 
concentrations has made it possible to distinguish 
four factors responsible for these phenomena. 
They are: first, the rapid increase in the speed of 
diffusion through the volume of the oil disk of 
the polar molecules as their chain length de- 
creases ; second, the rapid increase in the rate of 
edge diffusion as the chain length decreases; 
third, the inertia of the oil in the upper layer of 
the oil in a spreading disk; and fourth, the 
viscosity of the oil in causing a dampening effect 
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on any inertia flow or surging in the upper oil 
layer of the disk. 

The first factor is evident, since the rate of 
diffusion of a long chain molecule should increase 
as the chain length decreases. The fewer carbon 
atoms per chain the fewer bonds to be made and 
broken with the hydrocarbon solvent and also 
the less the probability of collision with other 
molecules. To illustrate this, if solutions of the 
same molal concentration are used corresponding 
to w=0.90X10-5 stearic acid (i.e., 3.1810-° 
gram molecule per gram oil), then the time, Ty, 
required for the disk to spread to its maximum 
diameter on alkaline water free from polyvalent 
metallic ions is found to be approximately 3 
minutes for tetradecanoic, 5 minutes for hexade- 
canoic, 63 minutes for octadecanoic, 8 minutes 
for nonadecanoic, and 9} minutes for eicosanoic 
acid. Since the time required for the actual 
adsorption at the interface must be extremely 
small, and since for a given molal concentration 
the total spreading area on such a substrate is 
approximately the same independent of chain 
length (see Section 8 below), this progressive 
increase in Ty can only be due to a corresponding 
decrease in the speed of molecular diffusion with 
increased chain length in a homologous series. 
Again, using the same molal concentrations of 
various 18 carbon aliphatic acids, it was found 
that stearic, oleic, and elaidic acids had the 
same value of 7y=6.5 min., whereas ricinoleic 
and ricinelaidic acids each had the large value of 
Tu=50 min. Hence the double bond in the first 
case was of no influence, but the hydroxy group 
in the second case was of great influence, and, 
in fact, may have caused association. 

The second factor, edge diffusion, has already 
been described. It was always evidenced by the 
behavior of the smaller oil disks formed after the 
break-up, for the disks so formed had a greater 
contact angle, showed more skating, were smaller 
in diameter, and contracted more quickly, the 
greater the rate of edge diffusion. Where thin 
disks exhibited considerable edge diffusion, the 
disks gradually contracted in diameter as a 
beaded rim formed. This was a common phe- 
nomenon. 

The third and fourth factors were made most 
evident by observing the spreading of oil disks 
of different viscosities containing acids of 16 or 
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Fic. 1. Damped oscillation of a spreading oil disk. 


more carbon atoms per molecule and present in 
weight concentrations in the range 2X10-* to 
5X10-*. Using light petrolatum and this range 
of concentrations, the spreading oil drop required 
from 8 minutes to a few seconds to spread to its 
maximum diameter. Under these circumstances 
the inertia of the upper oil layer caused it to 
pile up at the circumference so that at maximum 
spreading the disk exhibited a uniform rim 
bounding its circumference. This rim was nar- 
rower and higher the greater the speed of 
spreading (i.e., the greater the value of w) and 
the less viscous the oil. If w=1 10-5, the rim 
could not be seen, whereas if w was around 
5X10-° it could not be missed. At the low 
range of spreading speeds after the rim formed, 
it slowly subsided so that the final result was a 
critically damped oscillation (see Fig. 1 (a), (b), 
(c) and (d). At the high range of spreading speeds 
the flow of oil from the center of the disk caused 
a hole to form there; the hole grew rapidly until 
only a thin oil ring remained. Finally, there 
developed a necking-in or narrowing at various 
points on the internal circumference, and finally 
at these points necking-in progressed until these 
had caused the formation of a number of oil 
disks, see Fig. 2 (a), (b), (c), (d), (e) and (f). 
If a much more viscous oil like Squibbs heavy 
petrolatum is used with w=5X10-°, spreading 
though slower, resulted in the formation of a 
low and wide oil rim. After the disk had reached 
its maximum diameter, if air currents did not 
distort the disk enough to start a breaking-up 
process, the rim gradually became transformed 
into a border of oil mounds which in time became 
fewer in number and shallower until eventually 
the disk was flat. 
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As an oil drop spreads, the thickness of the 
oil layer will decrease and the polar molecules in 
solution will have less oil through which to 
diffuse in reaching the oil-water interface. There- 
fore, they will arrive at the interface at a rapidly 
accelerated rate until the maximum spreading 
area of the disk is approached. Hence, Ty will 
increase with great rapidity as the initial con- 
centration of the polar molecules is increased. 
It was found in the case of stearic acid in mineral 
oil that Ty was 6.5 minutes when w was 0.85 
10-5, 1.5 minutes when w was 210-5, and 25 
seconds when w was 2.5X10-*. Again in the 
case of a ricinoleic acid, 7, was 40 minutes when 
w was 1X10-° and 7 minutes when w was 
2x10-*. These are in reasonable accord with 
Langmuir’s theoretical conclusion'® that the 
initial rate of spreading of the oil disk should 
increase exponentially with time and with the 
square of the initial concentration. 

At still higher spreading speeds the various 
factors operating are no longer readily dis- 
cernible. If the speed of diffusion is considerable, 
the spreading effect of the polar molecules in 
the oil of the inertia layer as it flows towards 
the circumference is an important part of the 
whole and tends to nullify the trend to form a 
high rim. Hence at concentrations in the range 
w=10-* to 10-* rims are not usually evident. 
Since under such circumstances the break-up is 
by a streaming of small disks away from the 
circumference, it is possible that a border of oil 
mounds formed and in breaking off were the 
origin of these small disks; however, that is a 
conjecture. At the highest speeds, such as are 
obtained with long chain acids when w exceeds 
10-*, there is an additional phenomenon, a 
bursting of the oil film in its center to form a 
large hole, as if due to the tearing of a tight 
membrane. 

In addition to these causes of film instability, 
there are two due to solubility effects. First, if 
the interfacial film is so soluble in one of the 
liquid phases as to be forced to leave the interface, 
the oil film will contract from its maximum area. 
Second, if any polar material is undissolved in the 
oil, the oil disk in spreading will develop numer- 
ous holes due to the increased local spreading 
occurring whenever particles of that polar ma- 
terial reach the interface. This phenomenon can 
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be used to indicate the attainment of the 
solubility limit of polar substances in oils. 
Measurements of this type will be described in 
Part II. 


6. METHOD OF MEASURING AREA PER MOLECULE 


Langmuir and Blodgett'*'’ found the maxi- 
mum area of spreading of an oil disk on alkaline 
water to be the same whatever the volume of the 
oil so long as the same total number of stearic 
acid molecules were present in the disk, and hence 
concluded that every stearic acid molecule had 
been adsorbed at the oil-water interface. Mercer”! 
came to the same conclusion upon finding that 
after spreading had occurred oil obtained by 
sucking up some of the upper oil layer of the 
disk and in turn placing it on alkaline water no 
longer had the ability to spread. Assuming the 
correctness of their conclusion, if the average 
area occupied by each adsorbed molecule at the 
oil-water interface is A, if M is the gram molecu- 
lar weight of the polar substance, present in a 
weight concentration w, and if the mass of the 
oil in the disk is m, then the total area S of the 
interface between oil disk and water is given by 
S=6.06 X 10*2(mw/M)A. For very thin circular 
disks of radius R evidently S=aR*. For thicker 
disks S>R? since it is no longer a good approxi- 
mation to treat the interface as a plane surface. 

The procedure used here was to calculate an 
apparent area per molecule A’ obtained by the 
use of the above formula and the assumption that 
S=7R?. If w was not too low, a plot of A’ asa 
function of m revealed that A’ increased with m 
until a certain value my was reached, after which 
A’ became constant and equal to the true area 
per molecule. For example, in the case of stearic 
acid (w=0.9 X 10-5) in oil spread on 0.02 N NaOH 
+0.02 N Na2COs, it was found that if mo ex- 
ceeded 4.6 X 10? gram, then A’ = 106 X 10-"* cm?; 
if m=2.5X10-? gram, then A’=95X10-"* cm?; 
and if m=0.8X10- gram, then A’=84X10-'* 
cm*. 

In handling the more dilute oil solutions of 
polar molecules adsorption upon the walls of the 
glass containers employed can be a source of 
error. A simple calculation shows that for solu- 
tions having a weight concentration of 10-5, the 
maximum possible error will be 2 percent; for 
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higher solution concentrations the maximum 
possible error will be proportionately less. 

In order to obtain accuracy in measuring A, it 
is essential to choose carefully the value of w 
employed. If w ranges from 210-* to 4X10-, 
the disk spreads in a circle but does it so rapidly 
and breaks up so quickly at the edges after com- 
pleting the spreading that a rapid measurement 
of the diameter is necessary. If w is less than 
110-5, the curvature error involved in as- 
suming that the interfacial area equals 7R? 
becomes serious and also the edge loss error 
increases, since under such circumstances 7’ is 
large. The best results obtained were with either 
of the ranges w=1X10-5 to 510-5 and 210-4 
to 4X10-. The effect of the polyvalent metallic 
ions on the maximum diameter of spreading was 
less the smaller 74 and so unless the water was 
very carefully freed from traces of metallic ions, 
it was found best to use the higher concentration 
range. However, using either range and properly 
purified water, the same values of A were 
obtained. 


7. THE EFFECT oF pH 


The spreading areas of oil drops containing 
long chain acids vary with the pH of the aqueous 
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solution used. In Fig. 3, A gives the change in 
the diameter with pH of oil drops containing 
w=0.890 X 10-* stearic acid, the pH having been 
varied by adding to distilled water in equilibrium 
with the CO, of the air a suitable amount of HCl, 
KHCOs, or NaOH. Practically all of the increase 
in the spreading occurred in the pH range from 
pH 9 to pH 10.5. The value of this upper limit 
agrees with Langmuir and Schaefer’s results” 


28 1. Langmuir and Schaefer, J. Am. Chem. Soc. 58, 284 
(1936). 





























































































































































































































542 W. A. 


on the effect of pH in converting to the calcium or 
barium soap a monolayer of stearic acid at the 
air-water interface. A similar study of the effect 
of pH on the diameter of oil drops containing 
w=0.865 X10-* oleic acid is given in Fig. 3, B, 
where it is seen that the maximum effect is also 
at pH 10.5. This is in agreement with the results 
of Danielli'® on the effect of pH on the interfacial 
tension lowering due to oleic acid adsorbed at the 
brombenzene-water interface. 

The obvious conclusion to be made is that of 
the molecules in adsorption equilibrium at the 
interface, some are ionized by the water and 
remain permanently adsorbed. As the pH rises, 
the higher concentration of hydroxyl ions in- 
creases the number of ionized molecules present, 
and for values of pH exceeding 10.5 all are 
ionized and maximum spreading results. The fact 
that the same pH range is effective in causing 
maximum spreading of either the oleic or the 
stearic acids is evidence for inferring that this 
also is true for all long chain acid molecules. It 
need not be assumed that all of the permanently 
adsorbed acid molecules remain ionized, but 
instead that they have been converted to the 
undissociated sodium soap. However, it is theo- 
retically more reasonable to assume practically 
complete ionization of the sodium soap. The 
electrical attraction which then must exist be- 
tween the two portions of the double layer 
consisting of the negatively charged ionized 
carboxy! groups adhering to the oil disk and the 
corresponding sheath of sodium ions in the 
substrate immediately beneath will give rise to 
forces of attraction responsible for the infinite 
lifetime of adsorption of fatty acids at high pH. 

Similar results on the effect of PH were ob- 
tained with the other long chain saturated 
straight chain acids, the unsaturated acids, oleic, 
elaidic, linoleic, linolenic, and erucic; the aryl 
stearic acids, phenylundecylic, phenylstearic, 
ethylphenylstearic, diethylphenylstearic, xenyl- 
stearic, and tetrahydronaphthylstearic; chaul- 
moogric and hydnocarpic acids; dibenzylacetic 
and di-n-heptylacetic acid; and finally with a 
variety of so-called naphthenic acids. 

Oil drops containing acids in a weight concen- 
tration of less than w= 10- and having less than 
13 carbon atoms per molecule showed edge 
diffusion mainly. As has already been described, 
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at higher concentrations they spread so rapidly as 
to seem to spread and burst into hundreds of tiny 
drops in the space of a few seconds, the speed 
being greater the shorter the chain length. 

In the case of such molecules the effect of 
varying the pH was difficult to study, but roughly 
speaking, it was found that there was less change 
in oil drop spreading with the pH the shorter 
chain length. Only a decrease in the rate of edge 
loss with increased pH was noticed with the 6, 7 
and 8 carbon acids. Evidently the film pressure 
due to these acid molecules is not sufficient to 
overcome the —11.3 dynes/cm compression of 
the oil disk unless w considerably exceeds 10-*. 
Therefore dodecanoic acid is the shortest molecule 
of the series readily to exhibit disk spreading. 
But the shorter the chain length the greater the 
film pressure for constant area of packing per 
molecule. Hence the conclusion seems inevitable 
that for straight chains of less than 12 to 14 
carbon atoms per molecule, the lifetime of 
adsorption is not infinite but is less the shorter 
the molecule. 

Acids containing more than one polar group 
such as hydroxy acids behaved quite differently 
with respect to changes in pH, and for the sake of 
a logical order of presentation will be discussed in 
Part III. 


8. THE Errects DUE TO POLYVALENT 
METALLIC IONS 


Langmuir and Blodgett!’ have described the 
remarkable ability shown by small traces of 
calcium ions in alkaline water to prevent any 
spreading of oil drops containing stearic acid. 
The simple expedient of probing such an oil disk 
with a clean thin glass rod or a platinum wire will 
show that the disk is floating as if on a rigid 
island. A systematic qualitative study was made 
of the effects of ions of Ca, Pb, Cu, Fe, Al, La or 
Th upon the interfacial films of the acids men- 
tioned above. The metallic salts were either the 
chloride or the nitrate and were employed in 
concentrations of 10-* to 10-* molal. The pH was 
used as high as possible without causing com- 
plete precipitation of the metallic ions. The 
procedure used was to observe the rate of 
spreading and subsequent contraction of oil 
drops containing the fatty acid while talc and 
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probe tests were also made for edge loss and film 
rigidity. Drops containing the acid in weight 
concentrations of 10-5, 10-4 and 10- were tested. 
As is to be expected, the amount of the metallic 
salt necessary to cause rigidity was found to 
depend on the pH but it also depended on the 
concentration of the acid in the oil. When there 
was doubt as to the rigidity of the interfacial film, 
the interfacial surface tension of the oil solution 
was studied in a Cenco duNouy tensiometer, and 
it was then easy to test the physical state of the 
film. In every case a simple probe test gave the 
same results. 

The following conclusions resulted from these 
experiments: 

1) Of the monobasic acids studied only those 
having one saturated straight chain formed rigid 
interfacial films. 

2) Of the saturated single straight chain 
acids, none having less than 14 carbon atoms per 
molecule was made rigid by any of the metallic 
ions, while those having 16 or more carbon atoms 
per molecule were made rigid by all the polyvalent 
metallic ions used. The acids having 14 or 15 
carbon atoms per molecule were made rigid only 
by metallic ions of high valence. Thus, only 
thorium and lanthanum made both acids rigid. 
However, these polyvalent metallic ions probably 
make many of the acid interfacial films studied 
more condensed or more viscous, and as evidence 
may be mentioned the effect of calcium ions in 
considerably decreasing the adsorption area per 
molecule for oleic acid (to be described below). 

3) Certain hydroxy acids, ricinoleic, ricine- 
laidic, and 13-hydroxy stearic acids have inter- 
facial areas per molecule which are unaffected by 
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PH or by the presence of any polyvalent metallic 
ions. For the sake of logical presentation these 
acids will be discussed in Part IIT. 

The effect of the highly charged ions on the 
spreading of acids in oil was simply studied and 
furnished results of interest. Thorium, being 
tetravalent, showed the most marked results as 
one would expect. In the following experiments 
thorium nitrate dissolved in distilled water in a 
concentration of 3X10-* mole per liter was used. 
Such a solution will be briefly referred to here as 
thorium water and it has a pH of approximately 
3.5. 

The shorter chain-length saturated acids spread 
much more readily on the thorium water than on 
water (pH 2.0-pH 6.0) which had been acidified 
with HCI. In all cases the hexanoic acid showed 
edge loss and disk skating for w=10-* and 
w=10-* and spread and burst in an instant if 
w= 10-*. On the other hand, octanoic acid in acid 
water acted approximately the same as hexanoic 
with slightly less violent edge loss and skating, 
while on the thorium water it spread as a thick 
disk when w=10~ and as a thin colored disk 
when w=10-. Decanoic acid on thorium water 
showed the same behavior but the spreading 
effects happened more slowly, were more pro- 
nounced, and if w=10-*, the disk did not burst 
but contracted rapidly immediately after spread- 
ing as a colored film. Dodecanoic acid behaved 
similarly but the initial spreading was slower and 
the maximum diameter attained was greater. In 
all these cases the interfacial films formed were 
not rigid. On the other hand, as mentioned earlier, 
tetradecanoic and its higher homologs became 
rigid on the thorium water. 

The rapid contraction of the decanoic and 
dodecanoic acids on the thorium water was 
studied quantitatively. A disk of w=1.023 K 10-4 
dodecanoic acid weighing 13.53X10-? gram 
spread to a maximum diameter of 8 cm and then 
rapidly contracted while a disk of w= 1.096 K 10-* 
decanoic acid weighing 7.26X10-*? gram spread 
to a diameter of 6 cm and then also rapidly 
contracted as shown in Fig. 4. Now the spreading 
was due to the presence of the thorium ions; 
hence it is assumed it occurred because the 
thorium reacted with the acid molecules reaching 
the interface to form a fatty acid salt. Assuming 
that each straight chain molecule while at the 
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interface and combined with the thorium occupies 
a cross-sectional area A at the time of maximum 
spreading, and that at the maximum spreading 
area every molecule in the drop was adsorbed at 
the interface, then a simple calculation leads to a 
value of A = 12 X10~'* cm? for the dodecanoic and 
A=10X10-'* cm? for the decanoic acid chains. 
However, the cross-sectional area of each straight 
chain cannot be less than 18.5 X 10—'* cm?. Hence 
the conclusion seems inevitable that every mole- 
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cule was not adsorbed at the time of maximum 
spreading, for in the brief time required for such 


spreading far too few polar molecules could have 
escaped by edge loss (as shown by talc test) to 
have accounted for such a large difference. 
Therefore, if each adsorbed chain occupied ap- 
proximately (18 to 20) X10-'* cm?, only 60 per- 
cent of the molecules in the drop were adsorbed 
in the case of the dodecanoic acid and 50 percent 
in the case of the decanoic. Since the interfacial 
film at maximum spreading was not rigid, a more 
reasonable assumption is that the adsorbed mole- 
cules occupied an area of around (25 to 35) 
X10-'6 cm? in which case around 34 percent and 
29 percent were adsorbed, respectively. 

The subsequent rapid contraction of the oil 
disk could not have been due to a decrease in the 
area per molecule as the thorium ions bound the 
adsorbed chains together, for even though the 
area had decreased to 5 percent of its maximum 
value, the finally contracted disk was not rigid. 
Talc tests showed that the edge loss could not 
have caused this rapid contraction. Hence the 
conclusion is inevitable that the contraction was 
largely due to the dissolving back into the oil 
drop of a large part of the undissociated thorium 
salt of the fatty acid formed. Such an adsorption 
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equilibrium due to the short life of adsorption of 
the undissociated salt is quite reasonable in terms 
of results to be presented in another connection in 
Part II. 

If a hydrophil tray containing thorium water is 
first covered with an oxidized indicator oil and an 
extremely small amount of any one of the various 
fatty acids of interest is spread in the center, a 
circular monolayer of the acid is formed and its 
boundary is outlined by the colored indicator oil. 
The indicator oil can be so chosen and calibrated 
as to exert on thorium water from 2 to 15 
dynes/cm pressure on the monolayer.* The 
rigidity of this monolayer due to its reaction with 
the thorium can be observed by cutting pie- 
shaped islands out of the circular monolayer. The 
persistence of the shape of each piece after cutting 
it away can be easily seen because its shape is 
well defined by the surrounding indicator oil. By 
employing this technique, it was found that at 
the air-water interface, the straight chain satu- 
rated acids from hexanoic up all contracted and 
then became rigid even when a low film pressure 
of 2 dynes/cm was used, the time required 
decreasing from 5 to 15 minutes as the chain 
length varied from 6 to 10 carbon atoms per 
molecule. 

The fact that on the same aqueous solution oil 
disks containing the same fatty acids under a 
greater film pressure (11.3 dynes/cm) and ex- 
posed to the chemical reaction a greater length of 
time do not become rigid can most simply be 
interpreted to mean that at the oil-water interface 
either the thermal agitation of the acid molecules 
caused by the nearby oil molecules, or the 
separation of the adsorbed acid molecules by 
molecules of the solvent to form a mixed film, 
prevents cross-binding forces from forming a 
rigid film. The latter mechanism makes it possible 
to give a simple explanation of why after the 
undissociated compound was formed there re- 
sulted a considerable oil disk contraction yet did 
not form a rigid film. The thorium compound 
formed, no longer being ionized, and unbound by 
cohesional forces due to adjacent and similar 
molecules had a short lifetime of adsorption and 
so dissolved back into the oil drop until an 
adsorption equilibrium was reached. The same 
conclusion of the existence of mixed films is 
reached later on the basis of other experimental 
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facts concerning acid and amine molecules, con- 
taining 14 or more carbon atoms each. 

The effect of calcium ions on the spreading of 
oil drops containing stearic or oleic acids has been 
studied in some detail. However, it is difficult to 
do accurate quantitative work here because of 
the unknown and variable amounts of polyvalent 
metallic ions present in the distilled water and 
glass apparatus used. If chemicals to eliminate 
these metals are employed they make the results 
regular but then prevent it from being possible to 
study the effects of adding small quantities of 
calcium. If buffers are used they affect the calcium 
concentration also. The repeatability of the 
results obtained in the following experimental 
work was good and although the true polyvalent 
metallic ion concentration is not accurately 
known, they help to establish certain qualitative 
points of interest. 

Alkaline solutions 0.02 N in NaOH were used 
to which had been added 10-* molal KHCQ3. 
Then oil disks containing w=0.890 X 10- stearic 
acid were spread on this solution, and the effect 
on the diameter vs. time curves of adding various 
amounts of calcium chloride was studied. In 
Fig. 5 are given the results obtained with oil 
disks weighing 15.5210-? gram, each curve 
being for a different calcium concentration (by 
calcium concentration is meant only the ap- 
parent concentration added since the value 
initially present in the water is not known). In 
Fig. 6 is given the area per stearic acid molecule 
as a function of apparent concentration. The 
dotted portion indicates the effect of eliminating 
the last of the polyvalent metallic impurities by 
chemical methods described in the next section. 

Since all the changes caused by thecalcium were 
made by less than 200 parts of calcium per 10%, 
it is of interest to calculate approximately the 
effect of the HCO;~ ions and OH~ ions on the cal- 
cium ion concentration. Since 2X 10-? gram mole- 
cule per liter of Ca(OH). will dissolve in one liter 
of water, it is easy to show that the concentration 
of hydroxide ions present will not complete any 
precipitation of as little as 200 parts of calcium 
per 108. However, since the second dissociation 
constant of KHCO; is 6X10-", a concentration 
of 10-* molal in HCO;- will furnish approximately 
2.4X10-4 molal CO3-- ions. [Cat+][CO;-—] 
=0.9X10-8; the calcium ion concentration 
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present in solution cannot exceed 150 parts per 
10%. Hence the KHCOQO; concentration used may 
somewhat affect the calcium concentration but 
certainly not to a serious extent. 

Returning to Fig. 5 and Fig. 6, certain con- 
clusions of interest can now be drawn. In the case 
of curve A the first acid molecules reaching the 
water have been close packed by the calcium 
ions, and have become rigid so preventing any 
further spreading being caused by the later arrival 
of acid molecules yet in the oil drop. On the 
other hand, in curve D, since much less calcium 
was present, every molecule was able to adsorb 
at the interface and the film formed was liquid. 
Here A = 105 X 10~-"* cm?. The further elimination 
of the calcium and other polyvalent ions present 
led to the value A =125X10~-'® cm*. The slight 
but definite decrease in diameter after reaching 
the maximum spreading was in part due to a 
slow rate of edge loss, which sulfur powder tests 
showed still continued, and the greater part of 
the contraction due to the condensing effect of 
the traces of calcium. In this case if the contrac- 
tion was followed long enough it was found that 
A approached the value (15 to 17) 10~-'® cm’. 
From curve D it is seen that about 6} minutes 
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were needed for all acid molecules to reach the 
interface. The longer time needed at C to reach 
maximum diameter is presumably due to the 
slow rate of migration of the later arriving mole- 
cules through the partially condensed interfacial 
film. The more rapid rate of decrease of curve C 
than of curve D is almost entirely due to the 
increased rate of condensation of the adsorbed 
film made possible by the higher concentration of 
calcium ions. 
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9. RESULTS OF ELIMINATING THE POLYVALENT 
METALLIC IONS 


Because of the important influence shown by 
polyvalent metallic ions on the spreading of the 
oil disks, and since the best distilled water 
obtainable has enough of such metals when used 
in a Pyrex tray to have appreciable influence on 
the results, in some experiments resort was made 
to chemicals to precipitate or render less effective 
these metals. Like Langmuir and Blodgett,'’ use 
was made of a solution of 0.02 N NaOH and 
0.02 N NazCOs, and it was found that such a 
solution aged overnight helped matters con- 
siderably. Using still higher concentrations of 
NazCQ; is of little advantage since adding the 
Na2CO; also brings with it an increase in metallic 
impurities. To illustrate the results obtained 
using 0.02 N NaOH alone, values of A of from 
(90 to 105) X 10-'* cm? were obtained for stearic 
acid while using 0.02 N NaOH+0.02 N Na2CO; 
solution resulted in values of from (105 to 110) 
X10-'* cm?. Following Langmuir and Schaefer® 
use was made of about 610-4 gram of Calgon 
per liter of the above solution to remove poly- 
valent metallic ions. The result was that A 
ranged from (120 to 125) X10-'* cm? (see dotted 
portion of Fig. 6). Such an alkaline solution was 
adopted as a standard and for convenience will be 
called Calgon water. 

In Table I are listed the values found for A 
employing Calgon water. Differences in the 
values between any two acids greater than +5 
percent need not be taken seriously since it is 
entirely possible that they are due to error or to 
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accidental metallic impurities. In the case of the 
dodecanoic and tridecanoic acids, to minimize 
the error due to edge loss, A was measured using 
a weight concentration of w=2X10-‘, thus 
making 7 short. 

Summarizing these results: with the exception 
of the long chain hydroxy acids, A is always in 
the range 100 to 125 square angstroms for all 
acids whose chain lengths are greater than 14 
carbon atoms per molecule. The apparent pro- 
gressive decrease in the area per molecule with 
chain length for values of N less than 14 is so 
rapid that one infers the existence of a corre- 
sponding decrease in the lifetime of adsorption, 
but the error in the value of A due to the effect of 
edge loss is unknown and hence this cannot be 
considered the only logical conclusion. Since 
most of these molecules occupy, when close- 
packed at the water-air interface, an area of 
around (20 to 35) K10-'* cm?, it is apparent that 
at the oil-water interface they occupy an area of 
from 33 to 6 times as great. 


10. Proor oF 100 PERCENT ADSORPTION 


An experimental technique was developed to 
directly prove the correctness of Langmuir and 
Blodgett’s conclusion regarding the infinite life- 
time of adsorption of stearic acid on alkaline 
water and where necessary to extend the result to 
other adsorption problems. When the oil disk had 
spread to its maximum diameter, a device named 
for convenience a “‘multiple dropper’’ was used to 
cause droplets of “piston oil’”’ (usually either 
castor oil or oleic acid) to touch the surface of the 
water simultaneously at 16 equally spaced points 
arranged in a circle external to but concentric 
with the oil disk. The piston oil spread simul- 
taneously from the 16 points of contact and 
squeezed the disk from 16 directions at once, so 
compressing it under a pressure equal to the 
equilibrium spreading pressure of the piston oil 
used. With careful centering the contracting disk 
was not driven off to one side against the 
hydrophil tray wall. The oil disk contracted 
under the equilibrium pressure of the piston oil 
and the area per molecule decreased to a smaller 
value. The advantage of such a device is that the 
oil disk was permitted to spread to its maximum 
area in the absence of any piston oil or external 
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polar monolayer, thus minimizing the danger of 
contamination of the interfacial layer. After the 
spreading was complete the disk compression was 
so rapid that a measurement of the new value of 
the area per molecule was obtained before any 
appreciable error from contamination had re- 
sulted. The method has the disadvantage of 
inflexibility if one wishes to obtain a series of 
measurements of A as a function of external 
piston oil pressure. 

To illustrate the results obtained, the experi- 
ment was performed with w=1X10-° stearic 
acid on Calgon water. This disk usually spread to 
a diameter of around 6 or 7 cm. Then the castor 
oil caused that disk to contract immediately to a 
diameter of around 3 cm corresponding to 
A=(27 to 32) X10-'* cm? and in 2 or 3 minutes 
more it had decreased to a minimum diameter of 
2.6 cm corresponding to A=(18 to 24) x10-'® 
cm*. In the course of another 3 to 5 minutes it 
stayed close to that minimum diameter and 
thereafter it gradually increased with time. 
Simple tests demonstrated that this subsequent 
increase was due to the dissolving of the castor 
oil in the oil disk. It was not due to the effect of 
the alkali on the castor oil since it was found that 
the piston pressure of the castor oil, initially 20 
dynes/cm on such an alkaline solution, decreased 
to 19 dynes/cm in 35 to 40 minutes. The effect 
was reversible to a certain extent for if the castor 
oil was scraped from the water surface the disk 
spread out again, while upon again adding castor 
oil it contracted to the same minimum area. If 
oleic acid was used for the piston oil its initial 
pressure on alkaline water was around 38 
dynes/cm, decreasing to 31 dynes/cm in 35 to 40 
minutes. This higher pressure caused the oil disk 
to contract at once to a value corresponding to 
(15 to 19)10-'® cm?, but immediately after 
that the disk continued to contract ; in short, this 
pressure caused collapse of the interfacial 
monolayer. 

A large number of such measurements resulted 
in a mean value for the compressed stearic acid 
interfacial film of A = 20.6 X 10-!* cm?+5 percent. 

When a solution of oleic acid in light petrolatum 
was compressed in the same way, the value of A 
under the pressure of castor oil decreased im- 
mediately to (30 to 34)x10-'* cm?*, which is 
nearly the same value found at the air-water 
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interface on acid water at that same film pressure. 
However, using oleic acid for the piston oil 
caused the disk to contract further to a minimum 
area corresponding to A=20.5X10-'® cm?+5 
percent. Hence, the higher piston oil pressure of 
oleic acid was needed to close-pack the oleic acid 
molecules, presumably as a result of its cis- 
configuration. 

From these experiments it is evident that all of 
the stearic acid or oleic acid molecules originally 
present in the oil drop had migrated and re- 
mained at the oil-water interface. From the fact 
that all of the other long chain acids of Table I 
showed the same range of values of A on Calgon 
water, one concludes that all these long chain 
acids have infinite lifetimes of adsorption on 
alkaline water. 


11. ORIGIN OF LARGE VALUES OF A 


Langmuir and Blodgett’ inferred from their 
work on stearic acid in Squibb’s petrolatum that 
their large value of A =90X10-'* cm? was due to 
a considerable hydration of the carboxy! group of 
the molecules. From experiments on compressing 


such disks they concluded that the relation 


TABLE I. Area per molecule of fatty acids at the interface 
between Calgon water and oil. 








No. oF CARBON 
ATOMS PER 
MOLECULE 


A X10'6 
ACID 





Dodecanoic (lauric) 12 
Tridecanoic 13 
Tetradecanoic (myristic) 14 
Pentadecanoic 15 
Hexadecanoic (palmitic) 16 
Heptadecanoic (margaric) 17 
Octadecanoic (stearic) 18 
Nonadecanoic 19 
Eicosanoic (arachidic) 20 
Octadecenoic (oleic) 18-cis 
Octadecenoic (elaidic) 18-trans 
13-docosenoic (erucic) 22 





12-hydroxy-9 octadecenoic 18-cis 
12-hydroxy-9 octadecenoic 


13-hydroxyoctadecanoic 





Chaulmoogric 

Di-n-heptylacetic 

Methyltetradecylacetic 

Propyldodecylacetic 

Phenylstearic 91 
Ethylphenylstearic 103 
Xenylstearic 114 
Tetrahydronapthylstearic 117 








* Edge diffusion loss too great for an accurate measurement of A. 
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F(A —A»o)=hkT held true for such an interfacial 
film and that Ay=57X10-'* cm?. It seems quite 
unlikely that enough water of hydration adheres 
to each carboxyl group so tenaciously as to cause 
A» to be that much greater than the value 
22X10-'*, particularly when it is remembered 
that the interfacial monolayer is under a com- 
pression of between 11 and 14 dynes/cm due to 
the large negative spreading coefficient F, of the 
oil disk. If the water of hydration were re- 
sponsible, the castor oil piston pressure experi- 
ment described above would cause the total film 
pressure to rise from between 11 and 14 dynes/cm 
to between 19 and 20 dynes/cm; and since the 
first pressure would not brush off all the water of 
hydration, one would not expect the second 
pressure to do so. Yet the fact is that A under the 
higher pressure is 20.6 X 10-'* cm”. 

In the above experiment with the castor oil, a 
definite time interval was found to occur before 
the interfacial film was compressed to the 
minimum area corresponding to close-packing. 
This suggests that the acid molecules are sepa- 
rated by hydrocarbon molecules of the oil and 
that under the external pressure of the piston oil 
the hydrocarbon molecules were forced up and 
out of the interfacial film, thus eventually per- 
mitting close-packing of the acid molecules. If 
this were true, then the results of Table I would 
mean that all these long chain acids when 
adsorbed form 2-dimensional complexes or solu- 
tions with hydrocarbon molecules from the oil. 
The large value of A and the fact that it is nearly 
the same for all acids are then quite under- 
standable. The existence of oil molecules in the 
adsorbed film is not an unimportant point since 
it is reasonable to assume that the hydrocarbon 
molecules present between the adsorbed acid 
molecules will have a certain amount of bonding 
with them. If this were not the case, thermal 
agitation would eventually cause these hydro- 
carbon molecules to be worked out of the 
interfacial film and eventually they would be 
replaced by adsorbed acid molecules, thus leading 
to a system of lesser energy. Entirely independent 
evidence for this conclusion, will be presented 
later in Part II. 

The above-mentioned value of Ap=57X10-"® 
cm? obtained by Langmuir and Blodgett may be 
given a simple molecular interpretation: The 
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electrical repulsion between adjacent dissociated 
stearic acid molecules acts counter to the net 
surface tension force of —11 to —14 dynes/cm 
and hence more than one hydrocarbon molecule 
can insert itself and be weakly bound between 
each pair of acid molecules. Their difference 
(90 —57=33) X10-'* cm? suggests that there are 
two hydrocarbon molecules per acid molecule in 
the complex film, and that under the action of 
small external pressures, the disk readily contracts 
until one of the two hydrocarbon molecules is 
forced out of the film. However, the remaining 
hydrocarbon molecule is much more strongly 
held in the film, and a greater external pressure 
than was used in their experiments is needed to 
force it out. One would then be led to the 
conclusion that the true area per acid molecule is 
57 X10-'* cm?. But it has here been shown that a 
greater external pressure will cause the area to 
decrease to approximately 21 10-"* cm?. It will 
be noted that (57 —21=36) X10~-!® cm? which is 
in reasonably good agreement with the above 
deduced value of 33X10-'* cm? for the area 
occupied by each of the hydrocarbon molecules in 
the complex film. 

The above hypothesis as to the relative abilities 
of 1 : 1 and 1 : 2 film complexes to stand external 
pressure, and also the resulting inferred area per 
hydrocarbon molecule are in good agreement with 
the results of the recent study by Davis, Krahl 
and Clowes?® of two-dimensional complex films 
formed at the air-water interface by various 
mixtures of sterols and aromatic hydrocarbons. 


12. BEHAVIOR OF MIXED FILMS OF OLEIC AND 
STEARIC ACIDS 


When drops of oil containing various pro- 
portions of both oleic and stearic acids were 
spread on Calgon water, the observed area per 
molecule was the same as for the pure acids; i.c., 
12010-'® cm*. This was to be expected since 


they each occupied the same area when pure and 


also, being so far, apart the acid molecules showed 
little interaction. When drops of the same mixed 
acids were placed on 0.02 N NaOH to which 
2.5X10-* molal CaCle was added, the results 
given in Fig. 7 were found. The method of 


29 W. Davis, M. Krahl and G. Clowes, J. Am. Chem. Soc. 
62, 3080 (1940). 
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obtaining the values given requires some expla- 
nation. The 100 percent oleic and the 50-50 
percent oleic-stearic mixture spread on the 
calcium water without trouble. Using the multiple 
dropper with castor oil as described above, it was 
possible to prove there was no serious error from 
edge diffusion and that in each case every acid 
molecule in a drop had been adsorbed at the 
interface. Hence within the experimental error 
indicated by the scattering of points on the 
graph, the value of A, the area per acid chain, 
was really the same for the 100 percent oleic or 
the 50-50 percent oleic-stearic mixture. In the 
case of the 25 percent oleic 75 percent stearic 
mixture, the disk required from 3 to 5 times 
longer to spread to its maximum diameter but 
eventually spread to an area per molecule of 
(16 to 19) K10-"* cm?. Allowing for the edge loss 
of so long an experiment, the correct value of A 
was assumed to have been around (18 to 22) 
X 10-'6 cm?, indicating close-packing. The spread- 
ing area of 20X10-'* cm? indicated on the graph 
for the 100 percent stearic acid is a reasonable 
assumption although the actual spreading is 
negligible, for that lack of spreading is due to the 
formation of a rigid interfacial form in which A 
must have a value between (18 to 22) X107!* cm?. 

Before interpreting these results, it is of 
interest to note that the calcium ions cannot 
cause close-packing in interfacial films containing 
pure oleic acid. In fact, thorium ions also are 
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unable to close-pack interfacial films of oleic acid 
at either the oil-water or the air-water interface. 
The value of A =57X10-'* cm? is much greater 
than that to be expected from the results of 
studies at the air-water interface for there under 
film pressures of 11 to 14 dynes/cm, A hasa value 
of around (32 to 34)X10-'* cm*. It might be 
thought that this value of 57 X 10-'® cm?, is due to 
the adsorption of the double bond of oleic acid at 
the oil-water interface. However, Marsden and 
Rideal*® have shown from their studies of oleic 
and elaidic acids at the air-water interface that 
the lack of close-packing for oleic acid is due to 
the cis-configuration. Only that conclusion made 
understandable their result that elaidic acid, the 
trans-isomer, formed much more condensed films. 
In agreement with their result, experiments 
performed here using oleic and elaidic acids at the 
air-water interface and employing thorium water 
showed that whereas the elaidic acid was made 
rigid by these highly charged ions, the oleic acid 
was not. 

The use of the hypotheses of Section 11 con- 
cerning the formation of mixed film complexes 
also permits a simple molecular explanation of 
these observations on oleic-stearic acid films. A 
calcium ion reacts with two dissociated acid 
chains, and because of the strong electrical 
binding effect of a divalent metallic ion, the 
resulting molecules are little dissociated. Hence, 
the acid chains are no longer repelled by electrical 
forces and hence, as in Section 11, the film is 
acted on only by the compression effect of the 
surface tension forces (from —11 to —14 
dynes/cm). This forces out some of the hydro- 
carbon molecules of the original mixed acid- 
hydrocarbon film until only one hydrocarbon 
molecule is left per acid chain. The result is a 
complex made up of four-leaf-clover-like aggre- 
gates of two acid chains and two hydrocarbons. 
Apparently the cross-binding forces due to Ca 
ions are not sufficient to force all the hydrocarbon 
molecules out of the film to form a rigid interfacial 
monolayer. The area per acid chain can reasonably 
be taken to be around (30-35) X10-'* cm?2, since 
that is the area per oleic acid molecule at the air- 
water interface in acid water at a compression of 
11 to 14 dynes/cm. This reasonable assumption 


90 J. Marsden and E. Rideal, J. Chem. Soc. London, 1163 
(1938). 
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leads to the conclusion that the area per hydro- 
carbon molecule still left is (57-—33=24) 10-16 
cm’, which is a good approximation to the area 
occupied by a straight chain paraffin hydro- 
carbon molecule or by many types of aromatic 
hydrocarbons. 

However, the 0-25 percent mixture causes 
interfacial film freezing because a large enough 
percentage of these salt molecules are calcium 
stearate for there to exist considerable cross- 
binding due to the added effect of the van der 
Waals forces between the stearic acid hydro- 
carbon chains. In the case of 25-50 percent 
mixture such freezing apparently goes on in small 
patches between which places that are not close- 
packed are left permitting the newly arriving acid 
molecules to reach the interface. Hence the 
adsorption is slow but complete. 


13. Errect oF HiGH Actp CONCENTRATION 


Oil disks containing more than w= 10- of long 
chain acids spread even on acid water and the 
mechanism and quantitative relations involved 
will be discussed in Part II. 


14. THE SPREADING OF AMINES WHEN IONIZED 


The long chain primary normal amines, like the 
long chain acids, spread to a remarkable extent 
when ionized and very much less when not. The 
behavior of solutions of heptadecyl amine in 
mineral oil is typical. In Fig. 8 is plotted the 
diameter of spreading as a function of pH of the 
water for oil disks weighing 1.7 X10-? gram and 
containing the amine in a weight concentration 
w=1.204X10-*. The solutions were made up 
with hydrochloric acid, KHCO; or NaOH. It is 
evident that for pH values exceeding 6 there is 
practically no spreading, while it is complete for 
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values more acid than pH 2. The small diameter 
attained for pH greater than 6 is precisely the 
same as that of the same volume of oil free from 
the presence of amines. This behavior is quite 
general for all amine concentrations below 
roughly 10-*. This result is qualitatively in 
agreement with that obtained by Peters" in 
connection with his study of the effect of PH on 
the benzene-water interfacial tension lowering 
caused by hexadecyl amine dissolved in benzene. 
The spreading behavior of amines on alkaline 
water for higher values of w than 10-* will be 
discussed in Part II in connection with the 
discussion of the spreading of acids in similar 
concentrations on acid water. 

As is well known, a primary amine such as 
heptadecyl amine reacts in dilute aqueous 
hydrochloric acid to form heptadecyl amine 
hydrochloride, which then dissociates to form a 
long chain positively charged complex ion and a 
chloride ion. Or in general, 


H H 4 
R:N:+HtcCl=|]R:N:H] +Cr. 


H H 


Hence the spreading of the oil disk on acid water 
was caused by the migration to the interface of 
heptadecyl amine followed by a dissociation 
reaction and adsorption which was complete if 
PH was less than 3. 

Since dissociated acid molecules were found to 
have infinite lifetimes of adsorption, one would 
expect the dissociated amine molecules to behave 
likewise. Such was found to be the case when the 
multiple dropper technique was used to close- 
pack the fully spread oil disk containing hepta- 
decyl amine, and a contraction to an area per 
molecule of (18 to 20) X10~-'* cm? was observed. 

Measurements of the area per adsorbed mole- 
cule A carried out by methods entirely similar tc 
those already described for acids, but using 
distilled water 10-2? N in HCl, furnished the 
results of Table II. Only a brief study was made 
of the effects of various negative polyvalent ions 
on the value of A, and although the ions CO; 
PO,-—~ and Fe(CN).~~—~~ did not cause rigidity 
of heptadecy! amine at the oil-water or air-watet 
interface, it is possible that some condensation of 
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the film resulted. The molecular meaning of the 
data of Table II cannot be stressed too much for 
that reason. However, it appears that the mutual 
repulsions between the dissociated amine hydro- 
chlorides is less than between the corresponding 
dissociated acids, and A increases linearly with 
N, the number of carbon atoms per molecule. 

The rate of edge loss was found to be greater 
for a primary amine than for the corresponding 
single chain aliphatic acid, and as a result the 
error in measuring A due to edge loss was more 
serious. The observed decrease of A with de- 
creased chain length may be partly due to this 
cause. 

The time, 7 y, for the disk to reach its maximum 
diameter was found to be approximately the 
same for heptadecyl amine as for heptadecanoic 
acid at the same molal concentration. Hence the 
rate of diffusion of these two types of polar 
molecules through the oil is approximately inde- 
pendent of the nature of the polar group. 

The fact that the value of A is so much greater 
than the true molecular cross-sectional area 

18.5 to 24) X10-'* cm? as in the case of the acids 
can also be interpreted to mean that the adsorbed 
amine film at the oil-water interface is a 2- 
dimensional solution in oil. 

A study was made of the spreading of oil drops 
containing the secondary and tertiary amines, 
di-n-octyl amine and di-n-amyl amine, as well as 
tri-n-octyl amine and tri-y-amyl amine. The 
octyl amines behaved like the primary amines 
with regard to the effect of pH but the oil drops 
spread with greater rapidity and showed even 
more edge diffusion. The amyl amines either 
spread so rapidly as to appear to burst, or simply 
skated over the surface of the water while 
exhibiting intense edge diffusion. Since these 
triamines have either 24 or 15 carbon atoms per 
molecule, respectively, it is again evident that 
the rate of edge diffusion is not controlled by the 


TABLE II. Area per molecube of amines on acid water. 


NUMBER OF CARBON 


AMINE ATOMS PER MOLECULE A X10'6 cm? 





Octadecanoic amine 18 
Heptadecanoic amine 17 
Hexadecanoic amine 16 
Tetradecanoic amine 14 
Octadecenyl and 

octadecendienyl amine 18 








total molecular weight but the length of the 
longest chain present. 

Oil solutions of diphenyl amine and triphenyl 
amine were tested similarly but neither showed 
any spreading variation with pH, and in fact 
spreading occurred only if w considerably ex- 
ceeded 10-*. On the other hand, dibenzyl amine 
and tribenzyl amine showed more spreading on 
acid than on alkaline water, but the amount of 
spreading was negligible unless w also consider- 
ably exceeded 10-*. These results agree with the 
well-known dissociation rule that a phenyl group 
attached directly to the nitrogen atom of an 
amine group cuts down its degree of dissociation, 
while the further away the point of attachment of 
the phenyl group from the nitrogen atom, the 
less it affects the dissociation constant. The fact 
that the pH has effect on only the benzyl amines 
is therefore understandable. The difficulty in 
causing disk spreading is merely due to the fact 
that the longest chain length present in the 
molecule is still so short that edge diffusion is the 
most pronounced spreading effect unless high 
concentrations are used. 
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A comprehensive thermodynamic theory is developed 
for the spreading of any liquid or solid b over the surface 
of any liquid a. Spreading is considered to occur by two 
types of processes: (A) duplex film or D spreading, and 
(B) non-duplex or M spreading. A duplex film is always 
unstable, and transforms into a non-duplex film and a lens. 
Every non-duplex film of known structure is a monolayer 
(M). Either type of film may spread if its formation in- 
volves a decrease of free energy. The free surface energy 
of mercury at 20°C is 476 erg cm™, and the spreading of 
water or any organic liquid on the surface as either a 
duplex film or a monolayer, gives a decrease of free surface 
energy, so all such liquids spread on mercury. On water 
all organic liquids spread as monolayers; an organic solid 
also spreads as a monolayer unless it is too non-volatile in 
the two-dimensional system. A liquid } will spread over 
the surface of a liquid @ as a duplex film if the initial 
spreading coefficient Stia= Wa— We, is positive, that is if 
the work of adhesion W,4 between the two liquids is greater 
than the work of cohesion (Wc,) in b. If b is saturated with 
a, then Syya is designated as the semi-initial spreading 
coefficient, which has almost the value of the initial coef- 
ficient. If the liquids are mutually saturated the final 
spreading coefficient, Spy, is always negative. Thus a 
liquid 6 will never spread over its own monolayer, and if 
the whole surface of water in any vessel is covered by a 
duplex film, it is certain that a monolayer has not spread 
over any part of it. After a definite period, the length of 
which is dependent on many factors, a duplex film begins 
to segregate into a monolayer, patches of thicker duplex 
film, and lenses, and finally into a monolayer and a lens. 
The film pressure z, of a monolayer of 6 in equilibrium with 
a lens of b is equal to the difference between the semi- 
initial and the final spreading coefficient, or me= Sw 


— Soar. Since Soya: is always negative, the equilibrium 
film pressure is always greater than the semi-initial coef- 
ficient. Hydrocarbons of not too high molecular weight 
are found to spread on water as either duplex films or 
monolayers. Thus many nonpolar oils spread on water to 
form either type of film. Benzene, which according to one 
of the most prominent theories does not spread at all, spreads 
readily, and exhibits the following values, which are almost 
the same as those for isopentane: initial coefficient 8.94, 
semi-initial 8.90, and equilibrium film pressure 10.6, in 
erg cm. Carbon disulfide (nonpolar) and methylene 
iodide (highly polar) have initial coefficients —7.3 and 
— 23.8, and therefore do not give duplex films, but they 
vaporize in the two-dimensional surface, and give equi- 
librium film pressures of 2.3 and 0.6 dyne cm~}, respec- 
tively. Presumably their monolayers are gaseous. In 
general for liquids which spread as duplex films, it takes 
less work (i.e., the increase of free energy is less) to pull a 
liquid away from its equilibrium monolayer than to pull 
it away from itself. For benzene the difference is 1.7 and 
for n-heptyl alcohol, 5.4 erg cm~*. An oil which is non- 
spreading in the sense of duplex film formation, may be 
caused to spread as a duplex film by the addition of an 
oil of sufficiently high initial spreading coefficient. This 
also increases the equilibrium film pressure of the mono- 
layer. The oils added for this purpose are, in general, of 
the polar-nonpolar type, such as stearic acid. With such 
mixtures the values of the thermodynamic spreading coef- 
ficients, or the free energy levels, are highly dependent 
upon the nature and the concentration of the oil added. 
If the aqueous subphase is alkaline, stearic, or another 
similar acid is converted into a salt, and if the positive 
ions have a charge of two or higher, this has a pronounced 
effect upon the spreading coefficients and the film pressure. 





SYMBOLS 


A, Helmholtz free energy. 

a, The liquid phase a. 

a’, Phase a saturated with b. 

b, The liquid phase b. 

b’, The phase 6 saturated with a. 

, Surface tension; or free surface energy for 1 cm, 
vs, free surface energy of the film covered surface, and 
Ye, equilibrium value of ;. 

Yab, Yow OF yi, Interfacial tension or free energy. 

d, Duplex film, D spreading as a duplex film. 

5, Non-exact differential. 

E, Internal energy. 

e, Internal energy for unit area (molecular potential 
energy). 


F, Gibbs free energy. 

f, (subscript) denotes a film. 

f, free surface energy per sq. cm. 

G, Gravitational potential, g, acceleration of gravity. 
h, enthalpy per cm’. 

L, liquid. /, (subscript) denotes a lens. 

m, Monolayer or mass. M, Spreading as a monolayer. 
o, Oil. 

p, Pressure in dyne cm™. 

x, Film pressure in dyne cm~. 

p, Density, p, density of vapor phase. 

S, Entropy. 

~, Summation of terms. 

Soa, Initial spreading coefficient of b on a. 
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Sv’ sa, Semi-initial spreading coefficient of b on a. 

Swja, Final spreading coefficient of mutually saturated b 
on a. 

T, Temperature, Kelvin scale. 

v, Volume. 

Wc», Work of cohesion for the liquid 6. 

Wa. Work of adhesion for the interface ab. 

w, Water. 


1. INTRODUCTION 


OME of the most fundamental relations of 
surface chemistry are those involved in the 
spreading of a liquid phase over a subphase to 
form a duplex film and a monolayer. Nevertheless 
there is no branch of the subject which is in a 
more confused state in both the earlier and the 
more recent literature. This is probably due to 
the fact that there have been not only contra- 
dictory theories of the spreading of duplex films, 
but also their further transformation to mono- 
layers has not been treated from the thermo- 
dynamic standpoint, nor from that of statistical 
mechanics. 

The term duplex, as used here, does not refer to 
any application of a double film theory to a 
monolayer, in which case the upper and lower 
surfaces of a single molecular layer are treated as 
having almost independent interfaces. According 
to the definition of the term as used here a duplex 
film is a phase so thin as to have properties which 
are practically independent of the force of 
gravitation, but which is sufficiently thick to 
allow the interfacial energy, both free and total, 
of each of its interfaces to be the same as when 
both liquids are present in bulk, so that each 
interface is independent of the other. In this 
particular sense of the term a monolayer cannot 
be a duplex film. 


2. THEORIES OF SPREADING 


All of the divergent views concerning spreading 
may be considered as included in the three 
theories outlined below. 

1. According to the theory of Lord Rayleigh,' 
as based upon the work of Marangoni,? Van der 
Mensbrugghe,* and Quincke,‘ (A) All liquids 


‘Lord Rayleigh, Scientific Papers, Vol. III (1887-92), 
pp. 351, 397, 513, 562, 572. 

* Marangoni, Pogg. Ann. 118, 348, 1871 (1865). 

* Van der Mensbrugghe, “Essai sur la Theorie Mécanique 
de la Tension Superficielle,”” Bull. Acad. roy. Belgique, (3) 
9, 12 (1885). 

* Quincke, Wied. Ann. 2, 145 (1877). 
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spread on water: (B) Lenses of a liquid do not 
form upon water unless the surface of the water 
is contaminated by some foreign substances. 
The Neumann triangle of forces is imaginary and 
there cannot be equilibrium between a lens and 
its own film or water. ‘“‘No such equilibrium, and 
no such triangle, is possible if the materials are 
pure; when it occurs it can only be due to the 
contamination of one of the surfaces.” 

2. The theory of Langmuir® is extremely 
simple, and on this account has exerted a wide 
influence. It may be stated as follows: Liquids 
whose molecules contain polar groups spread on 
water; if polar groups are not present the liquid 
will not spread. In a recent paper Langmuir 
applies this theory to the spreading of hydro- 
carbons. Thus he states: “‘Pure hydrocarbons 
and other substances whose molecules are wholly 
hydrophobic do not mix with water, nor do they 
spread, but remain on the surface as a globule. 
The only substances which spread are those 
whose molecules contain both hydrophylic and 
hydrophobic parts.”’ Again on p. 182: “‘A pure hydro- 
carbon does not spread.”’ The italics have been intro- 
duced to emphasize the points under discussion. 

3. In the present paper a comprehensive 
thermodynamic theory of spreading is presented. 
This is not restricted to oils on water, but is a 
general theory which includes the spreading of 
liquids on the surface of solids, and of solids on 
the surface of liquids. 

The relations of the phenomena involved are 
much more clearly understood if it is considered 
that there are two types of spreading. These are: 
(A) Duplex (or D) spreading: a liquid spreads as 
a duplex film. The thermodynamic relations 
indicate that this is unstable. (B) Non-duplex 
(or M) spreading: a non-duplex film may spread 
directly from a three-dimensional liquid or solid, 
or it may be formed by the transformation of a 
duplex film. All non-duplex oil films on water, 
whose structure is known, are monolayers. 

The thermodynamic theory developed in this 
paper establishes the conditions which determine 
the spreading of material, initially in the form of 
a liquid or a solid phase, over the surface of a 
liquid. It applies also to the spreading over the 


51. Langmuir, J. Am. Chem. Soc. 39, 1864 (1917). 
61, Langmuir, Cold Spring Harbor Symposia on Quanti- 
tative Biology, Vol. 6 (1938), pp. 173, 182. 
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Fic. 1. Free surface energy (f) diagram for the spreading of duplex films and of monolayers on the surface of water. 
The free surface energy of a duplex film is the sum of that for the interface (ywo) and for the surface of the oil 
(fa=Ywot7o). The lower of the two energy levels of a duplex film represents oil saturated with water, while the upper 
level is that when no water is in the surface of the oil. S designates the semi-initial spreading coefficient Sy-/4 which may 
be positive (arrow pointing downward) or negative (arrow pointing upward). S’, which represents the final coefficient 
(Sv-sa’) 1s always negative. Obviously if the free energy of a duplex film is above that of water, it will not form: if below 
that of water it may spread, but is always unstable with reference to the formation of a monolayer (and a lens). For the 
spreading of water on an oil the free energy levels for the duplex film are above that for water (usually by 40 to 110 
erg cm~?) so water does not spread on an oil as a duplex film. Thus, except for the greater height of the duplex level the 
diagram is like that for methylene iodide. : 





surface of a solid, insofar as this can be con- 3. FREE ENERGY DIAGRAM FOR THE SPREADING 

sidered to have a uniform, instead of a periodic, OF OILS ON WATER 

surface field. The amount of work which must be done to 
The only condition which determines whether form one sq. cm of the plane surface of water is, 

spreading of a particular type may or may not at 20°C, 72.8 ergs. Thus in a diagram (Fig. 1) the 

occur, is the change of free energy involved. surface of water may be represented as existing 
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at a free energy level of 72.8 erg cm~*. Any oil 
which can reduce this free energy level by spreading 
to give any particular type of film, may spread to 
give such a film. If the particular type of oil film 
lies at a higher free energy level than that of 
water, it will not spread spontaneously, and if it 
could be spread by any mechanical means, would 
be unstable. 

It is apparent from the figure that all four of 
the liquids represented spread on water to form 
monolayers. The spreading coefficient, that is the 
decrease in free energy involved when a lens of 
the oil spreads to form its equilibrium monolayer, 
is given by the value of the film pressure z,. The 
values in erg cm are as follows: n-heptyl 
alcohol, 44.2 ; benzene, 10.7 ; carbon disulfide, 2.3, 
and methylene iodide, 0.6. Thus this diagram 
appears to justify that part of Rayleigh’s theory 
which considers that all liquids spread on water. 
However the equilibrium film pressure of a non- 
polar liquid oil of very high boiling point, such as 
octadecane is presumably very low (at 30°C). 
The film pressure of a solid n-paraffin of very high 
molecular weight is extremely small and would be 
difficult to detect. 

With respect to their spreading as duplex films 
all liquids fall into two classes, those for which 
the free energy of the duplex film is (1) below, 
and those for which it is (2) above, that of the 
subphase. The free energy level of any duplex 
film is equal to the sum of its two interfacial 
energies, or f=YotYwo. Thus for benzene, satu- 
rated with water, f=28.8+35.0=63.8, while for 
dry benzene f=28.9+35.0=63.9. Thus the 
duplex film does not exhibit a single free energy 
level, but its different levels are close together. 

The decrease of free energy for the spreading of 
a duplex film of benzene on water is 9.0 erg cm~, 
and for isopentane 9.4 erg cm~?, which indicates 
that not only aromatic compounds, but also 
lower n-paraffin oils, spread readily on water to 
form duplex films. 

While carbon disulfide and methylene iodide 
spread as monolayers, it is evident that they do 
not form duplex films. 

The relations between spreading coefficients 
and spreading involves the thermodynamic re- 
lations presented in the next section. 

The theory of (A) duplex film spreading as 
presented here is an elaboration of the prelimi- 
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nary theory of the writer,’ which was based on 
the earlier considerations of Dupré* and Hardy.’ 

Non-duplex, presumably monolayer, spreading 
is found to occur with all liquids which are not too 
non-volatile with respect to the two-dimensional 
system which consists of the clean surface of the 
liquid used as the subphase. To this extent there 
is agreement with (A) of Rayleigh’s theory. It is 
found, however, that (B) of Rayleigh’s theory is 
not valid, since a lens of oil is found to be stable 
in the presence of its own monolayer. 

The lower hydrocarbons are found to spread 
both as duplex films and monolayers, while those 
of higher (but not too high) molecular weight give 
monolayers alone. 


4. THERMODYNAMICS OF THE SPREADING OF 
DuPLEX FILMS 


The maximum work (6W) done by a surface 
when its area (a) is increased by do is 


5W = —ydo+ pdr. . (1) 


The value of the free energy of Gibbs is defined 


by the equation 
F=E—ST+pv (2) 


and the Helmholtz free energy by 


A=E-ST. (3) 
It follows that 
F=A+pv (4) 
and 


dF=dA+pdv+udp. (5) 


If a process is reversible and isothermal, TAS=Q 
is the heat adsorbed and 


dA=-—dwW, (6) 


dF= —dW-+ pdv+vdp. (7) 
From (7) and (1) 
dF=ydo+vdp, (8) 
so that at constant pressure and temperature 


(OF /00)p,r=¥ (9) 


and for a saturated surface y is a function of p 


7 W. D. Harkins and A. Feldman, J. Am. Chem. Soc. 44, 
2665 (1922); W. D. Harkins, Sixth Colloid Symposium 
Monograph (1928), p. 31. 

8 Dupré, Theorie Mécanique de la Chaleur, p. 368. 

9W. Hardy, Proc. Roy. Soc. A88, 313 (1913). 
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and T only, so 
AF=vyAo, (10) 


fHy. (11) 


The condition that a liquid } shall by itself 
spread on a liquid a, at constant temperature, is 
thus that the summation of the terms (0F/dc)»,7r 
must be negative, that is less than zero. Thus, if 5 
spreads on a, a certain area of the surface of a 
disappears and certain areas of the surface of b 
and of the interface ab, appear. The increase of 
free energy is 


OF OF 
dF= (—) io+(—) dear 
00/ »,T O0ar/ »,T 


OF 
+(—) dog. (12) 
00 p, T 


For a film of appreciable thickness the effects of 
gravitation must also be considered; but they 
may be neglected here. 

The condition for spreading to occur starting 
from any given condition, is 


dF <0, 
and for non-spreading 


dF>0. 


(13) 


(14) 
If it is assumed that 


doy=dear= —doa, (15) 
which is practically valid for a very thin layer, 
then 


(dF/do) >», T= Yo Yab— Ya- 


(16) 


Let —dF/do be designated as the final spreading 
coefficient (Sy-/4) then 


(17) 


(p, T=constant, the surfaces are saturated as is 
designated by the primes). 


Soja’ = Vai — (Ye + Ya'0’) 


It will be shown later that 


Soja’ <0. (18) 


That is, the final spreading coefficient is always 
negative. 

Now if a drop of the pure liquid b is put on the 
surface of the pure liquid a, and the surfaces of 
these liquids are considered to be still uncon- 
taminated with each other, then the initial 
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coefficient of spreading Ss is given’ by the 
equation 


— (dF /do) p, r= Svja=Ya- (yor Ya'b’). (19) 


If 6 is saturated with a, but the surface of a is 
clean, the initial coefficient may be designated as 
the semi-initial spreading coefficient (S»-/a), 
which, on account of the small effect of water on 
the surface tension of an organic liquid, is not 
very different from Sia when these are the 
liquids involved. 

Now 


Sv ja= Ya—(¥o +Ya'r)- (20) 


The work (Wa) required at constant T and p 
to pull the interface ab. apart to give clean 
surfaces of a and 3, is, for unit area of interface, 
given by: 

WA=Yat7o— Ya'b’ (21) 
and the work of cohesion (Wc) for the liquid 6 is’ 


Wer= 270 (22) 


Wa—Weo=Ya- (yo a's"). 


From (19) and (23) the initial spreading of the 
liquid 6 over the clean surface of a is given by 


—(0F/d) p,r=Svja=Wa—Wes=0. (24) 


(23) 


Thus spreading to form a duplex film occurs if 
Soja is greater than zero, while } remains as a lens 
if it is less than zero. 


5. TRANSFORMATION OF A DUPLEX FILM TO 
Form A Non-DuPpLex THIN FILM 
(PRESUMABLY A MONOLAYER) 


If a non-duplex thin film or a monolayer of 0 is 
present on the surface of the liquid a, the pressure 
(x) of the film is defined as: 


T=Ya—-Vs (25) 


where 72 is the surface tension of the pure liquid a 
and 7; is its surface tension when covered with 
a film of 6. While there is nothing in the thermo- 
dynamics which indicates that such a thin non- 
duplex film is always a monolayer, certain re- 
lations, presented later, seem to indicate that 
this is true. 

What is now desired is the relation between the 
film pressure and the spreading coefficient. 
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TABLE I. Equilibrium film pressures (x), and initial (Svia), 
semi-inttial (Swja), and final (Spa) spreading 
coefficients of oils on water at 20°C. 

(H.0 =72.75.) 








+ OF | y OF 
Dry | Wer 
Ye Om | Ow ; Sp/a 


Methylene iodide |72.20 50.68 |50.52 —23.80 
Carbon disulfide {70.49 32.35 |31.81 ' — 8.23 
Benzene 62.36 28.88 |28.82 |35.0: 8.84 
n-Heptyl alcohol |28.53 26.97 |26.48 37.83 
Isoamyl aleohol {25.92 23.73 |23.56 | & 44.02 


So+/ar 


—24.19 
— 9.95 
— 1.49 
— 5.90 
— 2.64 
































In the treatment of the transformation of a 
duplex film into a monolayer it is essential to 
consider the semi-initial spreading coefficient 
(Soja; Eq. (20)). 

Equation (17) may be written 


(17’) 


where ¥, is the equilibrium value of y; and Eq. 
(20)—Eq. (17’) gives 


Soja’ = Ve— (Yo #+Ya'0'), 


(26) 


Thus z,., the pressure of the monolayer in 
equilibrium with a lens of 6, is equal to the semi- 
initial minus the final spreading coefficient. 

With an oil on water the final spreading coeffi- 
cient is always negative, so the film pressure (m,) is 
always larger than the semt-initial coefficient (Sx ja), 
which has nearly the value of the initial coeffi- 
cient (Sp/a). 


Sosa = Soja’ = Te. 


6. DISCUSSION OF THE THEORIES OF SPREADING 


As a basis for a discussion of the validity of the 
three theories of spreading, the values of the film 
pressure z,, and of the semi-initial and final 
spreading coefficients of five oils on water are 
presented in Table I. They are calculated from 
surface and interfacial tension values obtained 
earlier in this laboratory’ by a special modifi- 
cation of the drop weight method. The determi- 
nations of y,., the surface tension of the water 
covered by the equilibrium monolayer, were 
made on a water surface in the presence of the 
saturated vapor of the oil. The most likely error 
in such a method is that due to undersaturation 
of the film, which would give too high a value of 
Ye and too low a value of z,, which for benzene 
=10.8 dyne cm as determined by the ring 


 W. D. Harkins and B. Ginsberg, Colloid Symposium 
Monograph, Vol. 6 (1928), p. 20. 
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method by L. E. Copeland. In a preliminary 
experiment at room temperature (26.5°C) H. F. 
Jordan" by the use of the ring method, found z, 
for CS. to be 2.6 dyne cm. 

The five liquids are chosen in such a way that: 

A. Three of the liquids have positive initial 
and semi-initial spreading coefficients, and thus 
should exhibit duplex film spreading, as is shown 
to be true by direct experiments on spreading. 
Of the three liquids benzene is a nonpolar 
hydrocarbon, and two are polar-nonpolar alco- 
hols. Isopentane, a paraffin, has an initial spread- 
ing coefficient of 9.4, so its spreading is much the 
same as with benzene, 

B. Two of the liquids, methylene iodide and 
carbon disulfide have negative initial spreading 
coefficients, and so exhibit only non-duplex or 
monolayer spreading. Methylene iodide is chosen 
as a highly polar liquid which does not spread as a 
duplex film on water. 

Carbon disulfide exhibits very interesting re- 
lations, since it does not spread to a duplex film, 
but has a high vapor pressure at 20°. The 
molecules are adsorbed from the vapor upon the 
surface of the water. If this adsorbed film were to 
exhibit a higher equilibrium film pressure than 
that from the lens of the liquid, the material 
would at constant temperature and pressure 
continually evaporate to form a three-dimen- 
sional vapor, which would condense into a two- 
dimensional gaseous film on the water. This film 
would condense into the one-dimensional edge of 
the lens, and the carbon disulfide would complete 
its cycle by returning to the lens. Thus the system 
would exhibit perpetual motion of the second 
kind. Actually a gaseous film spreads rapidly 
from the edge of the lens until the equilibrium 
film pressure of 2.3 dyne cm™ is attained. With 
methylene iodide, which has a very large negative 
value for its initial spreading coefficient as a 
duplex film, the positive film pressure for the 
gaseous film is smaller (0.55 dyne cm~"). 


7. THE INSTABILITY OF DUPLEX FILMS AND 
THEIR TRANSFORMATION INTO MONOLAYERS 


Up to the present time no evidence has been 
presented to show that any stable truly duplex 
(uniform polymolecular) film has ever been 


1H. F. Jordan, Thesis, University of Chicago, 1927. 
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prepared. Such films have always changed into 
a more stable system which consists of lenses and 
a monomolecular film. The closest approach to 
a stable film of this type ever seen by the writer, 
was obtained by spreading an extremely old, and 
therefore partly oxidized, lubricating oil on 
water. For a period of about half an hour, white 
light reflected at about 15° from the plane of the 
surface, gave a uniform color, dependent on the 
thickness, over the whole area which was about 
4X10* sq. cm. Finally, however, the film de- 
veloped a large number of striations and islands 
of varying color. 

An interesting experiment on the spreading of 
any oil which has a positive initial spreading 
coefficient, may be made in a very simple way. 

A large crystallizing dish is made overfull with 
very pure water, and the surface is cleaned 
thoroughly by any one of several well-known 
methods. If the form of the surface of the water is 
correct for the purpose, a layer of the oil may be 
spread so that it covers the surface of the water 
near the edge of the dish, leaving a practically 
circular hole of any desired size, smaller than the 
area of the surface. The instability of a uniform 
layer of benzene of a thickness of 1 mm, or 
somewhat more, may be illustrated by a modifi- 
cation of this experiment. 

The benzene, initially present as a layer several 
mm thick, may be allowed to evaporate. At some 
certain thickness, usually over 1 mm, a hole 
appears in the center of the layer. The non- 
duplex film (monolayer) in this hole expands, 
causing a thickening of the layer, which is now a 
large lens. If the surface of the water is below the 
top of the crystallizing dish, the shape of the 
meniscus is such as to cause the hole, or holes, to 
appear close to the edge of the surface of the 
water, instead of in the center. Suppose the hole 
to appear in the center of a duplex film, of large 
area, instead of in a thick layer as described 
above. The duplex film. may be considered as 
exerting a film pressure (4a= S»/a= 8.9 dyne cm—, 
1™m=ca. 10.4 dyne cm-"). 

On account of its higher pressure the monolayer 
enlarges its area at the expense of the area of the 
duplex film, which changes into a lens. As the 
lens becomes thicker it gains gravitational energy, 
which causes a greater pressure against the 
monolayer. As a result of this the pressure (am) 
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TABLE II. Some lower hydrocarbons which spread on water to 
form duplex films. 
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Paraffins 
Isopentane 
n-Hexane 


n-Octane 
Di-isoamyl 


Aromatics 
8.40 
8.89 
9.23 
9.47 
6.80 


Benzene 


Liquid tt? Shia 
o-Xylene 20° 
m-Xylene 20° 
Mesitylene 20° 
p-Cymene 20° 


6.80 
6.14 
5.54 


Toluene 10.05 


of the monolayer is slightly increased until it 
becomes z,. The theory is treated later. H. F. 
Jordan" finds that a lens of benzene 24.9 cm in 
diameter (v=ca. 14.5 cm*) is 0.257 cm thick, 
while one 16.9 cm in diameter (v=ca. 6.35 cm*) is 
0.212 mm thick. Thus a large sheet of a benzene 
phase on water must be several mm thick before 
it becomes stable, and without the effect of 
gravitation it is not stable at greater thicknesses. 
Presumably a sheet 2 mm thick will not rupture 
unless there is considerable vibration. If a 
circular ring of circular wire is pulled from the 
surface of this hole, the pull in the case of benzene 
(Fig. 1) at 20° is found to correspond to a film 
pressure of 10.6 dyne cm~', or a free surface 
energy of 62.15 erg cm~*. A very large lens or 
layer of benzene on the water has a free energy of 
28.8 erg cm~ in its surface, and of 35.05 erg cm~ 
in the interface parallel to the surface, or a total 
of 63.85 erg cm~*. If this lens could be made so 
thin as to give a duplex film, then, since 
63.85 —62.15=1.7, there would be a decrease of 
free energy of 1.7 ergs for every sq. cm of area 
lost by the duplex film and gained by the 
monolayer. But this is the numerical value of the 
final spreading coefficient. Thus if it were not for 
the work done against gravitation in the piling 
up of a part of the material into lenses, the 
monolayer would grow continually at the expense 
of the duplex film. Actually, if the surface of the 
water is completely covered by a duplex film, this 
transforms itself into a monolayer while the 
remainder of the material of the duplex film 
changes to a layer, lens, or lenses. As the lenses 
thicken the rate of decrease of area of their upper 
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and lower surfaces becomes less than the rate of 
increase in area of the monolayer. This reduces 
the drop in total free surface energy until the rate 
of decrease becomes equal to the rate of increase 
of the gravitational work done, which, obviously, 
gives an equilibrium state. 

Thus if the final coefficient of spreading is 
negative for all oils, as it seems to be, the 
monolayer is always more stable than the duplex 
film. This explains why uniform polymolecular 
oil films on water have not been found to exist for 
any considerable time. 

Whether the velocity of spreading of the 
monolayer is greater or less than that of the 
duplex film is not a question of the free energy 
decrease, but of the rate at which each type of 
film spreads. In an experiment with oxidized 
Nujol it was found that the oil spreads first as a 
duplex film, and that later this breaks up into 
lenses and a monolayer. One thin duplex film of 
this oil which was observed had a life of 90 
seconds, after which it exhibited segregation. 
This life seemed to decrease as the duplex film 
was made thicker. 


8. ADHESION BETWEEN BENZENE AND A 
MONOLAYER OF BENZENE ON WATER 


It is important to learn how far the attraction 
due to molecules of water penetrates into benzene 
at an interface between water and benzene. The 
work of adhesion between the two pure liquids at 
20° is 66.60 erg cm~’, or if the benzene is satu- 
rated with water, 66.55 erg cm~*. However, if the 
water is already covered by the equilibrium 
monolayer, the work of adhesion between this 
and the benzene, is only 56.20 erg cm. The fact of 
interest is that this work of adhesion is 1.7 
erg cm~* Jess than the work of cohesion in the 
benzene. Thus the presence of the water lowers 
the increase in free energy when benzene is pulled 
apart between the first and second layers. In the 
case of heptyl alcohol it requires 6 erg cm~ less 
work to pull this from its monolayer on water 
than to pull heptyl alcohol apart. .This is pre- 
sumably due to the orientation of the heptyl 
molecules in the film so that the paraffin chains in 
the first monolayer, are oriented toward the 
alcohol. Thus the intermolecular forces between 
the hydrocarbon chains and the alcohol, are 
weaker on the average than those in the alcohol. 
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Unfortunately the increase in molecular po- 
tential energy (€4=increase in internal energy) in 
the separation of a liquid from its first monolayer 
on water is not known, since, although the free 
energy has been determined, the latent heat of 
the process has not been obtained. 


9, THE SPREADING OF HYDROCARBONS ON WATER 
AS DUPLEX FILMS 


It has been shown that thermodynamics pre- 
dicts that liquids with a positive initial spreading 
coefficient may spread on clean water to form 
duplex films, and that these transform themselves 
into lenses and a monolayer. In Table II it is 
shown that hydrocarbons, either paraffin or 
aromatic compounds, exhibit positive initial 
spreading coefficients, provided their molecular 
weight is not too high, and that these liquids 
actually spread is shown if they are tested as 
described in the last section. 

This table shows definitely that hydrocarbons 
of not too high molecular weight spread on 
water, so Langmuir’s theory (3) that polar 
groups are necessary for spreading, is not ade- 
quate, since it often leads to conclusions which 
are not in accord with the facts of observation. 

The idea that hydrocarbons do not spread on 
water has part of its origin in the fact that, prior 
to the publication of this paper, there seem to 
have been no values of their film pressures in the 
literature. This is partly because the ordinary 
horizontal type of film balance cannot be used to 
determine this pressure, since the volatility of the 
compound causes it to be adsorbed on the water 
behind the float. 

The explanation of the fact that film pressures 
have always been determined by the use of oils 
whose molecules contain at least one polar group, 
is that this group anchors the molecule to the 
water and lowers the volatility. This not only 
allows the ordinary film balance to be used, but 
what is more important, it makes it possible to 
know how many molecules are present in the 
surface; which would be true of the higher 
hydrocarbons also, but they exhibit negative 
initial spreading coefficients and so will not 
spread. Even tetradecane, a non-spreading oil, is 
bound so lightly to the water and in the film, 
that it escapes too rapidly, as from mixed films, 
to give good measurements. 
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10. LowER HyDROCARBONS FORM MONOLAYERS 
ON WATER 


The geometrical relations of thin non-duplex 
films of oleic acid and other polar-non-polar oils, 
indicate that they are monomolecular, when 
spread from a liquid or a solid three-dimensional 
phase. It is assumed by the writer that the thin 
non-duplex film formed by transformation from 
the duplex film of a lower hydrocarbon, is also 
monomolecular. 

According to Table I the free energy level for 
the monolayer of n-heptféson water is 28.6, 
while that of the duplex ‘film is 34.5 erg cm~?. 
Thus the duplex film will not spread (Sy /a 
= —5.9) over the non-duplex film, which is a 
monolayer. The relations for benzene differ in 
magnitude only. 

The values show that the work required to pull 
n-heptyl alcohol from its monolayer on water is 
5.9 erg cm~? less than that required to pull the 
alcohol from itself. This difference is evidently 
due to the orientation of the alcohol in the first 
monolayer. As a result of this the alcohol has to 
be pulled away from only the hydrocarbon chains 
of the first monolayer, so the work involved is 
less than when the alcohol is pulled apart, since 
some polar groups of the alcohol must be 
separated from each other. While the orientation 
of the benzene does not involve so large an 
amount of energy, nevertheless the most active 
part of the benzene molecule is turned toward the 
water, with its more powerful field, that is the 
least active part of the benzene molecules in the 
first layer is turned toward the second layer. 

As a result of this the smallest molecular 
interaction in the water-benzene system is that 
between the first and the second monolayers of 
benzene. 

From our knowledge of the range of inter- 
molecular forces it is evident that above the 
third benzene layer there can be only a negligible 
attraction of the water at the interface. If this is 
true there is no source of energy which can keep 
the material of the benzene lens, which may have 
a thickness of as much as 2.7 mm. from dowing 
over a film of four molecules in thickne: 3 by the 
action of gravitation. However from the dis- 
cussion in the preceding paragraph it is evident 
that the greatest difference in the energy of 
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interaction in the benzene phase is not between 
the layers 3 and 4, or 2 and 3, but between 1 and 
2. Thus it seems practically certain that the 
stable hydrocarbon film is monolayer. 

It may be assumed that if a monolayer of 
benzene, in equilibrium with a lens, is present on 
the surface of water, a second layer may begin to 
deposit upon it. A molecule in the second layer 
will not be bound as tightly as if in the first layer. 
Therefore (t2) evap < (t1) evap, Where tg and ¢, are the 
life periods if molecules can only escape from each 
layer by evaporation. However, a molecule in the 
second layer can also escape into the first layer, 
where it is held more firmly. This makes 
(te) <(te)evap, Where te is the mean life with re- 
spect to both methods of escape. A molecule of 
the second layer which escapes into the first layer 
replaces a molecule in this layer, and causes a 
molecule in the first layer to move into the lens. 
This last step is accompanied by an increase in 
free energy, but this increase is less than the 
decrease of free energy which is involved in the 
movement of a molecule from the second into the 
first layer. Thus a double film of benzene is not 
stable and only the monolayer remains. 

The “hole” in a hydrocarbon lens disappears 
in a slightly supersaturated vapor of the hydro- 
carbon, not by a deposition of a second and 
a third layer on the monolayer, but by an 
increase in the amount of the hydrocarbon liquid 
phase, and a consequent closing in of the walls 
of the hole. 


11. CHANGES OF FREE ENERGY INVOLVED IN 
THE SPREADING OF A LENS INTO A DUPLEX 
FILM AND A MONOLAYER 


If a drop of an oil which spreads is put on the 
surface of water it forms a transient lens. If the 
area of the water surface replaced by the lens is 
o,, the area of the surface of the oil by o, and 
that of the interface by ow, then the decrease in 
free surface energy for this area a, when the lens 
spreads completely to form a duplex film is* 


—AF 5, r= Yo(Go— 1) FY wo(Fwo — 61) (27) 


* If the lens, which has an area o for its oil surface and 
¢wo for its interface with the water spreads to form a 
duplex film of total area o, the increase of free energy may 
be listed as follows: 


Oil: Yo(a— a0) 
Water: — wo 
Interface: Ywo(a—Gwo). 
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or for unit area 
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a (Yo Ywo)- (28) 


C1 Cl 
There is also a decrease of gravitational potential 
(G), which may be expressed for the whole 
conversion into a duplex film, as 


G=g-f(m, Yay Yb’, Ya'b’s Pay Pb’ Pv), (29) 


The oil of the lens, other than that which 
remains on the area a;, may also be considered as 
spreading to form a duplex film. The decrease in 
the free energy per unit area of the surface of the 
subphase covered in this process is, for the 
spreading of b’ on a 


—AF/(ca— 01) =Ya— (Yor + Y0'v') = — Soja, (30) 


where oq is the area of the duplex film. Thus the 
lowering of the free energy per unit area is given 
by the semi-initial spreading coefficient. 

If a duplex film is formed over the whole of a 
restricted area, the next process which occurs is 
the formation of “holes” in the duplex film. As 
has been shown earlier, the surface of the water in 
these holes is covered by a monolayer of the oil. 
In the formation of the monolayer the free energy 
is again decreased, as follows: 


—AF/om= — Svja’. (31) 


If the total available area is sufficiently small 
the duplex film of area oq changes completely 
into a monolayer of area o,, and a lens of area o/’. 
This area is not that (¢,) of the initial lens, but 


(32) 


In an experiment of limited duration either one 
lens or many lenses may be formed. It is obvious 
that, in the process in which a duplex film is 
changed into a monolayer and a lens, there is an 
increase of gravitational potential (Eq. (29)). 
The free energy of the surface of the water 
covered by the monolayer is per unit area 


(F/dam) tlt ee (Soja — Soja’) =YVw— Te, (33) 


where x, is the film pressure of the equilibrium 
monolayer. 


Thus 


onto =Ca. 


—AF= 06 — Yo(G— 0) —Ywol(t—Fwo). (27+30) 


Unless the lens is extremely small ¢, and ow. are negligible 
in comparison with the total area a, so 


—AF/a=~yw— (Yor + Yw'o’) = — Soriw (30’) 
which is Eq. (30) except for the minute difference in area. 
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* (1941). 
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If the area available is just that occupied by all 
of the oil in a monolayer at the equilibrium 
pressure, no lenses are left, and no work is done 
against gravitation. 

If a very much larger area than this is available 
the monolayer spreads over it, with a con- 
comitant decrease of film pressure and increase of 
vs. While the free energy per unit area of the film- 
covered surface is thus increased in a spontaneous 
isothermal, isopiestic process, the free energy (F) 
of the area as a whole is decreased. That is, while 
the free energy of that part of the surface of the 
water covered by the film is increasing, the free 
energy of the clean water is decreasing even 
more rapidly, as the film spreads over it. 

If a monolayer, with the free surface energy y. 
is not in contact with a duplex film or lenses, the 
rate at which ¥, increases with increase of area at 
constant temperature depends upon the state of 
the monolayer,'? which may be, for example, the 
liquid expanded (L;) state or the gaseous state. 
If in the former of these two states the monolayer 
at a certain stage of expansion undergoes a 
change of the first order in which there is a 
transition to a two-dimensional gas. Throughout 
the extensive area change through which this 
transition occurs the pressure remains constant as 
the area increases. It is obvious that the pressure 
(x;) of the gaseous film approaches zero, and the 
value of y; approaches that of water as the 
monolayer expands to very large areas. 

Thus the history of the expansion of a lens and 
a duplex film involves two steps in each of which 
the free surface energy decreases sharply (Fig. 1). 
The first, and usually the greatest, decrease is due 
to the spreading of the duplex film over the water 
surface, and the second, to the conversion of the 
duplex film into a monolayer and a lens, if the 
area is limited, or if sufficient area for its forma- 
tion is available, into a monolayer alone. Further 
spreading of the equilibrium monolayer produces 
an increase in the free energy of that part of the 
surface covered by the oil. 


12. THE SPREADING OF MIXED OILS 


A non-polar oil, which has a negative spreading 
coefficient (Sb/q or Sy/2) may be converted into an 
oil which spreads to form. a larger lens or to a 


#2 W. D. Harkins and E. Boyd, J. Phys. Chem. 45, 20 
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duplex film by the addition of a polar-non-polar 
substance, such as stearic acid. For tetradecane 
on water at 25° Langmuir" finds S»/_ to be — 7.6, 
while the value for a sample of Squibb’s petro- 
latum at 20° was —11.4. In the relation Stja= a 
—(yotYav) the values are: 


tetradecane: —7.6=72.0—(26.9+51.3) and 
petrolatum : —11.4=72.75—(31.1+53.0). 


If myristic acid or another polar-nonpolar 
material is present in the oil, it is partially 
adsorbed by the interface. This lowers Ya'», which 
gives a temporary value of the spreading coeffi- 
cient which is closer to zero, or may be positive. 
This causes the lens to expand with a reduction 
of the contact angle. As the lens expands the 
amount of acid adsorbed is increased on account 
of the increase of interfacial area, and this lowers 
the concentration of acid in the oil. This reduces 
the adsorption of the acid and thus increases the 
interfacial tension (decreases ya»). As a result the 
lens would attain a definite larger area with a 
value of Sy-/4 closer to zero, if no myristic acid 
were to spread from the lens to give a monolayer. 
If the surface of the water has a limited area any 
acid which escapes gives a film pressure, which 
depends on the concentration of the acid in the 
monolayer. This decreases ya, which, in turn, 
decreases the value of the spreading coefficient, 
and causes the lens to become smaller. 

In the non-equilibrium processes associated 
with changes in the area of the lens, the rate at 
which the monolayer spreads from the lens, 
through the linear region at its boundary, is a 
very important factor, which is not always kept 
in mind. This rate is often very small. 

When equilibrium is finally attained the rela- 
tive concentrations of the polar-non-polar sub- 
stance in the interface (ab) and the surface of the 
water (a) should be constant, or C,/Ca,=k. In the 
few cases for which values of k have been 
determined it has been found to be much larger 
than unity. Thus Harkins and King" found that 
many more molecules of butyric acid are adsorbed 
at the surface of water than at the interface 
water-benzene at a given concentration in the 


13], Langmuir, J. Chem. Phys. 1, 763 (1933). 

144 W. D. Harkins and H. H. King, J. Am. Chem. Soc. 41, 
982 (1919); F. Askew and J. Danielli, Trans. Faraday Soc. 
36, 785 (1940). 
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water, provided the film is not so tightly packed 
as to have its concentration determined by 
geometrical considerations. The value of k for 
trimethyl! cellulose or a-amino-palmitic acid at 
the water-brombenzene interface was also found 
to be much larger than unity (Askew and 
Danielli). In all instances in which k>1, it is 
evident that, afterequilibrium is attained between 
an oil lens and water, the concentration of the 
solute should be greater on the surface outside 
the lens than in the interface under it. 

If an acid, such as myristic acid, is present in 
the lens its spreading is highly dependent upon 
the pH of the aqueous subphase, and, if this is 
alkaline, upon the extent to which bivalent, 
trivalent, and quadrivalent positive ions are 
present. 

It was found by Mercer that a monolayer 
escapes from the edge of a non-polar oil in which 
stearic acid is dissolved, even when the subphase 
is alkaline. Langmuir" finds that if more than 
0.006 g stearic acid per g of Squibb’s petrolatum 
is present the oil spreads on water to a film so 
thin as to show interference colors, and that the 
molecular area of the stearic acid at the interface 
is 64A*. With more than 0.010 g of stearic acid 
the duplex film is too thin to exhibit interference 
colors. 

When duplex films are formed in this way the 
free energy diagram is of the type exhibited by 
benzene in Fig. 1. The general relations are 
unchanged, but the free energy level of the 
duplex film is now dependent upon the concen- 
tration of the polar-non-polar solute as well as 
upon the nature of the oil. The concentration of 
the solute depends upon the area to which the 
oil has spread. If the solution is basic, calcium 
ions reduce the molecular area of the acid at the 
interface to about a third of 64A®. Preliminary 
measurements in this laboratory by B. Adinoff 
show that Cat+ ions reduce the area on the 
surface of water from a value of 24.4 to 21.15A at 
zero pressure (extrapolated). 

The duplex mixed film is always unstable, and 
breaks up into the final state of a monolayer and 
a lens. A drop of oil of the composition of this 
equilibrium lens will not spread over the surface 
of the monolayer in equilibrium with the lens. 


15 E, H. Mercer, Proc. Phys. Soc. London 51, 561 (1939). 
16], Langmuir, J. Frank. Inst. 218, 165 (1934). 
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Thus the general relations for spreading as duplex 
films are the same as when a pure oil is involved, 
except that the free energies and the related 
spreading coefficients are made variable by the 
changing composition of the oil of the lens. 


13. NON-EQUILIBRIUM INTERFERENCE PATTERNS 
EXHIBITED BY DUPLEX FILMS, AND OTHER 
Non-EQUILIBRIUM PHENOMENA 


The phenomena associated with the non-equilibrium 
state of the duplex film and of its conversion into a mono- 
layer are very complex, and remarkably beautiful inter- 
ference patterns of extremely varied shapes and colors may 
be obtained in ordinary white light. On account of the fact 
that the light reflected by the surfaces of the layer is 
partially elliptically polarized, it is instructive to view the 
film through a Wollaston or a Nicol prism, when the elec- 
tric vector is parallel and when it is perpendicular to the 
plane of the surface. 

If a fine stream of benzene is allowed to flow slowly 
upon the clean surface of water in an open vessel, the 
benzene may be seen to advance over the surface as a 
duplex film, at the edge of which there appear narrow 
colored interference fringes, which exhibit as many as 
four or more orders of the spectrum. For a short distance 
beyond this the film is still duplex, but still beyond this 
spreads as a non-duplex film, presumably a monolayer. 
Islands of duplex film, which exhibit rapid, lifelike pul- 
sating movements are formed if evaporation occurs at the 
proper rate. In some instances the motions resemble 
those of amoebae. 

A duplex film too thin to give interference colors may be 
obtained by spreading one or two drops of an oxidized 
hydrocarbon oil (as oxidized Nujol) over a large water 
surface. On compression by a barrier of glass or metal, 
colors appear, and these are uniform over the whole area 
of the surface (about 4103 cm? in the experiments). The 
colors appear in the sequence pale yellow, yellow, gold, 
orange, brown, red, purple, blue, pale blue, green, pale 
yellow, etc. Finally, in the highest orders of interference, 
red and green are the predominant colors. 

On standing the duplex film of such a material may 
segregate into duplex films of different thicknesses. For 
example, a delicate lace pattern of interlocked fine curved 
lines of red on a field of green, has been observed fre- 
quently. Loose ends of the red lines are often coiled into 
spiral. The pattern seems to depend upon the chemical 
nature of the oil, its complexity of composition, and its 
viscosity. 

Often there is a segregation into flat duplex islands of 
uniform color, except for small colorless holes or lakes of 
monomolecular thickness. Between the islands there are 
canals of a colorless monolayer, studded with colorless 
islets, which are actually very small lenses. 

The patterns obtained with different oils are so manifold 
that it is impossible to describe them. 

If such a nonhomogeneous colored duplex film is com- 
pressed further, it gives apparently flat colorless duplex 


islands, with lakes, canals, etc., just as when the film is 
colored. Finally, when equilibrium is attained, a monolayer 
and a lens remain. 

Hardy" calls attention to an interesting phenomenon 
which he has observed with small lenses of carbon disulfide. 
This is exhibited by other highly volatile liquids whose 
lenses are imbedded in monolayers. Two small lenses, if 
they come within 1 or 2 cm, of each other, act as if 
violently attracted and move rapidly toward each other 
until close together, when they are just as violently re- 
pelled. However, if the velocity of approach exceeds a 
critical value, the two lenses fuse into one. 

Monolayers may in general be compressed to higher 
values of x, that is lower values of ys, than corresponds to 
equilibrium with a lens or crystal of the material at the 
given temperature. For example, at 20°C an ordinary 
crystal of stearic acid gives an equilibrium film pressure 
of only 3.1 dyne cm™, but by rapid compression it is 
possible to obtain the pressure-area (x-o) relation without 
much error up to film pressures as high as 50 dyne cm™, 
that is, toa free surface energy level as low as 22.8 erg cm~. 
In fact, a considerable fraction of the data in the literature, 
which have been used to give r-o diagrams, have been 
obtained in the regions of instability. 


14. THE Non-SPREADING OF WATER ON 
ORGANIC LIQUIDS 


From the relation for the initial spreading 
coefficient 


Svja= Wa— Wee=va-— (Yo Ya'b’) 


it is found that water will not spread on any 
organic liquid for which the relations are known, 
and presumably not on any organic liquid. Here a 
is an oil and 0 is water. For isopentane S./.= 13.7 
— (72.8+49.6) = — 108.7, an extremely large neg- 
ative value. For carbon disulfide S,,.=31.4 
— (72.8+-48.6) = —90. For water on primary iso- 
amyl alcohol, Ss;¢=23.7 —(72.8+5.0) = —54.1. 
However the greater mutual solubility of a pair 
of liquids such as these, reduces greatly the value 
of the final spreading coefficient. Thus for this 
alcohol Sy ja =23.6—(25.9+5.0) = —7.3, which, 
while still negative, is very small for water on an 
oil. However, it is obvious that the smallness of 
this value is due to the great reduction of the free 
surface energy of the water which is caused by the 
oil. 


15. THE SPREADING OF LIQUIDS AND SOLIDS ON 
MERcurRY (FIG. 2) 


The free surface energy of mercury in erg cm~? 
is: 480 at 0°C; 476 at 20°C, and 467 at 60°C. The 


17W. Hardy, Proc. Roy. Soc. A88, 560 (1913). 
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Fic. 2. Water and all 
organic liquids have posi- 
tive initial spreading coef- 
ficients on mercury, and 
thus may spread to form 
duplex films. The mono- 
layers lie at a lower free 
energy level, so the duplex 
film, if formed, transforms 
into a monolayer and a 
lens. The free energy 
levels of the duplex films 
are known and are prop- 
erly placed in the dia- 
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total surface energy (h or e) is 540. Thus if a free 
energy diagram, similar to Fig. 1, is constructed 
for the spreading of liquids or solids on mercury, 
the level for the mercury surface at 20°C is very 
high (476). The result of this is that all duplex 
and all monolayer levels for water and for oils lie 
far below the mercury level. There is a large 
decrease in free energy if water or any liquid oil 
spreads as a duplex film on mercury. Thus water 
and all organic liquids spread on mercury. 


nn ae yi 


Monolayer 


gram. Those for the mono- 
layers have not been 
determined, so only the 
sign of the final spreading 
coefficient (.S’) is known. 
While the film pressure x, 
is larger than the initial 
spreading coefficient Sp). 
it is probable that it is not 
very much larger, so the 
diagram gives the general 
location of the free energy 
level for each monolayer. 





Free Energy Diagram 
for the 
Spreading of Liquids on Mercury 


There are no reliable determinations of the 
final spreading coefficient S»/a or of the equilib- 
rium film pressure for any liquid on mercury. 
However the lenses which form on the non- 
duplex film, presumably a monolayer, indicate 
that this film pressure is somewhat greater than 
the initial spreading coefficient, Soa. In ergs per 
cm~? this quantity has the values: water, 32: 
acetone, 60; hexane, 79; propyl alcohol, 85: 
carbon tetrachloride, 93; benzene, 90; chloro- 
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form, 97; octyl alcohol, 102; ethylene bromide, 
116; ethyl iodide, 135; and methyl iodide, 137. 

Many organic solids are found to spread on the 
clean surface of mercury. 


16. THE SPREADING OF LIQUIDS ON THE SURFACE 
OF CRYSTALLINE SOLIDS 


Insofar as the periodic surface field of a solid 
may be considered as a uniform field, the free 
energy relations involved in the spreading of a 
liquid over the surface of a solid are the same as 
when the subphase is a liquid, except that with 
inorganic solids of high melting point W, is very 
much higher. The phenomena, however, are very 
different. In the atmosphere practically all 
inorganic solids are covered with a film of water, 
and if this is the case the spreading of another 
liquid or of more water is often prevented. Both 
the rigidity and the periodic electrical field of the 
surface of the solid, have a great influence on the 
rate of spreading. 

The spreading of a pure liquid over the clean 
surface of a solid does not seem to occur, as it 
may on a liquid, by the pulling out of a duplex 
film from the lens, but by the spreading of a 
monolayer. The rate of formation of the mono- 
layer by the process of spreading may, however, 
be less than its rate of formation from the vapor 
of the liquid. It is obvious that this condensation 
may form a polymolecular film on the surface of 
the solid, and, on account of the lack of mobility, 
neither the surface of the solid nor the first 
monolayer will move to make a place for the 
molecules of a second layer. Hardy'*® gives an 
interesting example of spreading on glass. A drop 
of anhydrous acetic acid remains drawn up as a 
lens on a clean dry glass plate, but if the atmos- 
phere is admitted for even a moment, the slight 
amount of moisture causes the drop to flash over 
the whole surface. 

It is now well known that films adsorbed on 
solid surfaces may be mobile, as shown by Volmer 
and collaborators,'® or immobile. It is not the 
purpose of this paper to deal with immobile 
layers. It has been shown recently” that certain 


‘*W. Hardy, Collected Papers, p. 713. 

_'* Volmer and Estermann, Zeits. f. Physik 7, 13 (1921); 

Volmer and Adhikari, Trans. Faraday Soc. 28, 359 (1932). 
** W. D. Harkins and F. M. Fowkes, J. Am. Chem. Soc. 

60, 1511 (1938); F. M. Fowkes and W. D. Harkins, ibid. 62, 

3377 (1940). 
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short chain organic acids, salts, alcohols, and 
amines, form the 
interface between water and a solid, such as 
stibnite, graphite, tale, or paraffin. 

The conclusion which may be drawn from the 
above relations is that the spreading coefficient 
relations presented in this paper are applicable to 
the spreading of mobile films over solids. If 
duplex films do not spread directly, it is obvious 
that the steps in which they are involved should 
be omitted. 

With solids such as BaSO,, TiOs, ZrSiO,y, ZrO>, 
and all other solids of high cohesion, the initial 
spreading coefficient is undoubtedly positive. 
That is, the spreading of the liquid film over the 
surface of any solids of this type gives a decrease 
of free energy. On solid paraffin however, the 
coefficient for water, as is well known, has a large 
negative value. This is true for many other 
organic solids. 


gaseous adsorbed films at 


17. THE PRINCIPLE OF INDEPENDENT SURFACE 
ACTION AND THE SPREADING OF CARBON 
TETRACHLORIDE ON WATER 


The extent to which the principle of inde- 
pendent surface action is correct is exhibited in 


TABLE III. Increase in molecular potential energy when a 
liquid is pulled apart (ec), or from another liquid (e,4) (erg: 
for a bar of liquid 1 cm* in cross section, for the separation of 
the upper from the lower liquid) (t=20°C).* 








n-Octane n-Octyl Alcohol n-Heptylic Acid n-Heptane 





100 


n-Heptane 


105 
n-Heptylic Acid 


100 
n-Octane 


101 
n-Octyl Alcohol 


n-Octyl Alcohol 








n-Octane n-Heptylic Acid n-Heptane 





107 165 151 106 





Water Water Water Water 
Effect of the Length of the Hydrocarbon Chain on the 


Cohesional Energy 


n-Butane n-Butyl alcohol n-Propylalcohol Propane 





92 97 96 68 





n-Butane n-Butyl alcohol mn-Propylalcohol Propane 








* The shortening of the chain a paraffin hydrocarbon produces its 
first very marked effect in lowering the cohesional energy in the decrease 
from 4 to 3 carbon atoms in the hydrocarbon chain. Shortening the 
chain has little effect with the alcohols, which give 97 and 98 erg cm™* 
for ethyl and methyl alcohols, respectively. 
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the early work of the writer and his collaborators?! 
on the free energy relations of surfaces and 
interfaces, but is shown even more perfectly by 
the magnitude of the increase of molecular 
potential energy when liquids are pulled apart,”* 
at a plane inside the body of liquid. The relations 
are illustrated (Table III) by the numbers which 
give the increase in potential energy in ergs when 
a column of liquid 1 cm? in cross section is pulled 
into two parts. 

It is remarkable that the increase of molecular 
potential (internal) energy which accompanies 
the separation of the polar-non-polar liquid 
(alcohol or acid) is only 1 or 5 percent higher 
than when the non-polar saturated liquid 
(n-octane or n-heptane) is involved. Since this 
amount of energy, of the order of 3 ergs per cm of 
interface, is surprisingly small for the energy of 
orientation of the molecules of octyl alcohol or 
heptylic acid, it is evident that the polar groups 
exert no appreciable effect on the energy of 
separation of the liquid when this separation 
occurs between the hydrocarbon chains only, 
that is after the molecules are oriented at the 
plane of separation (Fig. 3). This is not a sur- 
prising result, as the distance between the polar 
groups under these conditions is of the order of 
22A and it is known that for a given relative 
orientation of a pair of dipoles, the interaction 
energy decreases as the cube of the distance 
between them. 

Even though the dipoles in water molecules 
give a strong interaction with other permanent 
dipoles it is evident that only six or seven percent 
more energy is required to separate either heptane 
or octane from water than from itself. This 
indicates that the dipoles induced in the paraffin 
by the water are weak. Obviously the London 
dispersion forces are not the same when water 
constitutes one of the phases at the phase 
boundary, as when an oil is on both sides of an 
imaginary boundary. 

If the acid or alcohol is pulled from water at 
the interface the close proximity of the polar 
groups results in a large increase in the energy of 


21W. D. Harkins, F. E. Brown, and E. C. H. Davies, 
J. Am. Chem. Soc. 39, 354 (1917); W. D. Harkins, E. C. H. 
Davies, and G. L. Clark, ibid. 39, 541 (1917). 
( sank D. Harkins, Zeits. f. physik. Chemie A139, 670 
1928). 
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Fic. 3. Idealized diagram illustrating the separation of a 
bar of liquid alcohol into two parts. The molecules are 
represented as already oriented at the interface, since the 
break is well under way. This diagram indicates that in such 
an ideal case the energy of separation of a moderately long 
chain alcohol should be practically equal to that of a 
hydrocarbon of the same molecular size, except for the 
relatively small energy of orientation of the molecules at 
the interface before the actual separation has occurred. 


separation to 165 ergs for the alcohol and 151 
ergs for the alcohol. 

The principle of independent surface action 
indicates merely that the behavior of surfaces 
and interfaces is in accord with the idea that the 
range of intermolecular forces is very small. 

By the use of this principle Langmuir®* reaches 
the conclusion that ‘‘in general the properties of 
the upper surface of an adsorbed monomolecular 
film are nearly the same as those of the substance 
in bulk.” 

As an example illustrating the correctness of 
this point of view he states : “Carbon tetrachloride 
vapor forms an adsorbed monomolecular film on 
water which tends to become nearly complete 
before the vapor becomes saturated. The prop- 
erties of the upper surface of the film thus 
approximate to that of liquid carbon tetrachloride, 
so that the life of the molecules of CCl, which 
may condense on the film is about the same as 
that of molecules on liquid CCly. Thus the film 
grows in thickness which means that F,=0.”’ 


237. Langmuir, J. Chem. Phys. 1, 775 (1933). 
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LIQUID SPREADING 


In contrast with this, all determinations of the 
initial spreading coefficient of CCl4, —(0F/de) p,r 
=Sra=Fs, give a small positive value, and a 
direct experiment shows that the film which 
spreads on a clean water surface from a lens of 
the liquid ‘‘pushes”’ talc powder in front of it to 
form a circular area free from talc. After the film 
has spread any excess of liquid remains as a lens, 
almost flat above and quite curved below. The 
stability of the lens indicates that the final 
spreading coefficient is negative. 

The fact that the lens remains on the surface in 
a saturated atmosphere of the vapor indicates 
that the mechanism postulated by Langmuir is 
incorrect. Thus the film does not thicken beyond 
a monolayer by condensation of the vapor, and 
the spreading relations of CCl, are not changed, 
by its low value of Sj,/a, from those postulated 
earlier in this paper for liquids which have higher 
positive initial coefficients. Even if the initial 
coefficient were exactly zero this would not allow 
the monolayer to become thicker, since the final 
coefficient (S/a.) would still be negative. As 
pointed out earlier, this causes any molecules, 
which cannot be accommodated in the monolayer, 
to move into the lens. 


18. THe SPREADING OF A Souip As A Non- 
DupLEX FILM 


Acids, alcohols, amines, and other substances 
of the polar-non-polar type spread readily on 
water if the molecular weight is not too large. 
This fact has given rise to the idea that the 
presence of a polar group or groups in the 
molecule always gives rise to spreading, but this 
is incorrect. What the polar group does, in 
general, is to make the adhesion of the oil for the 
water relatively greater, when compared with the 
work of cohesion of the oil itself, than if the 
polar group were absent. 

However, if the nonpolar part of the molecule 
is a sufficiently large fraction of the whole the 
cohesion may become high enough to inhibit 
spreading. In the case of an oil a sufficient in- 
crease in cohesion causes the material to take on 
a solid, rather than a liquid state. The solid state 
does not, of itself, prevent spreading, but does 
eliminate spreading as a duplex film. Thus the 
films formed by the direct spreading of a solid are 
monolayers. 
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A monolayer formed by the direct spreading of 
a solid (e.g., an acid, alcohol, amine, etc.) may 
be in the solid state, or any one of the known 
liquid states. If the equilibrium film pressure is 
too low to give one of these more or less condensed 
states, the solid may, nevertheless, vaporize in 
the surface to give a two-dimensional gaseous film. 

Crystals of myristic acid exhibit non-duplex or 
monolayer spreading on the surface of clean 
water. The type of monolayer which is in 
equilibrium with the crystals is dependent upon 
the temperature. Thus at 12.5°C a condensed 
liquid (Lz) monolayer is formed, at 16° the film 
is in the intermediate, and at 23° in the liquid 
expanded (L;) state. At 30°C octadecyl alcohol 
spreads to give at equilibrium a monolayer which 
is in the solid state, with a film pressure (7,) of 34 
dyne cm. The equilibrium pressure of an 
alcohol is much higher than that of an acid, if the 
hydrocarbon chain is the same in both. 

As the length of the hydrocarbon chain is 
increased, the equilibrium film pressure, at a 
given temperature, decreases. As a result of this, 
acids with sufficiently long chains spread to give 
gaseous monolayers only. This seems to be true 
of arachidic acid, the normal acid with 20 carbon 
atoms per molecule, which, at 20°C, gives an 
equilibrium pressure of only a few tenths of a dyne 
per cm. If the number of carbon atoms is 19 or 
21 the equilibrium pressure is considerably higher. 

It is evident that if the length of the hydro- 
carbon chain of a polar-non-polar or a non-polar 
substance is increased sufficiently, the equilib- 
rium pressure of the gaseous film will be so low 
that it cannot be measured or detected by a film 
balance of any sensitivity now employed. From 
the practical standpoint such a substance can be 
considered as non-spreading. 

Thus non-spreading organic substances are not 
limited to those whose molecules are wholly non- 
polar. If a polar group is presetit,spreading may 
be attained by an increase in size of the non- 
polar part of the molecule. 


19. THE GENERAL EFFECT OF THE TYPE OF 
MOLECULE UPON THE SPREADING OF A 
LIQUID ON WATER 


Table IV shows how the stability of the 
duplex and monomolecular types of film depends 
upon the molecular nature of the liquid which 





G. KANE 


TABLE IV. Effect of molecular constitution upon the formation of duplex films and of monolayers. 








ATTRACTION 
FOR WATER 
MOLECULES OR GROUPS Wa 


ATTRACTION 
FOR LIKE 
MOLECULEs Wo 


STABILITY OF 
MONOLAYER 


DupLex Fim! 





Weak 


60 
Moderately 
strong 


Hydrocarbons of high molecular 
weight 


—CH.I, —CH.Br, —CH.,I 

Hydrocarbons of low molecular weight 

Esters 

—CH,0OH, —COOH, —CN, 
—CONH:2, —CsH,OH, 


— NHCONH,, etc. 
—SO;H, —CsH,SO;H, —SO.H 


Very strong 
100 


Extremely 
strong 


Extremely 
strong 


Soaps (solids) 


Weak 


Strong 


Weak 


Moderately 


Moderately 


None None or gaseous 

65 

None None or gaseous 

100 

Unstable Gaseous 

Stable. Relatively high 
compressibility 

Stable liquid, solid and 
gaseous 


Unstable 


Unstable 


Forms if material Form monolayers 
liquid. 
Unstable 


None 


weak 


Form monolayers 
weak 














‘If the spreading material is a solid, a duplex film does not form, but only a monolayer spreads. See Section 12 for the spreading of a duplex 
film of an oil which, of itself, does not spread in this way, by the addition of a polar-non-polar oil as a second component. 


forms the spreading or non-spreading phase. The 
free energy values listed represent typical 
magnitudes. It is obvious that these vary con- 


siderably with chain length and other factors 
which are not indicated with any completeness in 
the table. 
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The Equation of State of Solid Helium 
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Manhattan College, New York, New York 


(Received June 6, 1941) 


HE idealized crystal is one in which par- 
ticles form a perfectly ordered array and 
execute small oscillations about equilibrium posi- 
tions. This hypothetical picture may be used as 
a starting point to determine the equation of 
state of a solid and to establish a theory of fusion. 
The simplest case to consider is that of frozen 
rare gases. The only forces existing between two 
rare gas atoms are the Heitler-London forces of 
repulsion and the van der Waals forces of attrac- 
tion. The potential energy of the crystal is deter- 
mined by making a summation over the whole 
lattice of the potential energy of interaction be- 
tween pairs. The elastic constants may be deter- 
mined from the interatomic forces and the fre- 
quency of oscillation from the elastic constants. 


The total energy of the crystal is 
U=N/2+3Nhi, 


being made up of the lattice potential energy and 
the zero-point energy. The minimum of total 
energy should give the heat of sublimation at 
absolute zero and the lattice distance between 
atoms. The shift from minimum of potential 
energy to minimum of total energy, which is 
produced by the zero-point energy, causes a 
change in both lattice distance and heat of sub- 
limation. This change is most noticeable with 
the lighter elements. For temperatures above 
absolute zero, the free energy is 


A=U-TS. 





‘EQUATION OF STATE OF SOLID 


TABLE I. Energy ts in ktlojoules, rin angstroms, V in cm, 
p in atmospheres. @=hv/k in degrees Kelvin. Figures in 
column containing minimum pressures are italicized. 








EXPERI- 
MENTAL 
VALUES 


2.95 ‘ 3.57 
10.93 of 17.4 

0.584 

0.793 

0.209 0.050 
62.5 5 32 

1600 : 25 


CALCULATED VALUES 


2.85 2.93 
9.86 10.39 10.66 
0.852 0.682 0.630 
0.880 0.839 0.816 
0.028 O.157 0.185 
91.0 72.8 67.4 
2470 = 1440 830 





2.90 


2020 








From this, the equation of state may be estab- 
lished since 


p=—(0A/dV)r. 


Each isotherm drawn on a pressure-volume dia- 
gram has a minimum. Since any state is unstable 
where pressure increases with volume, the crystal 
must break down at this minimum pressure. The 
temperature of the isotherm whose minimum 
pressure is zero is interpreted as the melting 
point. 

Following this line of procedure, the calcula- 
tions were made for frozen argon and helium by 
Herzfeld and Goeppert Mayer.! For helium, they 
used the formula derived by Slater and Kirk- 
wood ;? 


¢= {770e-*- 6" — 1.49 /r®}10-” erg. 


They considered helium as a cubic face-centered 
crystal and calculated the lattice distance, melt- 
ing point, and Debye temperature. The agree- 
ment with experimental results was poor. For 
argon, they used the form 


g=ABe"—A/r*, 


in which the coefficients A and B were to be 
determined so as to make the minimum of the 
total energy curve coincide with the experimental 
interatomic distance and heat of sublimation at 
absolute zero. Two values were chosen for p, 
one deduced theoretically for two neon atoms 
(p=0.2091A), and the other for alkali halide 


' K. F. Herzfeld and M. Goeppert Mayer, Phys. Rev. 46, 
995 (1934). In Table V the values for the potential energy 
and zero-point energy should be interchanged. 

“9st Slater poe J. G. Kirkwood, Phys. Rev. 37, 682 
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crystals (9p =0.345A). Kane* made the calcula- 
tions for neon, krypton, and xenon and recalcu- 
lated the results for argon. The agreement with 
experimental melting point was fairly good but 
not remarkable. He showed also that the use of 
second-order van der Waals terms produced only 
slight changes. 

In the present paper the equation of state 
is deduced for solid helium. For the potential 
energy of interaction between two helium atoms, 
the values determined by Margenau‘ have been 
used 


y= {770e-* — 560e-5-33 
— 1.39/r'—3.0/r°}10-” erg. 


The crystal has been considered as hexagonal 
close-packed with an ideal axial ratio.’ The sum- 
mations for potential energy and those to be used 
in the frequency formula have been made by 
Kane and Goeppert Mayer.* The formula for 
the frequency has been derived directly from 
the equations of motion of a particle. It is diffi- 
cult to interpret the terms of the formula as 
elastic constants,’ but there is no doubt about 
the validity of the result. The other rare gases 
have cubic lattices and therefore the symmetry 
relations make it quite easy to identify the terms 
as elastic constants. In a previous paper*® the 
author stated that the frequency formula that 
Herzfeld and Goeppert Mayer had used was 
valid only at the minimum of potential energy. 
It now appears that their formula was correct 
and, consequently, that Kane’s was wrong. The 
error arose from using Born’s® formula for fre- 
quency which does not hold for cases of displace- 
ment from equilibrium. Any change resulting 
should be greatest for neon. Therefore, the cal- 
culations were repeated for this element using 
one value of p (p=0.2091). It was found that 
the error was slight and that deviations in the 
pressure-volume curves were negligible. 


3G. Kane, J. Chem. Phys. 7, 603 (1939). 

4H. Margenau, Phys. Rev. 56, 1000 (1939). 

5 W. H. Keesom and K. H. Taconis, Proc. K. Acad. 
Amst. 41, 95 (1938); Comm. Leiden No. 250e. 

6G. Kane and M. Goeppert Mayer, J. Chem. Phys. 8, 
642 (1940). Summations for hexagonal lattice are slightly 
larger than for cubic face-centered. 

7M. A. Biot, Phil. Mag. 27, 468 (1939). 

8M. Born, Dynamtk der Kristallgitter (Leipzig, 1923), 
p. 648. 
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The results for the calculations on helium 
using Margenau’s formula have been tabulated 
in Table I. The agreement with experimental 
values is rather poor in all cases. The difference 
between experimental interatomic distance and 
the positions of minimum total energy may be 
accounted for to some extent by the fact that 
what is measured experimentally is not the equi- 
librium distance but the time-average position. 
Since the vibrations of helium are very large 
and anharmonic it follows that there should be 
considerable difference between the two. There 
seems to be no method available for correlating 
these two positions, since the crystal breaks down 
completely shortly beyond the equilibrium point 


due to the fact that the elastic constants become 
negative. 

The method of Lennard-Jones® with assump- 
tions of order-disorder transitions, may be ap- 
plied to the equation of state of helium. This 
has been done by Nagamiya,'® using Margenau’s 
potential energy formula. His agreement with 
experimental results is remarkably good. How- 
ever, in using Margenau’s results, he changes the 
sign before 560e—5-** from minus to plus. This 
change of sign invalidates his results. 


9 J. E. Lennard-Jones and A. F. Devonshire, Proc. Roy. 
Soc. A169, 317 (1939). 

10 0) Nagamiya, Proc. Phys. Math. Soc. Japan 22, 492 
(1940). 
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The Restricted Rotation Potential and Frequency 
Assignment in Propylene 


D. TELFAIR AND W. H. PIELEMEIER 


Physics Department, The Pennsylvania State College, 
State College, Pennsylvania 


June 5, 1941 


HE assignment of the fundamental vibrational fre- 

quencies given by Pitzer,! and a potential restricting 
rotation of the methyl group in propylene of about 800 
cal./mole are not in disagreement with the heat capacity 
measurements reported by Kistiakowsky and Rice.? How- 
ever, Wilson and Wells’? frequency assignment and a 
potential restricting rotation of about 2000 cal./mole are 
also in accord with the reported values of the heat 
capacities. 

The temperature range of heat capacity measurements 
has recently been extended in this laboratory by supersonic 
methods with results which make possible a choice between 
the two suggested potentials and frequency assignments. 
As shown in Fig. 1 the measured heat capacities at the 
higher temperatures (between 360 and 500°K) are still in 
excellent agreement with the Wilson-Wells assignment and 
a 200 cal./mole potential, but values calculated from 
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Pitzer’s assignment and an 800 cal./mole potential deviate 
considerably more than the limits of accuracy of the 
measurements. 

The entropy of propylene calculated on Pitzer’s basis 
agrees with the measured values if one assumes an entropy 
of mixing due to the inability of the solid to distinguish 
between the two ends of the propylene molecule. However, 
Aston and his co-workers have recently found that cyclo- 
pentane is non-planart by comparison of the third law 
entropies with the spectroscopic values. The five corners, 
although not all lying in the same plane, must resemble 
each other more closely than do the two ends of propylene. 
Aston’s results show no entropy of mixing in the solid 
phase, and, therefore, an ordered arrangement in the 
crystal. Had the crystal not distinguished between the 
CH: groups in cyclopentane, the third law entropies would 
have appeared to correspond to a planar ring in spite of 
its non-planar structure. 

Thus an entropy of mixing cannot be expected in pro- 
pylene. Calculated on this basis, the Wilson-Wells fre- 
quency assignment and a 2000 cal./mole potential agree 
well with the measured entropy, whereas Pitzer’s frequency 
assignment and an 800 cal./mole potential do not. 

A description of the apparatus and justification of the 
methods used in the heat capacity measurements will 
appear later. We are indebted to the Pennsylvania State 
College Cryogenic Laboratory for furnishing the purified 
sample of propylene used for these measurements. 

1K.S. Pitzer, J. Chem. Phys. 5, 473 (1937); Chem. Rev. 27, 39 (1940). 

2 G. B. Kistiakowsky and W. W. Rice, J. Chem. Phys. 8, 610 (1940). 


3 E. B. Wilson and A. J. Wells, J. Chem. Phys. 9, 319 (1941). 
4 J. G. Aston, private communication. 





The Exchange Reaction between Gaseous 
and Combined Nitrogen 


Yosuto NisHina, TAKEO Irmori, Hipeo Kuso 
AND HiromMt NAKAYAMA 


The Institute of Physical and Chemical Research, Komagome, 
Hongo, Tokyo, Japan 


June 2, 1941 


Y using heavy nitrogen gas Burris and Miller' have 

recently concluded that no exchange takes place 
between gaseous and combined nitrogen during the process 
of nitrogen fixation by Azotobacter. A possibility, however, 
cannot be excluded that in some other cases the exchange 
reaction of nitrogen may happen between nitrogen gas and 
some nitrogenous constituents of the plants. In fact, the 
following results of our experiments show quite clearly 
that the exchange reactions of nitrogen do occur to a 
certain extent between nitrogen gas and some nitrogen 
compounds, which are considered to exist in plant juices. 
Whether or not these exchanges are strong enough to be 
taken into account for the experiments of nitrogen fixation 
will be seen in the future. 

Radioactive. nitrogen N™® gas (half-life 10.5 min.) was 
obtained by bombarding amorphous carbon with 3-Mev 
deuterons from our cyclotron, the carbon being then heated 
in a closed vessel containing air. In order to get the gas 
completely freed from impurities NO, NO2 and NHs, it 
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TABLE I. 








RELATIVE 
INTENSITY 
OF RADIO- 
ACTIVITY 


CONCEN- 
TRATION SUBSTANCE 
SUBSTANCE TESTED MoLe/! PRECIPITATED 





NaNO: 1/10, 1/20 _nitron nitrate 

NaNOs nitron nitrate 

NH:OH -HCl nitron nitrate 

NH.Cl (NH4)2PtCle, Hgx-ON Hel 
KCN AgCN 

NeHa-H2SO, 


benzalazine 
NasCo(NO2)6 K2NaCo(NOs2)e 
8-hydroxyquinoline 
acetate TiO-hydroxyquinolate 
CO(NH2): (evaporated up) 








was vigorously shaken repeatedly with caustic soda solu- 
tion. containing potassium permanganate, and with dilute 
sulphuric acid. For the study of the exchange reaction 
mentioned above, the gas thus purified was introduced into 
a bottle containing aqueous solutions of various nitrogen 
compounds in question, and shaken for a few minutes. 
After the gas was carefully driven out of the solutions, the 
bottle was thoroughly evacuated in order to get rid of 
nitrogen gas absorbed by the solutions. Nitrogen com- 
pounds under consideration were then precipitated by the 
addition of appropriate reagents and the radioactivities of 
the precipitates were examined by a Geiger-Miiller tube 
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counter. Nitrite and hydroxylamine, before being pre- 
cipitated, were oxidized to nitrate by adding excess of 
acidic permanganate and hypochlorite, respectively. The 
results are shown in Table I. 

As is seen from the results, the exchange reaction of 
nitrogen occurs between gaseous nitrogen and _ nitrite, 
hydroxylamine, etc. The rate of exchange reaction in the 
case of nitrite was rather rapid, though not immeasurable. 
The concentration of the solution was proved to have a 
remarkable effect on the rate of the reaction. In the cases 
of nitrite or nitrate, for example, this rate was constant in 
dilute solutions, but decreased with the increase of concen- 
tration above 1/10 mole/I and no exchange was observed 
even in the solution of 2 moles/I. As it is difficult to con- 
ceive that the amount of nitrogen gas absorbed by such 
dilute solutions would differ according to compounds in 
such a marked degree, the effects observed are hardly to 
be ascribed to the phenomena of gas absorption by solu- 
tions. Some of the solutions, moreover, are all converted 
into nitrate before precipitation and the great difference 
still observed among these compounds confirms our 
conclusion. 

The absolute rate of reaction together with the details 
of experiment will be published in the Scientific Papers of 
this Institute. 


1R. H. Burris and C. E. Miller, Science 93, 114 (1941). 
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